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Foreword 

The first Workshop on Exotic Physics with Neutrino Telescopes was held in Up- 
psala, Sweden, between September 20-22, 2006. The main idea behind the organi- 
zation of the workshop was to discuss the possibilities of using large-scale neutrino 
telescopes to study physics topics beyond the search for astrophysical point sources 
of neutrinos, which constitutes the main rationale for the construction of such de- 
tectors. The label 'exotics' is therefore to be understood in this spirit. Indeed many 
of the topics discussed in the workshop were not 'exotic' at all in a generic sense. 

The workshop was based on several review talks giving an overview of relevant 
subjects, and contributed talks from experiments and theorists touching upon more 
detailed developments. 

The smaller neutrino detectors, MACRO, Baksan and Super-K and the large- 
scale neutrino telescopes already in operation, Baikal and AMANDA, reviewed the 
status of their analyses and results on topics ranging from particle searches to tests 
of fundamental laws, proving that their detection techniques can produce competi- 
tive limits in comparison with present accelerator results. The next-generation ex- 
periments, IceCube at the South Pole and ANTARES and NESTOR in the Mediter- 
ranean (with their km'^ successor, KM3NET, already under R&D), will be able to 
improve the current limits by a few orders of magnitude, although we all expect 
that such detectors will lead to a discovery rather than limit improvements. 

The theory talks provided guidelines for experimental searches, and covered the 
signatures expected in neutrino telescopes from different dark matter candidates, 
micro black holes, models with extra dimensions, non-standard neutrino oscillation 
scenarios and new neutrino interactions. A wealth of possibilities for neutrino de- 
tectors to search for. 

The view of the 'exotic' universe is a comprehensive one, with results from satel- 
lite, accelerator and laboratory experiments complementing each other. The pro- 
gram included reviews from accelerator, underground labs and space-based projects 
whose research topics overlap with those of neutrino telescopes. 

The workshop was held in the Polhem auditorium of the Angstrom Laboratory, 
the complex housing all physics departments of Uppsala University, and brought 
together about 50 scientists from 16 countries. 

The support from the Swedish Research Council (Vetenskapsradet) and the Fac- 
ulty of Physics of Uppsala University are kindly acknowledged. 

Carlos P. de los Heros 
Editor 

Uppsala, October 2006. 
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Abstract 

The main goal of this article is to give an update on the experimental status of direct 
dark matter searches. Before doing that the evidence for dark matter will be very briefly 
reviewed along with a few words about the candidate solutions to the dark matter prob- 
lem. The experimental requirements for detection of the most favoured candidate, the 
neutralino, will be given, which then leads into a review of the current status and future 
plans for experiments across the world. 



1. The Evidence for [Cold] Dark Matter 

The evidence for the existence of dark matter is now extensive and present on all 
scale-lengths from galaxies up to the Universe as a whole. The direct observational 
evidence has steadily accumulated since the first observations by Oort and Zwicky 
in the 1930s and includes galaxy rotation curves, galaxy cluster velocity dispersion, 
cluster x-ray gas, gravitational lensing, and cluster streaming motions. In addition 
we now have a powerful suite of information being used to produce models of the 
Universe as a whole, which start with primordial density perturbations, as seen from 
microwave anisotropics measured by COBE and WMAP, and then look to see how 
these evolve under the action of gravity to form the structures we see today. Part 
of this process requires a knowledge of the structure and evolution of space-time 
itself, which is informed by high-redshift populations of standard tracers, such as 
supernova type la. The evolution process itself is then modeled using sophisticated 
N-body simulations and the resulting predicted galaxy distributions are then con- 
fronted with real data from large survey projects. Typical models of the Universe 
use a dozen or so parameters, amongst which arc the total matter density, ^m, the 
baryonic matter density, fif,, and the hot dark matter density, The quality of the 
fits between predictions and data, such as the microwave background anisotropy an- 
gular power spectrum and the matter density power spectrum (as seen from galaxy 
distributions) is now extremely good and they provide strong constraints on com- 
binations of parameters used in the modeling. For example Tegmark et al. [1] use 
microwave anisotropy data from WMAP and large scale structure data from SDSS 
to deduce that the total mass density in the Universe is Q.m = 0.24 ± 0.02. The 
baryonic mass density, which contains all the normal matter, can only contribute 
ilfc = 0.042 ± 0.002. Whilst this is nicely consistent with big-bang nucleosynthesis 
models, it clearly leaves the need for some other dominant mass component at the 
several sigma level. This 'dark matter' is usually divided into two classes, hot and 
cold, depending on whether it was relativistic in the earl Universe. Neutrinos are a 
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natural candidate for hot dark matter and the modeling gives fl^ < 0.024, which 
then leaves most of the unseen dark matter in some non-baryonic form which is not 
neutrinos or the like; i.e 'cold dark matter'. 

2. The Galactic Candidates 

N-body simulations of structure growth in the Universe show that galaxies them- 
selves are also cold dark matter dominated which produces a consistent solution for 
the ubiquitous flat rotation curves. In looking for candidate explanations for the 
cold dark matter in galaxies, and hence in our Milky Way, it is natural to first focus 
on possibilities which are already motivated for other reasons and which can provide 
dark matter both for the Universe as a whole and for galaxies. Two prime candidates 
exist. These are new particles called the axion and the neutralino. The first of these 
was postulated in the 1970s to suppress strong CP violation, which should otherwise 
have been much prevalent. Recently [2] there has been much excitement caused by 
a positive claim from a laboratory search for axion-like interactions. Unfortunately 
the inferred particle properties are many orders of magnitude away from what was 
expected for cosmologically relevant axions and additional observations and theo- 
retical development are needed to consolidate the result and its relevance to dark 
matter. The second well motivated particle candidate is the so-called neutralino, 
which is the lightest new particle coming out of supersymmetry. The properties of 
supersymmetric weakly interacting massive particles (WIMPs) can be such that 
they were produced in the early Universe and that the lightest of these remains 
with us today in cosmologically important quantities. In addition it turns out that 
their weak scattering cross-sections make them potentially detectable in laboratory 
based direct searches. 

Before discussing neutralino detection in detail below it is useful to also mention 
the possibility of resolving the dark matter problem not be adding new mass but by 
modifying the way in which gravity and/or the inertial response to gravity works. 
This was first suggested as Modified Newtonian Dynamics, MOND, by Milgrom [3] 
and, as a phenomenological approach, works very well [4]. Bekenstein [5] has since 
produced a more formal theoretical basis for MOND. including the relativistic sec- 
tor. This formulation received a tough test recently through the observation of two 
galaxy cluster which have passed through each other causing the main baryonic 
matter component, hot cluster gas, to separate from the cold dark matter compo- 
nent as seen through its lensing effects [6]. Paramos and Bertolami [7] summarise 
the current attempts to recover this behaviour from within MOND. 

3. Neutralino Detection 

The dominant interaction between neutralinos, X: and any target medium will be 
its elastic scattering from a nucleus. The energy imparted to the recoiling nucleus 
is then the means of direct detection of these particles. Indirect detection through 
measuring the annihilation products (neutrinos, 7-rays and other particles) from 
XX '^ill be covered elsewhere. 

When a neutralino scatters elastically from a target nucleus it recoils with a ki- 
netic energy typically below ^ 50 keV. As shown in figure 1 this energy deposit 
can result in three measurable effects; photon emission, charge production and/or 
phonon generation. Experiments can attempt to measure any of these and many 
experiments try to measure combinations of them simultaneously to obtain better 
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WIMP elastic nuclear 
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of energy at a rate 1 0-^ to 
1 event/day/kg 
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range of complementary 
techniques 



Fig. 1. Illustrating the possible techniques and diversity of direct searches 



diagnostic information of the cause of each event as the relative efficiency of chan- 
nehng energy into the three signal channels depends on dE/dx and hence particle 
species. In the following each technique is briefly reviewed in turn before presenting 
a summary of current status. 



3.1. Scintillation detectors 



In the class of scintillation detectors the experiments compare the shape of the 
background energy spectrum to that expected from WIMP nuclear recoils. In ad- 
dition they often make use of other information contained within the scintillation 
signal, such as an annual modulation in the rate/spectrum or pulse shape discrim- 
ination. The DAMA experiment used ^ 100 kg of Nal and searched for. and found, 
an annual modulation of the expected form at 6.3(t [8]. No supporting evidence has 
yet been seen by any other experiment and the same team has a new experiment, 
LIBRA, with a 300 kg array [9]. Similar experiments with Nal include ANAIS and 
NAIAD. ANAIS [10] is in development as a new 107 kg experiment at CanFranc. 
NAIAD [11,12] is a completed experiment run at Boulby which used pulse shape 
discrimination but not annual modulation. ELEGANT V [13] also used a large 
mass Nal searching for annual modulation but did not achieve sufficient sensitivity 
to check the DAMA result. They also tried looking for a spectral line from an in- 
elastic scattering process. Another crystal scintillator being developed into a high 
mass instrument is Csl. The KIMS experiment has just completed a pilot run with 
a 6.6 kg detector giving an upper limit which just reaches the top of the DAMA 
positive detection area [14]. 
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Another class of scintillator is the liquid noble gases. These exhibit particle dis- 
crimination through the scintillation time-constant. The light emission comes from 
a radiative dissociation of an excited dimer. However the dimer can be in either a 
singlet or triplet excited state and the decay lifetimes in these two states differs. The 
relative populations ending up in the singlet and triplet states depends on dE/dx. 
This time constant discrimination in xenon has already been used by the ZEPLIN 
I experiment [15]. Argon and neon can also be used in this mode, and indeed show 
somewhat greater time constant ratios. DEAP(argon) [16] and CLEAN(neon) [17] 
are two instruments in advanced development. Liquid noble gases offer excellent 
prospects for scaling to much larger mass detectors. XMASS is intending to build 
an 800 kg instrument relying on self-shielding to provide an inner background-free 
volume, thus no longer needing particle discrimination. Currently the project has 
results from a 100 kg prototype [18] 

3.2. Hybrid scintillation /ionisation detectors 

All the experiments shown as scintillation/ionisation hybrids in figure 1 are based 
on the noble gases. In fact most use two-phase systems in which the main target 
mass is in liquid phase. Prompt scintillation is created at the point of interaction 
in the liquid. However any ionisation released is drifted through the liquid phase 
towards the surface by an applied electric field. At the surface the field is strong 
enough for the electrons to be extracted into the gas phase. Here they cither cause 
proportional electroluminescence giving another burst of light or they can be di- 
rectly measured as charge. The argon instruments are ArDM [19] and WARP [20]. 
Both of these build on the heritage from ICARUS and aim to develop tonne-scale 
targets. ArDM is building a 600 kg prototype using many innovative design features 
and using an existing containment vessel. WARP has built a 2.3 litre prototype and 
has aheady some encouraging first results. XENONIO [21], ZEPLIN II [22] and 
ZEPLIN III [23] are xenon based instruments. XENONIO has a target mass of 
~ 10 kg and is operating in Gran Sasso. ZEPLIN III has a target mass of ~ 40kg 
and is operating in Boulby. XENONIO has a somewhat better light collection for 
the scintillation signal as it has photomultiplers immersed in the liquid phase as 
well as in the gas phase. Both use bulk PTFE as a structural material with insula- 
tion and with good reflectivity at the UV xenon wavelength. Results are expected 
from both experiments imminently. ZEPLIN HI has completed its build and test 
phase [24] and is now ready for deployment in the Boulby laboratory. It has several 
novel engineering and physics design features which probably give it a better grasp 
than either XENONIO or ZEPLIN II. 

Finally a relative newcomer is SIGN [25] which is proposing to use high pressure 
neon gas in a conventional cylindrical wire chamber configuration. Csl coatings are 
to be used to convert the scintillation signal into charge. 

3.3. Hybrid cryogenic detectors 

The hybrid cryogenic detectors combine a bolometric signal channel with either 
ionisation or scintillation. The technology demands small segmented detectors and 
this is currently limiting the total target masses to modest values (< 10kg). 
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CDMS [26] and EDELWEISS [27] use ionisation detectors with thermal sensing 
elements bonded to them. CDMS uses Ge and Si, whereas EDELWEISS only uses 
Ge. CDMS is midway through a second generation experiment, CDMS II, and 
has published results from mid-term operations. Plans are well developed for the 
next phase SuperCDMS [28] to increase the target mass to 25 kg. EDELWEISS has 
completed a first phase experiment and is now readying its second generation for 
operation with a partial complement of detectors. 

CRESST II [29] and ROSEBUD [30] use scintillators with thermal sensing ele- 
ments. The scintillation light is proximity coupled into photodetectors. CRESST 
has operated a 300 g CaW04 crystal in Gran Sasso and has established useful lim- 
its from it. They now have a second generation instrument allowing up to 10 kg 
of target crystal waiting for full deployment. ROSEBUD has demonstrated a 46g 
prototype using a BGO crystal. 

3.4. Bubble chambers 

Figure 1 shows a few experiments in the phonon only segment. Two of these are 
worthy of discussion here as they are progressing rapidly. These are COUPP [31] 
and PICASSO [32], which are bubble chambers. Contrary to the bolometers they 
not specifically sensitive to the total energy deposit but rather to the local dE/dx 
in the track of the recoiling nucleus. 

COUPP is a heavy-liquid bubble chamber and which uses CF3I allowing a good 
sensitivity to both spin-dependent and spin-independent scattering cross-sections. 
A 2 kg prototype is currently operating underground at 300 m water equivalent at 
Fermilab. 

PICASSO use superheated droplets of C4F10 and looks for bubble nucleation 
caused by elastic WIMP scattering from the F nuclei. This makes the experiment 
particularly effective for searching for spin-dependent scattering. The experiment is 
operating at 6000m water equivalent depth in SNOLAB. An earlier stage has been 
completed with some 19 g of F. The current phase is scaling up to 2 kg active mass 
of F. 

3.5. Directional gas time projection chambers 

A critical diagnostic to bring to play once a detection appears secure is to study 
the directional signature on the WIMP velocity distribution imparted by the motion 
of the solar system as it orbits around the Galactic centre. This will be modulated 
by both the orbit of the Earth around the Sun and, in detector coordinates, but 
the Earth's rotation. These will be very powerful signatures to confirm the Galactic 
origin of the signals and, ultimately to probe the local WIMP velocity distribution 
as an astrophysical tool. To be able to determine a WIMP direction requires the 
nuclear recoil track to be extended allowing the track direction and sense to be 
recovered. This points to low-pressure gas detectors with very good position sensi- 
tivity and fine energy resolution to allow the track to be sampled along its length 
for head-tail determination. Two projects are pursuing this goal; DRIFT [33] and 
NEWAGE [34]. DRIFT uses Im^ modules containing CS2 gas at 40torr pressure. 
It is now in second generation prototype implementation operating in the Boulby 
mine at 3000 m water equivalent depth. The NEWAGE project is developing a 
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micro-TPC (10 x 10 x 10 cm'^) using CF4 gas at ISOtorr. Current tests arc aimed 
at characterizing the response of this and Ar-C2H6 to nuclear recoils. 



4. Current Status 

Figure 3.4 summarizes the current status of direct dark matter searches in the 
context of a standard model of the Galaxy, which fixes the local dark matter density 
and velocity distribution and allows each experiment to be assessed in terms of 
the spin-independent WIMP scattering cross-section normalized to an individual 
nucleon. The solid curves correspond to the experimental upper limits. Drawing 
a vertical line at a mass of 100 GeV the lowest line is that from CDMS II. Going 
upwards the next curve reached is from ZEPLIN I, then EDELWEISS and CRESST 
II. These last three curves all pass through the DAMA positive detection allowed 
region, with ZEPLIN I just making a marginal contact. CDMS II is well below the 
DAMA region and hence excludes it in the context of the model assumptions. The 
DAMA region is shown as a filled area, just above a cross-section of lO""'^ cm^ and 
between masses of ^ 40 to ~ 80 GeV. Other filled areas are regions populated by 
the many different supersymmetric models which are possible. These predictions 
span many orders of magnitude in cross-section and the dashed lines then indicate 
the extra reach the experiments will have as they progress through 100 kg to 1000 kg 
target masses. 

Figure 3.4 gives more hope to the experimentalists by concentrating on a con- 
strained minimal supersymmetry model (CMSSM) parameter space [35]. Using a 
Markov Chain Monte Carlo simulation the authors have assigned a probability rat- 
ing to each region of parameter space. Although this does not rule out the need for 
tonne-scale targets to reach down to the 10"^" pb level it does indicate that future 
improvements in experiments will rapidly start to probe new space with significant 
probabilities of a positive detection and a number of such experiments are indeed 
poised to make inroads into the space. 
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Summary of search for exotic phenomena with 
Super-Kamiokande 

Jun Kameda*^, 

'^Institute of Cosmic Ray Research, University of Tokyo, Japan 



Abstract 

In this paper, I present the results on exotic phenomena search in Super-Kamiokande. 



1. Introduction 

The Super-Kamiokande detector is a cylindrically-shaped water Cherenkov de- 
tector with 50 kiloton of ultra-pure water. It is located about 1000m underground in 
the Kamioka Observatory in the Kamioka mine in Gifu Prefecture, Japan. Super- 
Kamiokande is a multi purpose experiment: Neutrino physics using atmospheric 
neutrino and solar neutrinos, search for astrophysical neutrinos, nucleon decay 
search, exotic particle search (Q-ball, WIMPs, magnetic monopoles, etc.), and so 
on. In this paper I report the results of indirect WIMPs search and new physics 
search using atmospheric neutrinos (CPT violation, Lorentz invariant violation). 

The observed atmospheric neutrino events in Super-Kamiokande can be cate- 
gorized into four types:(l) All visible particles are contained in the detector, we 
call these Fully-Contained (FC) events, (2) at least one visible particle escaped 
from the detector, called Partially-Contained (PC) events, (3) The neutrino in- 
duced muons entering and passing through the detector from below, we call these 
upward-through going muons (UT) (4) The neutrino induced muons entering and 
stopping the detector from below, called upward stopping muons (US). Details of 
the event selections and reconstructions are described in [1]. 

2. Indirect WIMPs search using neutrino-induced muons 

Weak Interacting Massive Particles (WIMPs) are one of the plausible candidates 
of dark matter. The lightest supersymmetric particle (LSP) of supersysmmetric 
theories is the most theoretically well developed WIMP candidate. Many super- 
symmetric theories predict LSPs with mass from several 10s GeV to TeV region. 
The WIMPs orbiting in the galactic halo will be captured by celestial objects like 
the Sun and the Earth, and will annihilate in pairs to the secondaries (cc, 66, rr,..). 
High energy neutrinos will be produced as a decay product of the secondary par- 
ticles. A good channel to observe such high energy neutrinos is a entering muons 
induced by the neutrino interactions in the rock surrounding the detector, because 
of their large target volume. The backgrounds are the produced high energy muons 
induced by the atmospheric neutrinos. 

Figure 1 shows the zenith angle of the 1679.6 live days of the upward-going 
muons. No significant excess was found in the signal region. 
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Fig. 1. Zenith angle distribution of upward 
tlirougli-going muons with respect to the cen- 
ter of the Earth. The circles show data, Boxes 
and solid line show the MC prediction with- 
out and with neutrino oscillation, respectively. 
(Am2,sin2 2e)=(2 X 10"^, 1.0) was assumed. 



Fig. 2. Supcr-K WIMP-induced upward-muon 
flux limits from the Earth as a function of 
WIMP mass along with those from the other 
experiments [9—11]. 



The signal window size whieh ineludes 90% of the neutrino-induced muons is 
a function of neutrino energy, because of the finite scattering angle of neutrinos, 
muon multiple scattering, and the finite size of the WIMP annihilation region in 
the Earth. The neutrino energy from WIMPs annihilation depends on the WIMPS 
masses, so the flux limit depends on the WIMPs masses. We employed the result 
of the Monte Carlo simulation done in Ref. [3] to calculate the angular windows 
which contain 90% of the signal for various WIMPs masses. Figs 2-4 show the flux 
limits as a function of the WIMP mass from the Sun, the Earth, and the Galactic 
center, respectively. 

The limit of the indirect WIMPs search can be translated and compared to the 
results from the direct search experiments using the calculated ratio of the direct- 
to indirect detection rates [4]. The calculated ratio has uncertainties in the model 
dependence of the WIMP annihilation, and we used the maximum value of the ratio 
to get conservative limits. 

Figs. 5 and 6 show the 90% excluded regions in WIMP parameter space from 
Super-Kamiokande and the results of the other experiments for scalar coupling and 
axial vector coupling, respectively. The sum of the flux from the Sun and the Earth 
has been used for scalar coupling, and only the flux from the Sun was used for axial 
vector coupling. We had excluded a significant portion of the favored region from 
the DAMA experiment in scalar coupling, and largely pushed the limit on axial 
vector coupling. 



3. Search for new physics on neutrino using atmospheric neutrinos 

The observed deficit pattern of the atmospheric neutrino is perfectly explained by 
the —^ i>T neutrino oscillation due to neutrino masses [1]. But the finite masses 
of the neutrino is not the only possible source of neutrino oscillation. In general, 
survival probability of the neutrino oscillation can be written as: 
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Fig. 3. Super-K WIMP-induccd upward- muon 
flux limits from the Sun as a function of WIMP 
mass along with those from the other experi- 
ments [9,10,12]. 



Fig. 4. Supcr-K WIMP-induccd upward-muon 
flux limits from the Galactic center as a func- 
tion of WIMP mass along with those from the 
other experiments [9] 
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Fig. 5. Super-K 90% C.L. excluded region 
in WIMP parameter space (solid line) for 
WIMP with scalar coupling. DA-MA 3fT al- 
lowed region(filled) [13], 90% C.L. excluded 
region from CDMS(dot-dashed) [14], EDEL- 
WEISS(dashed)[15], and ZEPLIN(dotted)[16] 
arc also shown. 



Fig. 6. Super-K 90% C.L. excluded region in 
WIMP parameter space (solid line) for WIMP 
with axial vector coupling. The 90% C.L. ex- 
cluded region from UKDMC(dashed)[17], EL- 
EGANT(dot-dashed)[18] are also shown. 



P{y^ -^v,,)^\~ sin' 26 sin' [P-^)' (1) 

where is a mixing angle of neutrinos. (3 is the frequency of the neutrino oscillation, 
and the power index ri is a parameter which depends on each theory. Examples of 
these theories are: violation of Lorentz invariance [5] (n = -1), violation of the 
equivalence principle [6](n = -1 ), CPT violation [7](n = 0), and coupling to space- 
time torsion fields [8](n = 0). In the case of 'standard oscillation', the index n is 1. 
We carried test in three dimensional parameter space (sin^ 29,(3,n) to find the 
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Fig. 8. The 68%,90%,99% C.L. allowed region 
on {Aa,Am'^) space. 



best value of n. 1144 days FC + PC, 1138 days of UT, and 1117 days of US were 
used for this test. Fig. 7 shows the distribution as a function of n. The is 
minimized by optimizing the sin^ 29, [i and the systematic parameters. We fitted 
a parabola function near the minimum region, and obtained n — 1.14 ± 0.11, 
which is consistent with the neutrino oscillation induced by finite neutrino masses 
(n = 1.0). 

The new physics cannot be a dominant source of the neutrino oscillation, but 
there are possibility of the sub-dominant source of the neutrino deficit. We tested 
two scenarios as a sub-dominant mechanism. The survival probability of the neu- 
trino with neutrino masses and CPT violating term can be written as: 



1 



sin^ 29 sin^ 



Aa L 



(2) 



where Aa is the CPT violation parameter. Fig. 8 shows the allowed region of the 
parameter space using 1489.2 live days FC-I-PC, 1645.9 live days of upward-going 
muons. The definition of the is shown in [1]. The result is consistent with no 
CPT violating term, and we got the limit of the Aa < 10"^^ at 90% C.L. 

A type of Lorentz violating term can also introduce neutrino oscillations. The 
survival probability of the neutrino with neutrino masses and a CPT violating 
term can be written: 



P{Vf, sin^ 29 sin^ (n) 



(3) 



where tan^ 6 = (l + {E/Ecf sin26l„) / {{E/Ecf cos26'„), Ec = y/Am'^/2c^T^ ^ ^ 
1.27^(Am2L/£:)2 ± 4c^^ sin26'„L£; + ^{c^'^LEf, c^^ is the difference of the at- 
tainable velocity of neutrinos, 9^ is a mixing angle of the velocity eigenstate, respec- 
tively. The plus/minus sign in 17 corresponds to the relative phase of the mixing of 
velocity eigenstate and mass eigenstate, and we consider both signs. We scanned 
the (c'^'^,sin 29y) space assuming Am? = 2.4 x lO^'^eV^, and full mixing of the mass 
eigenstate, which is the best fit value from the neutrino oscillation due to neutrino 
masses. Figs. 9 and 10 show the allowed regions of the parameter space for both 
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Fig. 10. 68%, 90%, and 99% allowed region of 
the parameter space for the sign in Q is minus 
case. 



sign. The result is consistent with no Lorentz invariant violating term, and we got 
the hmit of c^^ < 10^24 99% c.L. for both cases. 
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Abstract 

MACRO was a multi-purpose experiment that took data from 1989 to 2000, at the under- 
ground Laboratory of Gran Sasso (Italy). MACRO gave important results/limits concern- 
ing: (i) the oscillation of atmospheric neutrinos, also in the non-conventional scenario of 
violations of Lorentz invariance, (ii) the searches for exotic particles (supermassive GUT 
magnetic monopoles, nuclearites, WIMPs), (iii) muon physics and astrophysics. A sum- 
mary of the MACRO results will be presented and discussed, focusing the attention on 
the exotica searches. 



1. Introduction 

MACRO was a large area multipurpose underground detector [1] designed to 
search for rare events and rare phenomena in the cosmic radiation. The detector 
was located in Hall B of the underground Gran Sasso Laboratory (Italy). It was 
optimised to search for the supermassive magnetic monopoles [2,3] predicted by 
Grand Unified Theories (GUT). The experiment obtained important results on 
atmospheric neutrino oscillations [4-8] and performed neutrino astronomy studies 
[9], indirect searches for WIMPs [10], search for low energy stellar gravitational 
collapse neutrinos [11], studies of the high energy underground muon flux (which 
is an indirect tool to study the primary cosmic ray composition and high energy 
hadronic interactions [12]), searches for fractionally charged particles (LIPs) [13] 
and other rare particles that may exist in the cosmic radiation. 

The detector started data taking in 1989 and it was running until December 2000. 
The apparatus had global dimensions of 76.6 x 12 x 9.3 m'^ and was composed of 
three sub-detectors: liquid scintillation counters, limited streamer tubes and nuclear 
track detectors. Each one of them could be used in "stand-alone" and in "combined" 
mode. It may be worth to stress that all the physics and astrophysics items listed 
in the 1984 proposal were covered and good results were obtained on each of them. 

2. Atmospheric neutrino oscillations 

MACRO detected i^^-induced muon events in 4 different topologies. 

The upthroughgoing muons come from interactions in the rock below the 
detector; the f^'s have a median energy of ^ 50 GcV. 

Fig. 1 shows the zenith distribution of the measured 902 upthroughgoing muons 
(black circles) compared with two MontcCarlo (MC) predictions: the Bartol96 [14] 
flux with and without oscillations (the dashed and solid lines, respectively) and 
the Honda2001 flux [15]. The FLUKA MC predictions [16] agree perfectly with the 
Honda2001. 
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Fig. 1. Comparison of the MACRO up- 
ward-throughgoing muons (black cir- 
cles) with the predictions of the Bar- 
tol96 and of the Honda2001 MC oscil- 
lated and non oscillated fluxes (oscilla- 
tion parameters Am? = 2.3 ■ 10"'^ eV-^ 
and sin^ 29m = 1). 
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Fig. 2. Ratio Data/MCno osc as a function of the es- 
timated L/Ev for the upthroughgoing muon sample 
(black points). The solid line is the MC expectation 
assuming Am^ = 2.3 ■ 10"'' eV-^ and sin^ 29m = 1. 
The last point (black square) is obtained from the lU 
sample. 



For a subsamplc of ~ 300 upthroughgoing events, wc estimated the muon energy 
through Multiple Coulomb Scattering in the rock absorbers in the lower apparatus 
[7]. The evaluated resolution on is ^ 100%. The parent neutrino path length 
is L ^ 2Re cosQ, where Re is the Earth radius. Fig. 2 shows the ratio data/MC 
as a function of the estimated L/Ey for the upthroughgoing muon sample. The 
black circles are data/Bartol96 MC (assuming no oscillations); the solid line is 
the oscillated MC prediction for Am^ = 2.3 • 10"^ eV^ and sin^ 20„j = 1. The 
shaded region represents the simulation uncertainties. The last point (black square) 
is obtained from semicontained upward going muons. 

The low energy events (/[/, ID + UGS [5]) are produced by parent t'^; interacting 
inside the lower detector, or by upgoing muons stopping in the detector. The median 
energy of the parent neutrino is ~ 3 — 4 GeV for all topologies. In both cases, 
the zenith distributions are in agreement with the oscillation prediction with the 
optimised parameters [8]. 

In order to reduce the effects of systematic uncertainties in the MC absolute 
fluxes we used the following three independent ratios [8] : 

(i) High Energy data: vertical/horizontal ratio, i?i = Ny^rt/Nhor 

(ii) High Energy data: low energy/high energy ratio, i?2 = Nio^/Nhigh 

(iii) Low Energy data: i?3 = {Data/MC)iu/iData/MC)iD+UGS 
Combining the three independent results, the no oscillation hypothesis is ruled out 
at the ~ 5a level. 

To evaluate the hypothesis of oscillation for different values of Am^ and sin^ 20,„, 
the Feldman- Cousins [17] procedure was used and the corresponding 90% C.L. 



region for the 



Am2 = 2.3 • 10-3 eV2 and sin^ 29,^ = 1. 



oscillation is given in ref [8]. The best fit is reached at 



2.1. Search for exotic contributions to atmospheric neutrino oscillations 

MACRO searched for "exotic" contributions to standard mass-induced atmo- 
spheric neutrino oscillations, arising from a possible violation of Lorentz invariance 
(VLI), using two different and complementary analyses. The first approach uses the 
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low energy {Ei^ < 28 GeV) and the high energy {E^ > 142 GeV) samples. The mass 
neutrino oscillation parameters have the values given in Seet. 2 and we mapped the 
evolution of the estimator in the plane of the VLI parameters Af and sin^ 26y . 
No improvement was found, so we applied the Feldman-Cousins [17] method to 
determine 90% C.L. limits on the parameter: \ Av\ < 3 • 10~^^ [18]. 

The seeond approach exploits a data subsample characterised by intermediate 
neutrino energies. It is based on the maximum likelihood technique and considers 
the mass neutrino oscillation parameters varying in the 90% C.L. border [8]. The 
obtained 90% C.L. limit on the Av parameter is also around 10^^^ [19] . 

3. Neutrinos from astrophysical sources 

3.1. Search for astrophysical HE muon neutrinos 

High energy i/^'s are expected to come from several galactic and extra-galactic 
sources. An excess of events was searched for around the positions of known sources 
in 3° (half width) angular bins. The 90% C.L. upper limits on the muon fluxes from 
specific celestial sources were in the range 10~^^ 10"^'' cm~^ s~^ [20]. A search 
for time coincidences of the upgoing muons with 7-ray bursts was also made. No 
statistically significant time correlation was found [9] . 

A different analysis was made for the search for a diffuse astrophysical neutrino 
flux, using a dedicated method to select higher energy upthroughgoing muons. The 
flux upper limit was set at the level of 1.5 • 10"^'^ cm~^ s~^ [21]. 

3.2. Indirect searches for WIMPs 

Weakly Interacting Massive Particles (WIMPs) could be part of the galactic dark 
matter; they could be intercepted by celestial bodies, slowed down and trapped in 
their centres, where WIMPs and anti- WIMPs could annihilate and yield neutrinos 
of GeV or TeV energy, in small angular windows from their centres. One WIMP 
candidate is the lowest mass neutralino. 

To look for a WIMP signal, we searched for upthroughgoing muons from the 
Earth centre, using 10° 15° cones around the Nadir; the 90% C.L. muon flux 
limits are 0.8 1.4 • 10^^"* cm^^ s~^ [10]. These limits, when compared with the 
predictions of a supersymmetric model, eliminate a sizable range of parameters used 
in the model. 

A similar procedure was used to search for z/^ from the Sun: the muon upper 
limits are at the level of about 1.5 2 • lO"^"' cm~^ [10]. 

3.3. Neutrinos from stellar gravitational collapses 

A stellar gravitational collapse of the core of a massive star is expected to produce 
a large burst of all types of neutrinos and antineutrinos with energies of 5 ~ 60 MeV 
and with a duration of ~ 10 s. No stellar gravitational collapses in our Galaxy were 
observed from 1989 to 2000 [11]. 

4. Search for exotic particles 

4.1. Search for GUT magnetic monopoles (MMs) 

Supermassive magnetic monopoles predicted by Grand Unified Theories (GUT) 
of the electroweak and strong interactions should have masses tum ^ 10^^ GeV. 
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Fig. 3. 90% C.L. upper limit obtained by 
MACRO for an isotropic flux of GUT MMs 
with g = go compared with direct hmits given 
by other experiments. 



Fig. 4. 90% C.L. upper hmits versus mass for 
downgoing nuclearites with /3 = 2 ■ 10~^ at 
ground level. The MACRO limit for nuclearite 
masses larger than 5 ■ lO'^^ GeV/c^ has been ex- 
tended and corresponds to an isotropic flux. 



MACRO was optimised to search for an isotropic flux of GUT MMs in the cosmic 
radiation. The three sub-detectors had sensitivities in different P regions, covering 
the velocity range 4 • 10~^ < /3 < 1. They allowed multiple signatures of the same 
rare event candidate. No candidates were found by any of the three subdetectors. 
Fig. 3 shows the 90% C.L. flux upper limits for g = go poles (one unit of Dirac 
magnetic charge) plotted versus (3 [2] together with direct limits set by other ex- 
periments [22]. The MACRO MM direct limits are by far the best existing over a 
very wide range of /?. 

The interaction of the GUT monopole core with a nucleon can lead to a reaction 
in which the nucleon decays, M +p 71/ + e+ +7r". MACRO dedicated an analysis 
procedure to detect nucleon decays induced by the passage of a GUT MM in the 
streamer tube system (a fast e"*" track from a slow (/3 ~ 10~^) MM track). The 
90% C.L. flux upper limits established by this search are at the level of 3 • 10~^^ 
cm~^ s~^ sr^^ for 10~* < /? < 0.5 • 10~^; they are valid for catalysis cross sections 
5- 102 < a^at < 10^ mb [3]. 

4.2. Search for nuclearites, Q-balls and LIPs 

Strangelets should consist of aggregates of u, d and s quarks in almost equal 
proportion [23] and would have typical galactic velocities (3 ~ 10^'^. The MACRO 
90% C.L. upper limits for an isotropic flux of nuclearites with 10~^ < /3 < 1 was at 
the level of 1.5 • lO'^^^ cm'^ s'^ sr"! [24]. 

MACRO searched also for charged Q-balls (aggregates of squarks, sleptons and 
Higgs fields) [25], giving an upper limit of ~ 10~^^ cm"^ s~^ sr"-'^ [26]. 

Fractionally charged particles could be expected in Grand Unified Theories as 
deconfined quarks; the expected charges range from Q = e/5 to Q = 2/3e. LIPs 
should release a fraction (Q/e)^ of the energy deposited by a muon traversing a 
medium. The 90% C.L. flux upper limits for LIPs with charges e/3, 2/3e and e/5 
are at the level of lO^^^ cm^^ sr^^ [13]. 
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5. Conclusions 



Standard atmospheric oscillations: no-oscillation hypothesis ruled out at 5a. 
VLI: |A„| upper limits of the order of 10~^^. 

MM search: upper flux limit of 1.4 • 10"^^ cm~^ s^^ sr^^ for 4 • 10"^ < /3 < 1. 
Nuclearite search: upper flux limit of 10"^^ cm~^ s^^ sr^^ for /? ~ 10~^. 
Charged Q-balls search: upper flux limit of ^ 10^^^ cm^^ s^^ sr^"'^. 
WIMP search: upper flux limit of ~ 10^^"' cm^^ s^^ sr^^. 
LIP search: upper flux limit of 6.1 • 10^^^ cm^^ s^^ sr^^. 
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Abstract 

The obtained exposure of the Baksan Underground Scintillation Telescope (BUST) is now 
one of the best for astrophysical challenges. We present a performance of the sample of 1102 
upward through-going muons collected by the BUST during 21.15 years of live time. We 
looked for an excess in arrival directions relatively to the expectations from atmospheric 
neutrinos and for annihilation signal of Dark Matter weakly interacting massive particles 
(WIMPs) inside the Sun. 



1. Introduction 

Nowadays searches for exotic physics serve mainly to resolve a " Double Dark" [1] 
picture of the Universe evolution. The neutrino telescopes may contribute by evi- 
dence of new neutrino properties or a new hypothetical heavy particle either neutral 
or charged, like WIMPs, monopole, or Q-balls. 

Recent data from WMAP [2] as well as the combined observations of 7, and 
X-rays in various experiments [3] point out with high precision on the existing of 
non-luminous, invisible and unknown form of matter filling in all space of the flat 
universe with only 5% of barions. 

We are focusing on searching for one of a favorite Dark Matter WIMP candidate - 
the neutralino (for a review see [4]). The probabilities for WIMP to be gravitation- 
ally trapped by the Sun and by the Earth were calculated theoretically [5]. It has 
been shown that WIMPs can scatter off a nucleus in the medium, sink to the central 
regions and accumulate there. During the evolution of the Universe the rate of neu- 
tralinos captured by the Sun could be large enough to provide a thermal equilibrium 
with their annihilation rate in the centre. Decays of generated fermions and bosons 
in neutralino annihilation channels produce high energy leptons from which only 
neutrinos survive on their way to detector. With benchmark points of mSUGRA 
we analyze their survival in comparison with our upper limits on annihilation rates. 



2. Experiment, Site and Triggers 

The Baksan Underground Scintillation Telescope is known as high energy neu- 
trino telescope of early generation in the big underground laboratory of Baksan 
Neutrino Observatory (BNO). The duration of continuous measurements covers 26 
years since 1978, while effective exposure amounts to 1.6 x 10^^ cm? s. 

The BNO site is in North Caucasus mountains at 1700 m above sea level and at 
the latitude of Mediterranean Sea: 42°41'E and 43°16'N. The BUST is located at 
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effective depth of ShQhg/cm? which is relevant to the 5x 10'^ reduction of background 
from the downward-going atmospheric muons. Trajectories of penetrating particles 
are reconstructed by the positions of hit tanks, which put together a system of 3,150 
liquid scintillation counters of standard type (70cm x IQcrn x 30cm) in configuration 
of parallelepiped (17m x 17m x 11m). The counters entirely cover all its sides and 
two horizontal planes inside at the distances 3.6 m and 7.2 m from the bottom. 
Each plane is separated from another one by concrete absorber of IQQg/cm?. The 
configuration provides 1.5° of angular accuracy. The time resolution measured with 
downward-going muons is equal to 5 ns. The arrival direction is determined by time- 
of-flight method, so that the measured value of (3 [6] is expected to be negative for 
upward-going particles. More details about the detector can be found in [6]. 

While expected ratio of downward-going muon flux to upward-going one is about 
10^, there were two hardware triggers used to select upward-going muons. They 
reduce initial rate of downward-going muons approximately by factor of 10^. Trigger 
1 covers the zenith angle range 95° ~ 180° while trigger 2 selects horizontal muons 
in the range 80° ~ 100°. The hardware trigger efficiency of 99% has been measured 
with the flux of atmospheric muons. These two triggers give about 1,800 events per 
day for further processing. In the year 2000 the telescope data acquisition system 
was upgraded and allowed to simplify the trigger system down to one general trigger 
with the rate of 17 Hz. All raw information is than undergone further selection using 
off-line reconstruction code. 

There are two simple requirements in data selection to be satisfied: only single 
trajectory of penetrating particle and a negative measured value of /3. However 
all events with negative values of (3 have been scanned by eye to check possible 
misinterpretation. 

It was found finally 1439 events survived these cuts and collected from December 
of 1978 until April of 2006 with 21.15 years of live time. 

3. Sample of upward through-going muons 

The reconstructed muon trajectories are required to have more than 7 m of 
measured length inside the telescope to reject upward-going muons that could be 
mimicked by downward-going atmospheric muon interactions or muon groups. So 
far we cut off particles with energy below 1 GeV [7]. In total 1102 events survived 
all cuts. From Fig.l one can conclude about continuous and stable rate of mea- 
sured upward through-going muons during all years of observations. The obtained 
zenith angular distribution of these events is shown in Fig. 2 in comparison with two 
hypotheses of simulated events: with and without two-flavours oscillating neutrino. 

The complete Monte Carlo simulations of both atmospheric neutrino interactions 
and the detector response on induced upward through-going muons have been done 
on base of author codes. In part of neutrino interactions we used the Bartol atmo- 
spheric neutrino flux [8] , the GRV-94 parton density distributions [9] and energy- loss 
parameterization [10] for muon propagation in surrounding rock. Induced muons in 
the generated neutrino interactions have to pass the same hardware trigger require- 
ments and the same set of cuts as real data. In the simulations of detector response 
the spread in threshold setting (w 10%), PMT's gain (w 10%), and time off-set of 
PMT's (w 2ns) as well as the part of dead tanks (sa 1%) have been accounted. By 
varying the detector parameters relevant to observed upward-going muon events 
the systematic uncertainties have been evaluated. 
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Cos(0) 

Fig. 1. Integral number of measured upward Fig. 2. Zenith distributions of upward through— 
through-going muons versus years. going muons: data (soHd) and MC simulations 

with oscillations (dash) and without ones (dot). 



Finally we have found the averaged number of 1268.5 events survived all cuts, 
corresponding to 22 runs of real data taking (in total of 460 years). The comparison 
of observed and simulated distributions in number of hit tanks as well as in thickness 
of absorber crossed by trajectories shows good agreement. Mean energy of generated 
neutrinos for muons going out through the telescope is about 50 GcV. 

All-in-all, the ratio of observed total number of events to expected one without 
oscillations is found to be 0.87 ± 0.03(stat.) ± 0.05(syst.) ± 0.15{theor.). 



4. Performance, Astrophysics Aspects and Limits 

The analysis of upward through-going muons at the BUST for searching for 
high energy neutrinos other than atmospheric origin has been done previously, 
when a live time of observations were approximately twice less. Corresponding flux 
limits have been obtained for astronomical point-like sources [6] and for WIMPs 
annihilations in the Sun and the Earth' core [7] . Now with selected upward through- 
going muons for 21.15 years of live time we have analyzed arrival muon directions 
in sky map by bootstrap method in searching for specified and unidentified sources 
and by optimal cone-size method for WIMP search from the Sun. 

For neutrino energy range of reconstructed events mentioned above the neutrino- 
induced muons deviate from parent neutrino direction inside cone of 4 degree with 
90% efficiency. In Fig. 3 the event number distribution of upward through-going 
muons is shown, where the sky map is in galactic coordinates and sizes of bin are 
4° X 4°. To get a sky map in standard deviations for each point, we obtained the 
expected noise by swapping the event time and event direction randomly and using 
Poisson statistics with sigma of 4 degree per bin. The obtained sky map in sigma 
deviations is shown in Fig. 4. There is no bin with sigma exceeding the value of 
3.6 and it is the evidence of the absence of statistically significant excess in any 
direction at 90% c.l. We have set new upper limit on muon flux from the Galactic 
Centre to be equal to 6.3 x 10~"'^^cto~^s~^. The limits for some known sources 
among pulsars, SNRs, AGNs, etc will be published later. 
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Fig. 3. Arrival directions in galactic coordinates Fig. 4. Standard deviations in galactic coordi- 



with binning 4° X 4° 



nates with binning 4° X 4° 



According to latest theoretical estimates the Sun is most promising somce in the 
search for hidden processes of Dark Matter than others (the Earth or G.C.). We 
have analyzed the sample of 1102 muon arrival directions in correlation with the 
Sun position. The distribution is shown in Fig. 5 together with averaged histogram 
of angles between of muon arrival directions and 5 fake suns. In absence of any 
features and excess toward the Sun we have set upper limits on number of muons 
at 90%c.l. for each cone sizes < 16° and improve our previous numbers by factor 
of 1.5 approximately. 

The optimal angular windows to collect 90% of muons from annihilation source 
for different neutralino masses are shown in Table 1 in [7] . They have been found by 
using the code wSUSY4 [7] to model expected neutralinos, and by further Monte 
Carlo simulations down to the detector response. The effects from oscillating neu- 
trino propagations in absorption medium were not included yet. However we intend 
to clarify their consequences on sizes of the optimal angular windows. 

Due to low energy threshold of the detector we can look for signatures from 
WIMPs with masses starting from magnitude about 12 GeV e.g. much smaller 
than one from LEP limits, if we suppose WIMP could be non supersymmetric 
particle. 

We have compared the obtained upper limits at 90% c.l. on annihilation rates 
in the Sun with DarkSusy [11] calculations for twelve known benchmark points 
of mSUGRA, and also the lightest neutralino mass model ( 53 GeV) and declared 
phenomena signature in EGRET analysis ( 64 GeV) [12]. The results are plotted 
in Fig. 6 showing that all these models are still not excluded. 

In summary, the Baksan data do not show any excess in measured sky-binning 
muon fluxes, including correlation with the Sun. The reached exposure of 5 x 
lO'^'^km^ X yr allows to improve our upper limits on muon fluxes and to get them 
down to limits from Super-Kamiokande [13]. Clearly, the exposure of several square 
kilometers in year as well as decreasing of current uncertainties in nuclear physics 
and astrophysics are required to produce adequate results compared to ones in 
standard cosmology. 
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Fig. 5. Distribution of angle (eps) between 
upward through-going muons and the Sun. 
Shaded histogram is averaged one for five fake 
positions of the Sun. 
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Abstract 

This note gives a short review of the statistical issues concerning upper hmit calculation 
and claiming of discovery arising in the search for exotic physics with neutrino telescopes. 
Low sample sizes and significant instrumental uncertainties require special consideration. 
Methods for treating instrumental or theoretical uncertainties in the calculation of limits 
or discovery are described. Software implementing these methods is presented. The issue 
of optimization of analysis cuts and definition of sensitivity is briefly discussed. 



1. Introduction 

Searches for exotic physics in neutrino telescopes pose challenges to statistical 
inference. In neutrino telescopes often physicists have to deal with low event counts 
and significant instrumental uncertainties which need to be taken into account. 
The small sample size makes the use of asymptotic methods doubtful. Presence of 
significant instrumental or theoretical uncertainties introduce nuisance parameters 
which need to be considered in the estimation of physics parameter intervals or in 
claims for discovery. 

Neutrino astronomy has so far been in an infant state, the main aim of existing 
experiments being proof of principle with little hope of detecting neutrinos of cosmic 
and/or exotic origin. Analyses have therefore often been optimized for setting the 
most stringent upper limit. With the advent of the IceCubc detector this situation 
might change, in the sense that discovery should become more likely. Analysis 
optimization methods should therefore be reviewed. 

In the next section, we remind the reader of methods to calculate confidence 
intervals and to claim discovery. In section 3 we present two methods for including 
uncertainties, followed by a discussion on analysis optimization. 

For simplicity, unless otherwise stated, throughout the paper we will consider a 
statistical process with probability density function (PDF) P{n\s), i.e. the proba- 
bility to observe n given the parameter s. 

2. Claiming discovery and calculating confidence intervals 

Though often mixed up [1,2], claiming discovery and calculation of confidence 
intervals are generally presented as two different cases in statistical literature, see 
e.g. [3] for an exhaustive introduction. 

Claiming discovery Claiming discovery is an example of hypothesis testing. In the 
search of exotic physics the null hypothesis, Hq, is usually "No signal present" and 
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the alternative hypothesis, Hi, is "Signal of exotic physics present". Discovery is 
claimed if the measurement result would be very unlikely under the condition that 
the null hypothesis is true. This is quantified by calculating the p-value: 

p^PiT>Tobs\Ho) (1) 

where T denotes a test statistics, which is a function of the observation and the 
involved hypotheses (see below) , and P is a probability density function describing 
the distribution of the test statistics under the null hypothesis. If the p-value is 
smaller than some predetermined threshold, asign, then the null hypothesis is re- 
jected. A particular common choice of the threshold is the probability corresponding 
to 5ct in an equivalent Gaussian distribution, i.e. asign ^ 3 • 10"^. 

The test statistics which in a statistical sense is most optimal in distinguishing two 
alternative hypotheses is according to the Neyman-Pearson lemma the likelihood 
ratio: 

^ C{n\Ho) 

C{n\Hi) ^ > 

One useful property of the likelihood ratio is that under the assumption that the 
null hypothesis is true, asymptotically the following statement holds: 

-21nr~x' (3) 

i.e. minus twice the logarithm of the likelihood ratio follows a distribution. 
Though this property is often used in particle physics it is not always true. Except 
for normality, there are a series of other requirements, some of which are violated 
even for quite common cases [2,4]. An experimenter intending to use the 'x^ dis- 
tribution should therefore be convinced that this is justified and in case of doubt, 
calculate the distribution of the test statistics under the null hypothesis for example 
by using Monte Carlo simulations. 

The threshold probability in order to decide if the null hypothesis should be 
rejected is based on the number of false detection one is willing to accept and on 
the number of searches which are performed. For example, in the context of one 
of the large CERN experiments, there might be about 50000 independent channels 
[5], i.e. 500 channels, 1000 resolution elements each. This implies a false positive 
detection rate of about 1.5 % if one requires a 5(t detection. It is not obvious, that 
a similar reasoning would lead to the same requirement in case of searches with 
neutrino telescopes. 



Calculation of confidence intervals Consider s a fixed (unknown) parameter. In 
frequentist statistics^ , one defines a confidence interval [si,S2] as a member of a 
set of intervals for which: 

P{s e [si, S2]) = 1 - a , for all s (4) 

These sets are to be constructed using only P{n\s), s being a fixed but unknown 
parameter and not P{s\n). The latter would make the method Bayesian, since s is 
treated as a random variable. The frequentist construction has been introduced by 
Neyman [6]. A special case, which is very common in high energy physics has been 



^ the Bayesian calculation is somewhat simpler, since here the interval can be found by integrating 
in the probability density function P(s\n), where the parameter is considered a random variable 
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proposed by Fcldman & Cousins [7]. Other common approximative methods exist 
which make use only of the hkchhood function, e.g. [3,8]. 

An frequentist methods which fulfiU the condition in equation 4 arc said to have 
coverage ^ . Methods which do not have coverage are not vahd or need correction . 
In case of doubt, coverage has to be demonstrated for example by using Monte 
Carlo simulations. 

3. Treatment of nuisance parameters 

Nuisance parameters are parameters which enter the data model, but which 
are not of prime interest. The probably most common example is the expected 
background in a Poisson process. In a usual physics experiment, only confidence 
intervals on the parameter of primary interest (for example signal flux or cross- 
section) are of interest. Thus, ways have to be found to marginalize the nuisance 
parameter. There are two common approaches: 

In the first method, the PDF without uncertainty in nuisance parameters is 
replaced by one where there is an integration over all possible true values of the 
nuisance parameter (integration method): 



Here btme is the true value of the nuisance parameter and best is its estimate. Since 
the integrated PDF is describing the probability of the true value given its estimate 
(and not vice versa) this method is Bayesian. Some prior probability distribution 
of the true value of the nuisance parameter has to be assumed. 

In the other common method, the PDF is replaced by one where for each s the 
PDF is maximized with respect to the nuisance parameters (profiling method) 



with notation as above. This method is completely frequentist, since it never treats 
btrue as a random variable. Therefore the argument of the maximisation is a likeli- 
hood function and not a PDF. 

Both approaches have recently been subjected to detailed studies regarding their 
coverage, i.e. with applications to confidence intervals [9-11]. For typical problems 
arising in neutrino telescopes they perform satisfactory. In the context of the LHC 
searches there are indications the Bayesian method under-covers badly, whereas the 
profiling method still seems to work fine [1] . 

Software has been developed which implements the integration and profiling 
method for typical problems arising in neutrino telescope analyses. 
pole++ is a C-I--I- library of classes which was developed based on the method 
presented in [12]. It allows calculating confidence intervals using the Feldman & 
Cousins method with integration of nuisance parameters and coverage studies. The 
signal process considered is a Poisson with known (possibly uncertain) background 
with different models for nuisance parameters. Several experiments with correlated 



^ To be more precise, in frequentist statistics in addition the probability has to be defined in 
terms of repeated identical experiments. 

^ In general, physicists tend to accept over-coverage (meaning, being conservative), whereas under- 
coverage is considered unacceptable. 




'true 



(5) 



P{n\s, btrue) > ma.xC{n\s, btrue) 



(6) 
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or uncorrelated uncertainties in the nuisance parameters can be combined. The 
pole++ library can be obtained from http://cern.ch/tegen/statistics.html 

TRolke is class which is part of the ROOT analysis package [13]. It treats a Poisson 
signal with background process with seven different models for nuisance parameters, 
which are marginalized using the profiling method. For extraction of the confidence 
intervals the likelihood function is used. For a description of the method see [10]. 

It should be noted that the MINUIT package [14] (in particular using the MINOS 
facility) also applies the profiling method and is completely general as long as the 
likelihood function can be written down. For the special case of a Poisson process 
TRolke applies some improvements. 

Though profiling and integration are here presented in the context of confidence 
intervals, both approaches can also be used in hypothesis testing. The substitutions 
in equations 5 and 6 then have to be applied to the likelihoods. The statistical 
properties to be studied in this case would be the distribution of the test statistics 
under the null hypothesis as well as possibly the power function (see next section). 

4. Analysis optimization 

Analyses are optimized defining some figure of merit (FOM), which will be max- 
imized (or minimized) with respect to some cut value t 

In searches with neutrino telescopes, it is often chosen to try to set the most 
stringent upper limit leading to the introduction of Model Rejection Factor [15]: 

MRF = ^'-^ ^ (7) 
ris 

Here, Si_q denotes the (1-a) confidence level upper limit on s. The mean is taken 
over the Poisson distribution with no signal. In case of presence of uncertainties in 
nuisance parameters the corresponding mean becomes: 

/>oo />oo 

< Si_„ / / dndbestSl-aP{n\best)P{best\btrue) (8) 

Jo Jo 

Unless the size of the uncertainties depends on the cut value or assymetric PDFs 
for the nuisance parameters have to be considered, the mean calculated using only 
the Poisson distribution should be sufficient for the calculation of the FOM. 

In case a decision on whether to set a limit or claim discovery beforehand is not 
desired, the figure of merit will have to be based on a sensitivity region which is 
meaningful both if the experimenter would report a limit and if the experimenter 
wants to claim discovery. Punzi [16] suggests to define the sensitivity region by: 

> 1-acL (9) 

Where asign is the significance which the experimenter requires to claim discovery 
and acL is the confidence level required in case the experimenter wants to calculate 
a limit. P(s) is the so called power function. Power is a concept arising in hypothesis 
testing. The power of a test is the probability to reject the null-hypothesis given 
the alternative hypothesis is true. The concept can also be applied to confidence 
intervals [11]. This definition of sensitivity fulfills several desirable properties: for 
example if s is inside the region defined by equation 9, then there is a probability 
of at least 1 — acL that it will be discovered. From the sensitivity a FOM can be 
calculated. Simple expressions of the FOM for common problems can be found in 
[16]. 
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5. Conclusions 

Major statistical challenges to be faced in searches for new physics with neutrino 
telescopes are low statistics, large instrumental uncertainties and the desire to de- 
tect an unknown process while at the same time being able to set stringent upper 
limits. 

In this note we point out and briefly discuss different solutions to the above chal- 
lenges. There are several methods existing to calculate confidence intervals in the 
presence of instrumental uncertainties which have been tested and behave well also 
for the small sample sizes and relatively large uncertainties usually encountered in 
neutrino astronomy. Code performing the necessary calculations is readily available. 
Sensitivity regions linking the power of a hypothesis test with the confidence level 
of a confidence interval yield figure of merits which can be used to optimize analysis 
with respect to both discovery and stringent limits. 
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Abstract 

In this talk, I will discuss the potential for the direct and indirect detection of Kaluza-Klein 
dark matter found in models with universal extra dimensions. Although the prospects for 
direct detection are somewhat daunting, those for future kilometer-scale neutrino tele- 
scopes and satellite-based cosmic ray detectors are very encouraging. Experiments such 
as IceCube, PAMELA and AMS-02, will be quite sensitive to dark matter in this model. 



1. Introduction 

Although there exists an enormous body of evidence for the existence of dark 
matter, its identity remains unknown [1]. Weakly Interacting Massive Particles 
(WIMPs) are, perhaps, the most well motivated class of candidates for dark matter. 
Among these, the lightest ncutralino in models of supcrsymmctry is the most widely 
studied. 

Models with extra spatial dimensions can provide an alternative candidate for 
dark matter, however. In particular, in models in which all of the Standard Mod- 
els fields arc free to propagate in the bulk, called universal extra dimensions, the 
Lightest Kalzua-Klcin Particle (LKP) may be stable and a potentially viable dark 
matter candidate [2,3]. 

The most natural choice for the LKP is the first Kaluza-Klein excitation of the 
hypercharge gauge boson, B^^\ sometimes called the "KK photon". I will simply 
refer to this state as "Kaluza-Klein Dark Matter" (or KKDM) throughout this talk. 
Previous studies of KKDM have found that the relic density predicted for such a 
state would naturally coincide with the measurements of WMAP for masses near 
about mLKP ~ 900 GeV if no other Kaluza-Klein states participate in the freeze- 
out process [2]. If other states are light enough to significantly effect this process, 
however, the LKP can be substantially lighter or heavier why still generating the 
observed dark matter abundance [2,4,5]. 

In this talk, I will review some of aspects of the direct and indirect detection 
of KKDM. The prospects for indirect detection are particularly promising in this 
model. 

2. Direct Detection 

KKDM particles can elastically scatter via scalar couplings with nuclei through 
the exchange of Higgs bosons and KK quarks. These contributions lead to a spin- 
independent KKDM-nucleon elastic scattering cross section approximately given 
by [6]: 
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where 

the KK quarks and the LKP. The various quantities, /, account for the quark and 
gluon content of the nucleon. 

Currently, the strongest constraints on the WIMP-nuclcon clastic scattering are 
on the order of 10~^ pb [7]. The range of cross sections predicted for KKDM, in con- 
trast, are far below the experimental sensitivity. Although future ton-scale detectors 
may be able to test such models, in the near future KKDM will go unexplored by 
direct detection experiments. 



3. Indirect Detection 

For the purposes of indirect detection, KKDM has several attractive phenomeno- 
logical features. Firstly, approximately 60% of KKDM annihilations are to charged 
lepton pairs (20% to each generation). 33% of annihilations produce pairs of up- 
type quarks and 3.6% produce neutrino pairs. The remaining fraction generate 
down type quarks and Higgs bosons. This is in stark contrast to neutralinos which 
do not annihilate efficiently to neutrinos, positrons, muons or other light fermions. 

Secondly, the total annihilation cross section for KKDM is given by 

95gf _ 1.7 X 10-26 cm3/s 

Notice that this consists entirely of an a-term in the expansion, < av >= a + bv^ + 
0{v'^), thus the low velocity cross section is naturally the maximum possible for a 
thermal relic. 

Furthermore, as we will see in the following section, KKDM's spin-dependent 
elastic scattering cross section with protons can be quite large, leading to the ef- 
ficient capture of such particles in the Sun, and correspondingly large neutrino 
fluxes. 



3.1. Indirect Detection with Neutrino Telescopes 



Dark matter particles travelling through the Galactic halo can occasionally scat- 
ter and become trapped in deep gravitational wells, such as the Sun or Earth. 
Within these bodies, they accumulate and their annihilation rate is enhanced, po- 
tentially providing an observable flux of high-energy neutrinos [8] . The capture rate 
of KKDM particles in the Sun is given by [9] 

where crH,SD is the spin-dependent, elastic scattering cross section of KKDM off of 
hydrogen. This expression assumes a local dark matter density of 0.3 GeV/cm'^ and 
a RMS velocity of 270 km/s. The elastic scattering cross section is given by [10] 
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Fig. 1. The rate of neutrino-induced muons above 50 GeV predicted in a kilometer-scale neutrino 
telescope, such as IceCube. Curves are shown for Kaluza-Klein quarks 10%, 20% and 30% heavier 
than the LKP. 
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(mg(i) — mLKp)/™LKP is the fractional shift of the Kaluza-Klein quark 
masses over the LKP mass and the A^'s parameterize the fraction of spin carried 
by a constituent quark q. Inserting numerical values for the A^'s, one finds 

1 o in-6 , /lOOOGeVy / 0.1 V 
cth.sd « 1-8 X 10 "pb I — ) . (5) 



\ n^LKP / \rq(i) / 

For the annihilation and elastic scattering cross sections of KKDM, the annihilation 
rate in the Sun should reach (or nearly reach) equilibrium with the capture rate. 
These annihilations can produce neutrinos directly or in the decays of tau Icptons 
or quarks [9] . 

Muon neutrinos which reach the Earth from the Sun can scatter in charged 
current interactions with nucleons to produce high-energy muons. These muons 
produce observable "tracks" as they propagate through the medium of a neutrino 
telescope, such as the Antarctic ice of the IccCube experiment. The rate of muon 
tracks generated in a kilometer-scale neutrino telescope from KKDM annihilations 
in the Sun in shown in figure 3.1. Notice that these results depend strongly on the 
LKP's mass and the mass of the Kaluza-Klein quarks. Calculations of the radiative 
corrections to the Kaluza-Klein spectrum estimate values of r^(i) roughly in the 
range of 0.1 to 0.2 [11]. For an 800 GeV LKP. about 5 to 50 events per year are 
predicted in IccCube over this range. For a lighter LKP of 500-600 GeV, up to 100 
events could be observed. 



3.2. Indirect Detection with Positron Experiments 

Dark matter annihilating in the Galactic halo can produce several potentially ob- 
servable species of cosmic rays, including gamma-rays, anti-protons, anti-deutcrons 
and positrons. Kaluza-Klein Dark Matter has characteristics which are favorable 
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Fig. 2. The ratio of positrons to positrons plus 
electrons with a contribution from annihilating 
KKDM, compared to the HEAT data. Results 
for several choices of diffusion parameters are 
shown. The annihilation rate was normalized 
to the HEAT data. 



Fig. 3. The sensitivity of PAMELA and 
AMS-02 to cosmic positrons generated in 
KKDM annihilations. Contours are shown for a 
completely homogeneous distribution (BF=1) 
and a mildly clumped distribution (BF=5). 
The downward sloping curve represents the an- 
nihilation cross section predicted for KKDM. 



for detection with cosmic positrons: a large low-velocity annihilation cross section 
and a large fraction of annihilations which produce energetic positrons (such as the 
modes e+e" 



To calculate the observed positron spectrum, the positrons injected in dark mat- 
ter annihilations must be propagated through the Galactic magnetic fields, includ- 
ing scattering with starlight and the CMB. These effects are taken into account by 
solving the diffusion-loss equation: 
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+ Q{E,^,x), (6) 



where K{Ep+,x) is the diffusion constant, b{E^+,x) is the energy loss rate and 
Q{Ep+,x) is the source term. 

The HEAT experiment, during balloon flights in 1994-95 and 2000, observed an 
excess in the cosmic positron flux when compared to the electron flux. This excess 
peaks near 8-10 GcV and extends to above 30 GeV where the detector's sensitivity 
falls off [12]. It has been suggested that this excess could be the product of dark 
matter annihilations [13]. Neutralino dark matter, which does not annihilate sig- 
nificantly to light fermions, produces positrons inefficiently, however, and therefore 
requires a very high annihilation rate (i.e. a clumpy distribution) to account for 
the excess observed by HEAT [14]. KKDM, on the other hand, can produce such 
positron fluxes more naturally [15]. 

Satellite based cosmic ray experiments, such as PAMELA and AMS-02, will be 
capable of measuring the cosmic positron spectrum with much greater precision 
and to much higher energies than HEAT [16,17]. The first results from PAMELA, 
which began its mission in June of 2006, arc eagerly awaited. Experiments such as 
these will be very sensitive to the presence of KKDM in our Galaxy [16]. In figure 3, 
the sensitivity of these experiments to KKDM are shown. 
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4. Summary 

Kaluza-Klcin particles within the context of models with universal extra dimen- 
sions constitute a viable candidate for dark matter with a number of attractive 
and interesting phenomenological features. Although the prospects for the direct 
detection of such a state, at least within the next few years, arc not particularly 
encouraging, indirect detection techniques can be particularly powerful within such 
models. 

The prospects for the indirect detection of Kaluza-Klein Dark Matter (KKDM) 
are very promising for both future neutrino telescopes and cosmic positron experi- 
ments. Experiments such as IceCube, PAMELA and AMS-02 will provide powerful 
probes of KKDM in the coming years. With PAMELA currently taking data, the 
near future may well indeed reveal to us new insights into the nature of dark matter. 
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Abstract 

We present the results of a search for nearly vertically upward going neutrino-induced 
muons with 1038 live days data from the neutrino telescope NT200 taken over the five 
year period 1998-2002. No excess from WIMP annihilation in the Earth above the expected 
atmospheric neutrino background has been found. Upper limits at 90% confidence level 
have been derived on the moun flux induced by neutrinos from WIMP annihilation. 



1. Introduction 

The Baikal Neutrino Telescope is operated in Lake Baikal, Siberia, at a depth of 
1.1 km. Lake Baikal deep water is characterized by an absorption length of Lafcs(480 
nm)=20-^24 m, a scattering length of Lg =30-^70 m and a strongly anisotropic scat- 
tering function f{6) with a mean cosine of the scattering angle cos(6') = 0.85 ~ 0.9. 
The first stage telescope configuration NT200 [1] was put into permanent opera- 
tion on April 6th, 1998 and consists of 192 optical modules (OMs). An umbrella-like 
frame carries 8 strings, each with 24 pairwisc arranged OMs (see central part of Fig. 
1 left panel). All OMs face downward, with the exception of the second and eleventh 
pairs on each string which face upward. Three underwater electrical cables connect 
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Fig. 1. Left; The upgraded Baikal Telescope NT200+ : The old NT200 surrounded by three 
external long strings at 100 m radius from the center. Also indicated: external laser and DAQ 
center. Right: Detection area after WIMP cuts as a function of zenith angle. 

the detector with the shore station. Each OM contains a 37-cm diameter QUASAR 
- photo multipher (PM), which has been developed speciaUy for our project [2]. The 
PMs record the Cherenkov Ught produced by charged particles in water. The two 
PMs of a pair are switched in coincidence in order to suppress background from 
bioluminescence and PM noise. A pair of OMs defines a channel. The light arrival 
time assigned to a channel is the response time of the OM with the earliest hit. 
The amplitude assigned to a channel is that recorded by one pre-selected PM of 
the two PMs in a pair. 

A trigger is formed by the requirement of > hits (with hit referring to a chan- 
nel) within 500 ns. A^ is typically set to 3 or 4. For these events, amplitude and 
time of all fired channels are digitized and sent to shore. 

Two nitrogen lasers are used for the calibration of the detector. The first one {fiber 
laser) is mounted just above the array. Its light is guided via optical fibers of equal 
length to each OM pair. The fiber laser provides the OMs with simultaneous light 
signals in order to determine the offset for each channel. The second laser {water 
laser) is arranged 20 m below the array. Its light propagates through the water. This 
laser serves to monitor the water quality, in addition to dedicated environmental 
devices located along a separate string. A full cycle of detector calibration running 
both lasers over a wide range of intensities is repeated every third day. 

The upgraded telescope NT200+ was put into operation on April 9th, 2005 [3]. 
This configuration consists of the old NT200 telescope, surrounded by three new 
external strings (sec Fig.l). The external strings are 140 m long and are placed at 
100 meter distance from the center of NT200. Each string contains 12 OMs grouped 
in pairs like in NT200. The upper pairs are at approximately the same level as the 
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Muon energy, lg(I^ /GeV) cos(9) 



Fig. 2. Left: Energy distribution of muons produced by atmospheric neutrinos, after cuts. 
Dashed/full histogram include/neglect oscillations. Right: Angular distributions of selected neu- 
trino candidates as well as expected distributions in case of and without oscillations (solid and 
dashed curves respectively). 

bottom OMs of NT200. For NT200+, calibration is done with a powerful external 
laser light source with up to 5 x 10^'^ photons per pulse and nsec-pulse duration, 
which is located between two outer strings, see Fig.l. 

Our earlier limits on muon flux induced by WIMPs annihilation at the center 
of the Earth, which have been obtained with NT200 using data collected during 
1998-1999, were presented elsewhere [4]. Here we present the new results of a search 
for nearly vertically upward going neutrino-induced muons with 1038 live days data 
from detector NT200 taken over the five year period 1998-2002. 

2. Search for Neutrinos from WIMP Annihilation 

The search for WIMPs with the Baikal neutrino telescope is based on a possible 
signal of nearly vertically upward going muons, exceeding the flux of atmospheric 
neutrinos. The method of event selection relies on the application of a series of 
cuts which are tailored to the response of the telescope to nearly vertically upward 
moving muons. The cuts remove muon events far away from the opposite zenith 
as well as background events which are mostly due to pair and bremsstrahlung 
showers below the array and to bare downward moving atmospheric muons with 
zenith angles close to the horizon {6 > 60°). The candidates identified by the cuts 
are afterwards fitted in order to determine their zenith angles. 

For the present analysis we included all events with >5 hit channels, out of which 
>4 hits are along at least one of all hit strings. To this sample, a series of 5 cuts 
is applied. Firstly, the time differences of hit channels along each individual string 
have to be compatible with a particle close to the opposite zenith (1). The event 
length should be large enough (2), and the center of gravity of hit channels should 
not be close to the detector bottom (3). The latter two cuts reject efficiently brems 
showers from downward muons. Finally, also time differences of hits along different 
strings have to correspond to a nearly vertical muon (cuts 4 and 5). 
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Fig. 3. Left: Limits on the excess muon flux from tlie center of the Earth versus half-cone of the 
search angle. Right: Limits on the excess muon flux from the center of the Earth as a function of 
WIMP mass. 

Fig. 1 (right panel) shows the dependence of the detection area on the cosine of 
the zenith angle 6. The applied cuts select niuons with -1< cos{9) <-0.65 and result 
in a detection area of about 1800 for vertically upward going muons. 

The expected (normalized) energy spectrum of muons produced by atmospheric 
neutrinos (Bartol flux [5]), which survive all cuts, with oscillations (using Super- 
Kamiokande parameter set =2.5T0~"^ eV^ with full mixing. 9m ~ 7i'/4 [6]) and 
without oscillations is shown in Fig. 2 (left panel). The energy threshold for this 
analysis is E^^ir ~ 10 GeV. We expect a muon event suppression of (25-30)% due 
to neutrino oscillations. 

From 1038 days of effective data taking between April 1998 and February 2003, 
48 events with -1< cos(6') <-0.75 have been selected as neutrino candidates, com- 
paring to 56.6 events expected from atmospheric neutrinos in case of oscillations 
and 73.1 without oscillations. The angular distribution of these events as well as 
the MC - predicted distributions without (dashed curve) and with (solid curve) os- 
cillations are shown in Fig. 2 (right panel). Within la statistical uncertainties the 
experimental angular distribution is consistent with the prediction including neu- 
trino oscillations. Note, also, that the theoretical uncertainty in the atmospheric 
neutrino flux calculation is estimated to be about 20% [7] . 

2.1. Limit on muon excess from WIMP annihilation 

Regarding the 48 detected events as being induced by atmospheric neutrinos, 
one can derive an upper limit on the additional flux of muons from the center of 
the Earth due to annihilation of neutralinos - the favored candidate for cold dark 
matter. The 90% C.L. muon flux limits for six cones around the opposite zenith 
obtained with NT200 {E^Yit >10 ^eV) in 1998-2002 are shown in Fig. 3 (left panel). 
It was shown [8-10] that the size of a cone which contains 90% of signal strongly 
depends on neutralino mass. The 90% C.L. flux limits are calculated as a ftmction 
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of ncutralino mass using cones which collect 90% of the expected signal [8] and are 
corrected for the 90% collection efficiency due to cone size. These limits are shown 
in Fig. 3 (right panel). Also shown in Fig. 3 are limits obtained by Baksan [8], 
MACRO [9], Super-Kamiokande [10] and AMANDA [11]. 

3. Conclusion 

The Baikal neutrino telescope NT200 is taking data since April 1998. A search 
has been performed for nearly vertically upward going muons based on NT200 
telescope data from the five years 1998-2002. The number of detected events, as 
well as their angular distribution, is compatible with expectation from atmospheric 
neutrinos. No statistically significant excess of nearly vertically upward going muons 
over atmospheric neutrino background was seen. 

Limits on an excess of the muon flux due to WIMP annihilation in the center of 
the Earth for various cone angles around nadir as well as muon flux limits for dif- 
ferent neutralino masses have been derived, and compared with previous estimates 
by other experiments. 
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Abstract 

Data taken with the AMANDA high energy neutrino telescope can be used in a search 
for an indirect dark matter signal. If non-baryonic dark matter exists in the form of neu- 
tralinos a neutrino flux at energies lower than about 1 TeV is expected from the decay 
of neutralino pair annihilation products. We present published results from searches for 
neutralinos accumulated in the Earth and the Sun, using AMANDA data from 1997-1999 
and 2001 respectively. Current efforts to improve neutralino trigger and selection efficien- 
cies in analyses of multi-year AMANDA data samples are discussed. Finally, expected 
sensitivities for a best-case 1 km'^ IceCube scenario are given. 



1. Introduction 

Cosmological observations have long suggested the presence of non-baryonie dark 
matter on all distance scales. The WMAP results [1] confirmed our current under- 
standing of the Universe, summarized in the concordance model. In this model the 
Universe contains about 23% non-baryonic cold dark matter, but nothing is pre- 
dicted about the nature of this dark matter. A massive, weakly interacting and 
stable particle appears in Minimally Supersymmetric extensions to the Standard 
Model that assume R-parity conservation. Indeed, the supersymmetric partners of 
the neutral electroweak and Higgs Standard Model bosons mix into an interesting 
dark matter candidate, the neutralino, whose mass is expected in the GeV-TeV 
range [2]. On their trajectory through the Universe these particles will scatter 
weakly on normal matter and lose energy. Eventually, the dark matter particles 
will be trapped in the gravitational field of heavy celestial objects, like the Earth 
and the Sun [3]. The particles accumulated in the center of these bodies will anni- 
hilate pairwise. The neutrinos produced in the decays of the Standard Model anni- 
hilation products can then be detected with a high energy neutrino detector as 
an excess over the expected atmospheric neutrino flux. In this paper we present 
the published results of searches with the Antarctic Muon And Neutrino Detector 
Array (AMANDA) for neutralino dark matter accumulated in the Earth (1997-1999 
data set) and the Sun (2001 data set). We also summarize current improvements and 
show preliminary results from ongoing analyses on higher statistics data samples 
accumulated during recent years (2001-2003). 

The AMANDA detector [4] at the South Pole uses the polar ice cap as a Cherenkov 
medium for the detection of relativistic charged leptons produced in high energy 
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Fig. 1. (a) Detection efficiencies relative to trigger level for the different filter levels in the ter- 
restrial neutralino analysis {mx=250 GeV, hard spectrum) for 1997-1999 data, neutralino signal, 
atmospheric muons and neutrinos, see Ref. [9]. (b) As a function of neutralino mass, the 90% 
confidence level upper limit on the muon flux coming from hard neutralino annihilations in the 
center of the Earth compared to the limits of other indirect experiments [10] and the sensitivity 
estimated for a best-case IceCube scenario. Markers show predictions for cosmologically relevant 
MSSM models, the dots represent parameter space excluded by CDMS [11]. 

neutrino interactions with nuclei. The detector was constructed between 1996 and 
2000. Now totaling 677 light sensitive devices distributed on 19 strings, the detec- 
tor is triggered when at least 24 detector modules are hit within a sliding 2.5 ^s 
window. Before 2000, the detector configuration consisted of between 10 and 13 
strings and consequently a lower multiplicity trigger condition was able to cope 
with the high rate of events produced by cosmic ray interactions with the atmo- 
sphere. Reconstruction of muons, with their long range, offers the angular resolution 
required to reject the atmospheric background and search for a ncutralino-induccd 
signal, which, due to the geographic location of AlVIANDA, yields vertical upward- 
going (Earth) or horizontal (Sun) tracks in the instrumented volume. Indeed, it 
is possible to eliminate the dominant background, downward-going atmospheric 
muons. However, upward-going atmospheric neutrinos will always contaminate the 
final, selected data sample. IceCube [5], the 1 km'^ successor of AIMANDA, is cur- 
rently being deployed and will continue searching for dark matter in the coming 
10-15 years. 

2. Signal and background simulation 

Wc have used the DARKSUSY program [6] to generate dark matter induced 
events for seven neutralino masses between 50 GeV and 5000 GeV, and two 
annihilation channels for each mass: the W^W~ channel produces a hard neutrino 
energy spectrum (t+t~ for < JTivy), while hh yields a soft spectrum. The cosmic 
ray showers in the atmosphere, in which downward-going muons are created, are 
generated with CORSIKA [7] with a primary spectral index of 7=2.7 and energies 
between 600 GeV and lO^'^ GeV. The atmospheric neutrinos are produced with 
NUSIIVI [8] with energies between 10 GcV and 10^ GcV and zenith angles above 80°. 

3. Search for neutralino annihilations in the center of the Earth 

A neutralino-induccd signal from the center of the Earth was searched for in 
AlVIANDA data collected between 1997 and 1999, with a total effective livetime 
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of 536.3 days. To reduce the risk of experimenter bias, the complete data set of 
5. Ox 10^ events was divided in a 20% subsample. used for optimisation of the selec- 
tion procedure, and a remaining 80% sample, on which the selection was applied 
and final results calculated. Similarly, the sets of simulated events were divided 
in two samples: the first for use in the selection optimisation and the second for 
the selection efficiency calculations. The simulated atmospheric muon sample con- 
tains 4.2x10^ triggered events (equivalent to an effective livetime of 649.6 days). 
The sample of atmospheric neutrinos totals 1.2x10* events, which corresponds to 
2.2x10* triggers when scaled to the livetime of the analysis. 

First, we try to suppress the dominant atmospheric muon background which is 
about 10^ times more abundant than the atmospheric neutrino background. This is 
partially done by selecting the events that are reconstructed as upward-going and 
that satisfy a cut correlated with reconstruction quality (filter level 3). However, 
only a 10"'^ reduction of the atmospheric muons is obtained this way (Fig. la) 
and more elaborate selection criteria are needed to reject downward-going muon 
tracks misreconstructcd as upward-going. Depending on the detector configuration 
and the neutralino model under study, the characteristics of the signal differ, which 
influences selection efficiencies significantly at this point. Therefore, all further cuts 
are fine-tuned separately for each neutralino model and year of data taking. At filter 
level 4, a neural network selection is applied, reaching another 10^^ reduction. Filter 
level 5 cuts sequentially on observables, with the goal of removing downward-going 
muons that resemble signal events. 

At filter level 5 the data sample is dominated by atmospheric neutrinos (see 
Fig. la). With no significant excess of vertical tracks observed, the final selection 
on reconstructed zenith angle (filter level 6) was optimized for the average most 
stringent 90% confidence level upper limit on the muon flux. From the number of 
observed events and the amount of (simulated) background in the flnal angular 
search bin, we infer the 90% confldence level upper limit on the number of signal 
events. Combined with the effective volume at the flnal cut level and the livetime 
of the collected data, this yields an upper limit on the ncutrino-to-muon conversion 
rate, which can then be related to the muon flux [9] (see Fig. lb). 

4. Search for neutralino annihilations in the Sun 

The AMANDA data used in the search for solar neutralinos consists of 8.7x10* 
events, corresponding to 143.7 days of effective livetime, collected in 2001. In con- 
trast to the search in Section 3, reducing the risk of experimenter bias in this 
analysis can be achieved by randomizing the azimuthal angles of the data. The 
advantage of this procedure is that it allows the use of the full data set for cut 
optimization. The azimuthal angles are restored once the optimization is flnal- 
ized and results are calculated. The simulated atmospheric background sample at 
trigger level totals 1.6x10* muons (equivalent to 32.5 days of effective livetime) 
and 1.9x10** neutrinos. 

The solar neutralino analysis suffers the same backgrounds as the terrestrial neu- 
tralinos, but the signal is expected from a direction near the horizon, due to the 
trajectory of the Sun at the South Pole. This analysis was only possible after com- 
pletion of the AMANDA-II detector, whose 200 m diameter size provides enough 
lever arm for robust reconstruction of horizontal tracks. 

A similar analysis strategy as in Section 3 was adopted. First, events were selected 
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Fig. 2. (a) Detection efficiencies relative to trigger level for the different filter levels in the solar 
neutralino analysis {'mx=500 GeV, hard spectrum) for 2001 data, neutralino signal, atmospheric 
muons and neutrinos, (b) As a function of neutralino mass, the 90% confidence level upper limit 
on the muon flux coming from hard neutralino annihilations in the center of the Sun compared 
to limits of other indirect experiments [10] and the sensitivity estimated for a best-case IceCube 
scenario. Markers show predictions for cosmologically relevant MSSM models, the dots represent 
parameter space excluded by CDMS [11]. 

with well-reconstructed horizontal or upward-going tracks. The remaining events 
are passed through a neural network that was trained separately for the neutralino 
models under study and used data as background. Although a data reduction of 
~10~^ compared to trigger level is achieved, the data sample is still dominated 
by misreconstructed downward-going muons. Finally, these are removed by cuts on 
observables related to reconstruction quality. 

There is no sign of a significant excess of tracks from the direction of the Sun in the 
final data sample. The expected background in the final search bin around the Sun 
is estimated from off-source data in the same declination band, which eliminates the 
effects of uncertainties in background simulation. Combining this with the number 
of observed events, the effective volume and the detector livetime, we obtain 90% 
confidence level upper limits on the muon flux coming from annihilations in the 
Sun for each considered neutralino mass [12], as shown in Fig. 2b. 

5. Discussion and outlook 

Figures lb and 2b present the AlVIANDA upper limits on the muon flux from 
neutralino annihilations into a hard channel in the Earth and the Sun respec- 
tively, together with the results from other indirect searches [10]. Limits have been 
rescalcd to a common muon threshold of 1 GeV using the known energy spectrum of 
the neutralinos. Also shown are the cosmologically relevant IVISSIVI models allowed 
(crosses) and disfavoured (dots) by the direct search from CDIVIS [11]. Compared 
to our search for a terrestrial neutralino signal in 1997 AIMANDA data [13], the 
limit has been improved by at least a factor of 2, more than expected from addi- 
tional statistics alone. This is due mainly to the separate cut optimization for each 
neutralino mass, which exploits the characteristic muon energy spectrum of each 
model. 

In 2001 an additional trigger was installed that lowered the energy threshold of 
the detector. This trigger takes into account spatio-temporal correlations in the 
event hit pattern. A preliminary analysis with data taken between 2001 and 2003 
shows an improvement of a factor of about 30 in the effective volume in the search 
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for 50 GcV ncutralinos (soft annihilation channel) from the Earth with respect 
to the analysis presented in this note. We are currently performing searches for a 
dark matter signal both from the Earth and the Sun with data taken from 2000 
and later. The increased detector exposure combined with improved reconstruction 
techniques and the new trigger setting will result in improved sensitivities (a 4- year 
exposure alone would already give an improvement of a factor of two) . 

As of February 2006 the first nine IceCube instrumented strings are operational, 
monitoring a volume as big as that of AMANDA (0.022 km'^). First results show 
that the performance of the system meets or exceeds design requirements [14]. The 
final IceCube detector, to be completed around 2011, will encompass an approxi- 
mately hexagonal instrumented volume of 1 km'^. The sensitivity to dark matter 
signals will not only be improved by the factor 50 increase in fiducial volume, but 
also (especially for solar ncutralinos) by the greater horizontal extension compared 
to AMANDA. Furthermore, the smaller trigger efficiencies for low energy muons 
{Ef^ <100 GeV), due to sparse instrumentation, will be partially compensated for 
by incorporating the denser AMANDA array in IceCube. Assuming a best-case sce- 
nario (10 and 5 km'^yr exposure for Earth and Sun respectively, no misreconstructed 
atmospheric background and without incorporating systematic uncertainties) we es- 
timated the IceCube sensitivity for a terrestrial and solar neutralino-induced muon 
flux, shown in Fig. lb and 2b. The IceCube results will be complementary to future 
direct detection experiments, and even have an advantage over these for low mass 
solar ncutralinos with large spin-dependent scattering cross sections. 
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Abstract 

In supersymmetric theories where the hghtest supersymmetric particle is the gravitino the 
next to hghtest supersymmetric particle (NLSP) is typically a long lived charged slepton. 
These NLSPs can be produced by high energy neutrinos interactions with nucleons in the 
Earth and be detected by km'^ neutrino telescopes. The signal, consists of two parallel 
charged tracks separated by a few hundred meters. This is compared to the main back- 
ground, coming from direct di-muon production. The distance between the background 
tracks is much smaller and allows for a clean separation from the NLSP signal. We conclude 
that neutrino telescopes will complement collider searches in the determination of the su- 
persymmetry breaking scale, and may even provide the first evidence for supersymmetry 
at the weak scale. 



1. Introduction 

One of the most attractive candidate theory for the extension of the standard 
model of particle physics is weak scale supersymmetry (susy) . The particle spectrum 
of this model is determined by the susy breaking mechanism {'s/F). When susy is 
broken at high scales such that a/F > lO^^GeV the LSP is typically the neutralino. 
If however susy is broken at lower scales, ^/F < 10^°GeV, the LSP tends to be 
the gravitino and the Next to Lightest Supersymmetric Particle (NLSP) is usually 
a charged slepton, typically the right-handed stau. Within these models, if VF is 
much larger than a TeV the stau NLSP decay is suppressed [1] and its lifetime can 
be very large. 

It was recently proposed [1] that the diffuse flux of high energy neutrinos collid- 
ing with the Earth can produce pairs of slepton NLSPs. The energy loss of these 
particles is very small and they travel for long distances before stopping. This com- 
pensates their small production cross section since they can be detected far away 
from their production point. 

The NLSP production typically has a high boost and they will not decay inside 
the Earth provided the susy breaking scale is VF > lO^GeV. Since the NLSP is 
charged; its upward going tracks can be detected in large ice or water Cerenkov 
detectors, such as IceCube. 

Wc performed a Monte Carlo simulation of the NLSP production and propagation 
through the Earth and determined their signature in km^ neutrino telescope [2]. 
The main background (di-muon events) was also simulated and compared to the 
NLSP signal. 
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2. NLSP Production and Propagation Through the Earth 

The susy processes for NLSP production is analogous to the standard model 
(SM) charged current (CC) interaction and involves a t-channel production of a left- 
handed slepton {II) and a squark (q) through a gaugino exchange. In the dominant 
process the gaugino is a chargino and in the subdominant process a neutralino. The 
II and q prompt decay results in two lighter right-handed sleptons (NLSPs) plus 
SM particles. The NLSP will always be produced in pairs and is typically the stau 

The energy threshold for the NLSP production is given by the II and q masses. 
We take the typical values of 250 GeV for both chargino and II masses, 150 GeV 
for the NLSP mass and set the q mass to three values of 300, 600 and 900 GeV. 

We compute the NLSP cross section including both dominant and subdominant 
processes. The susy cross section is about three orders of magnitude lower than the 
SM production cross section [2]. 

For the high energy neutrino flux reaching the Earth, we take the Waxman- 
Bahcall [3] limit 

dcf,,\ (1-4) X 10-« ^ _2 _i . 

Since the NLSP production is independent of the neutrino flavor, the initial neutrino 
flux contains both z/^ and z^g in a 2 : 1 ratio. 

The NLSP propagation through the Earth depends on its energy loss which at 
high energies is dominated by radiative losses. The main processes are bremsstrahlung, 
pair production and photo- nuclear interactions [2,4]. As a result the NLSP will 
travel distances much greater than a muon which will compensate for the lower 
cross section. 



3. Simulation of NLSP Production, Propagation and Detection 

In order to analyze the NLSP production, propagation and detection we de- 
veloped a Monte Carlo simulation. Assuming an isotropic incoming neutrino flux 
we generated 30K NLSP events for each q mass. The neutrino incoming energy is 
distributed in steps ranging from the stau production threshold to 10^^ GeV. 

The survival probability (Pg) for a neutrino traveling through a path dl is given by 
Ps = exp(J' nadl) where n is the Earth number density [5] and a is the interaction 
cross section. The probability of interaction is 1 — Pg. 

To simulate the NLSP production an interaction point is randomly chosen from an 
interaction probability distribution. If the interaction point falls within a distance 
from the detector that is smaller than the NLSP range the event is accepted. If 
this distance is greater than the NLSP range the event will only count for the 
normalization. 

The center of mass (CM) angular distribution was chosen based on the differential 
cross section distribution. We assume that the CM angular distribution of the right- 
handed slepton pair (NLSPs) is the same as the CM angular distribution of the II 
and q which is a good approximation [2]. The 4-momentum of the NLSP pair is 
determined in the CM from the angular distribution and boosted to the lab frame. 
The lab angular {&iab) distribution is determined from the lab 4-momentum. 

With this procedure we can determine not only the rate of events at the detector 
as well as the two NLSP track separation in the detector. The separation is simply 
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Fig. 1. Energy distribution of £ii pair events per km^, per year, at the detector. Curves that do 
not reach y axis; from top to bottom: rriq = 300, 600 and 900 GeV. Here, = 150 GeV and 

= 250 GeV. Also shown are the neutrino flux at earth and the fj, and the di-muon flux through 
the detector (curves that reach y axis; from top to bottom respectively). In all cases we make use 
of the WB limit for the neutrino flux. 

Olab times the distance from the neutrino interaction point and the detector. 

4. Background 

Muons produced in the Earth could be a potential source of background for the 
NLSPs events. They can however be eliminated by requiring two charged tracks 
in the detector. The remaining important background is di-muon events origi- 
nated from charm decay. Charm is produced from high energy neutrino interactions 
through the following process: — > fi"Hc ii~ fj,^ H^i/, where the charm hadron 
decays according to H^fi^v and Hj. can be a strange or non-strange quark. 

We computed the charm production cross section a{vN — > cX) from a d or s 
quark [2] . Although this cross section is around an order of magnitude greater than 
the one for NLSP production [2], muons lose much more energy while traveling 
through the Earth and therefore have to be produced very close to the detector. 
For this reason their track separation in the detector will also be much smaller. 

We performed a Monte Carlo simulation of the di-muon production, propagation 
and detection using the same procedure as for the NLSP. 

5. Results 

The energy distribution of NLSPs and di-muon background are shown in Figure 1. 
The number of events per km^ per year is shown in Table 1 where not only a 
neutrino incoming rate equal to the WB limit is assumed but also one equal to the 
Mannheim. Protheroc and Rachcn (MPR) limit [6]. 

Figure 1 and Table 1 show that km'^ neutrino telescopes are sensitive to NLSP 
detection. Although the number of di-muon events is larger than the NLSP rate 
there are many ways to reduce this background. 

First of all the lower integration limit to determine the number of events shown 
in Table 1 is different for the NLSP and the background. It is the NLSP production 
threshold energy for the number of NLSPs and 10^ GeV for the di-muon back- 
ground. The reason for this is that as the NLSPs lose less energy than muons, they 
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Tabic 1 

Number of events per km^ per year assuming the WB and MPR limits. The first column refers to 
di-muon events. The last three columns correspond to NLSP pair events, for three different choices of 
squark masses: 300 GeV, 600 GeV and 900 GeV. The number of di-muon events are given for energies 
above 10"^ GcV and of NLSPs above threshold for production. 
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Fig. 2. Arrival energy distribution of tlic Iji at the detector and for = 300, 600 and 900 GeV.. 
Here, m| = 150 GeV and = 250 GeV. Also shown is the arrival distribution for the di-muon 
background. The energy deposited in the detector by a stau traveling the average track length of 
800 m [7] is E^^"^ = 150 GeV, approximately the same for all the masses considered here. 

will look like lower energy muons in the detector. The NLSP energy deposition in 
the detector is 150 GeV (see figure ( in [2]). For the same reason, one can cut higher 
energy events removing events that deposit more than 300 GeV and reduce most of 
the dimuon background. Figure 2 shows the arrival energy at the detector for both 
NLSP and di-muon background where the efficiency of such a cut is clear. 

Another powerful way for background reduction is the track separation between 2 
NLSPs in the detector. This is shown in Figure 3. While a good fraction of the NLSP 
events are more than 50 m separated the dimuons are less than 50 m separated. 

6. Conclusions 

We conclude that Km^^ neutrino telescopes have the potential to discover the 
NLSP if susy models that predict the gravitino as the LSP and a charged slcpton 
as the NLSP are correct. This will be an indirect determination of the dark matter. 
If the NLSP is observed it will constitute a direct probe of the susy breaking scale. 
This search is complementary to the LHC. 
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Fig. 3. Top panel: Track separation distribution of In pair events. Bottom panel: Tlie track 
separation distribution of the di-muon background. The relative normalization corresponds to the 
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different binning between the two figures. 
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Abstract 

A number of signals involving charged cosmic rays and high-energy photons have been 
interpreted as being due to annihilating dark matter. This article provides an overview 
of the experimental evidence and discusses in particular detections of antiprotons and 
positrons in the cosmic radiation, the diffuse 7-ray emission between 10 MeV and 100 GeV 
from the Milky Way, and the 511 keV annihilation radiation and the flux of very high- 
energy photons (> 100 GeV) from the Galactic Centre. 



1. Introduction 

In contrast to laboratory-based dark matter searches (see [1] for a review of direct 
methods), indirect searches do not aim at the dark matter particles as such, but 
try to measure their annihilation or decay products that arrive from outer space. 
Various experiments have searched for antiprotons and positrons at the boundary of 
the Earth's atmosphere since those antiparticles should stick out among the wealth 
of ordinary particles in the cosmic radiation. Gamma-rays as annihilation products 
have the additional advantage that their arrival direction can be correlated with the 
dark matter density, and signals recorded by both space-borne and ground-based 
7-ray detectors have been searched for indications for dark matter. 

2. Antiproton and Positron Searches with Balloon Experiments 

Balloons carry particle spectrometers to heights of around 40 km above sea level 
but the reached exposure is somewhat limited since these experiments (like BESS 
or HEAT) are flown typically once a year for a duration of roughly 10 hours. The 
BESS experiment [4], for example, recorded around 2000 antiprotons in eight flights 
which took place in northern Canada over a period of 11 years. The small fraction 
of antimatter in the cosmic radiation {p/p ~ 10~^) poses a particular challenge to 
the balloon-borne spectrometers, and the large background of ordinary particles 
(i.e. protons and electrons) requires high precision tracking and excellent particle 
identification capabilities. Some of the associated experimental uncertainties were 
maybe not so well known in early experiments whose data suggested a much larger 
number of antiprotons and positrons than expected from purely secondary produc- 
tion. Meanwhile, both the understanding of the measurements and of the expected 
background of secondary antiprotons have improved, resulting in less room for ad- 
ditional sources of antimatter. 
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Fig. 1. Positron fraction measured by HEAT 
in 1994/95 and 2000. The lines indicate a con- 
tribution from annihilating Kaluza-Klein dark 
matter which has been normalized to the data. 
The various line styles correspond to differ- 
ent choices for the positron diffusion parame- 
ters [7]. 
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Fig. 2. SPI 511 keV flux spectrum that is ob- 
tained when modelling the Galactic Centre 
emission region as a Gaussian with a FWHM 
of 10° [16]. 



For antiprotons, both the absolute flux and the energy dependence of the antipro- 
ton-to-proton ratio seem to be broadly consistent with a picture where no dark 
matter contribution is needed to account for the observations. The BESS collab- 
oration reported a possible excess of antiprotons [4] at energies below 1 GeV, but 
the effect seems small given the limited statistics of the data and the uncertainties 
related to the proper description of solar modulation in this energy range. In De- 
cember 2004, the BESS collaboration managed for the first time to conduct a 9-day 
flight in Antarctica that should have doubled the available antiproton statistics and 
might give new insight into the reality of the observed effect. 

Positron measurements were conducted as well and analysed in terms of the 
so-called positron fraction (the fraction of positrons among all detected electrons 
and positrons). The energy dependence of the positron fraction has been measured 
between 1 and 30 GeV, and features of the observed spectrum have been interpreted 
as being due to dark matter. In particular the data of the HEAT experiment (High 
Energy Antimatter Telescope) have been discussed in this context. The HEAT 
spectrometer has been flown in two different detector configurations (HEAT-e* and 
HEAT-pbar). The flights took place near solar maximum (1994 and 1995, HEAT- 
e*) and near solar minimum (2000, HEAT-pbar), and at different geomagnetic 
cutoffs (around 1 GeV in 1995, and around 4 GeV in 1994 and 2000). The energy 
dependence of the positron fraction measured during the 1994/95 flights agrees 
with what has been found in 2000 [5,6] (Fig. 1), and the HEAT collaboration finally 
concluded that the intensity of positrons is consistent with a purely secondary origin 
due to nuclear interactions in interstellar space [6] . 

Given the experimental and model-related uncertainties, a contribution from dark 
matter cannot be ruled out and this possibility has been the subject of theoretical 
studies. A popular dark matter candidate motivated by elementary particle physics 
(the supersymmetric neutralino) couples in proportion to the mass of the final 
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state particles and is thus an inefficient electron-positron source. In this situation, 
particles expected in Kaluza-Klein theories ofter a viable alternative and have been 
studied in detail [7] (Fig. 1). 

3. Electron-Positron Annihilation Radiation from the Galactic Centre 

Another way of tracing positrons is the observation of the 511 keV radiation that 
is created when electrons and positrons annihilate at rest somewhere in the Milky 
Way. Energetic positrons suffer rapid ionization losses and can annihilate with elec- 
trons depending on the atom and electron density in the interstellar medium. Since 
the stopping distance of positrons is typically much smaller than their mean free 
path for annihilation, positrons do indeed thermalise before annihilation. The anni- 
hilation process itself seems to be dominated by positronium formation [8] creating 
a narrow line (25 % of the time) or a three-photon continuum (75 % of the time). 

Electron-positron annihilation radiation from the Galactic Centre was detected 
for the first time in the 70s [9-11]. After initial claims of a flux variability in time, 
consistent results for the absolute photon flux and the line shape were obtained by a 
number of experiments [12-15]. There was much less agreement on the morphology 
of the emission region. Some measurements suggested both a contribution from the 
bulge and the disk of the Milky Way, and an additional component from a location 
at positive galactic latitude (i.e. above the galactic plane) was also discussed. 

The most recent high-statistics observations of the 511 keV annihilation radiation 
have been obtained using the SPI detector onboard the INTEGRAL satellite. SPI 
is sensitive for photon energies between 20keV and lOMeV and has an energy 
resolution of 2 keV (at 1 MeV) . The arrival direction of photons is reconstructed 
with an angular resolution of around 2°. 

The SPI data [16-18] cover galactic longitudes of ±30° and galactic latitudes 
of ±10°. The data show a strong signal (Fig. 2) with a significance exceeding 20 
standard deviations. The measured flux and the line width are consistent with 
earlier measurements. Despite significant exposure out to galactic longitudes of 
±30° the morphology of the emission region is accurately modelled by a spherical 
distribution which is centred at the Galactic Centre and exhibits a Gaussian radial 
profile with a FWHM of « 10°. Upper limits on the flux contribution from the 
galactic disk were derived. 

The restriction of the annihilation signal to the Galactic Centre region seems 
difficult to explain in the framework of conventional astrophysical positron sources. 
Most of the discussed sources (like supernovae, micro-quasars, pulsars, stellar flares, 
Wolf-Rayet stars, cosmic-ray interactions with the interstellar medium and others, 
see for example [19-21]) would probably imply a non-negligible contribution from 
the galactic disk. To avoid such a contribution is has even been speculated that the 
large-scale magnetic field in the galaxy could be such that positrons are channeled 
from the disk into the bulge [20] . 

The SPI measurement suggests the galactic bulge as the positron source and this 
agrees with the notion that the density of dark matter particles should be highest 
in the Galactic Centre region, while a contribution from the galactic disk is not 
necessarily expected. This has resulted in attempts to model the SPI data as being 
due to dark matter annihilation [22] and it has been shown that the right relic 
density can be obtained. To evade observational constraints the dark matter should 
be light (1-100 MeV) to create the e+e^ final state almost exclusively. The exchange 
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Fig. 3. Diffuse 7-ray energy spectrum of tlie 
Inner Galaxy (longitude 330 — 30°, latitude 
|fe| < 5°). The EGRET data are sfiown as dots. 
The dashed line denotes the photon contribu- 
tion from hadronic interactions in the interstel- 
lar gas, the dot-dashed line denotes the pho- 
ton spectrum expected from neutralino annihi- 
lations into bb pairs. The solid line is the sum 
of the two components [24] . 



Fig. 4. Spectral energy density of 7-rays 
from observation of the Galactic Centre with 
H.E.S.S. in 2003 (red circles) and 2004 (black 
dots). The shaded area is the result of a pow- 
er-law fit dN/dE oc . The lines show the 
spectra expected from a MSSM neutralino with 
a mass of 14TeV (dashed), from Kaluza-Klein 
DM with a mass of 5 TeV, and a 10 TeV DM 
particle (solid) annihilating into t+t" (30%) 
and bb (70%). (From [44].) 



particle in the annihilation process could be a fermion (with suppressed Z couplings 
to evade collider constraints) or a new type of gauge boson {U boson). This dark 
matter interpretation has been checked against various experimental constraints, 
and it has been shown for example that already the presence of Bremsstrahlung 
photons in the annihilation process would violate COIVIPTEL data if the mass of 
the dark matter particle exceeds 20MeV [23]. 



4. Diffuse 7-ray Emission from the Milky Way 

Photons from the decay cascade of particles created by annihilating dark matter 
can potentially make a significant contribution to the diffuse 7-ray emission from 
the ]VIilky Way and this possibility has lately been investigated [24] using EGRET 
data. The EGRET (Energetic Gamma Ray Experiment Telescope) instrument was 
a detector sensitive to photons with energies between 20 IVIeV and 30 GeV that was 
flown between 1991 and 2000 aboard NASA's Compton Gamma-Ray Observatory. 
It allowed to reconstruct photon energies with an accuracy of 20%, and its angular 
resolution varied between 1.3° (at 1 GeV) and 0.4° (at 10 GeV). 

The EGRET observations discovered several hundred 7-ray sources. Diffuse emis- 
sion was searched for by subtracting those sources from the overall photon flux, and 
it was noted already in 1997 that the energy spectra of the diffuse component show a 
60 % excess compared with model predictions [25] . This excess was found at energies 
above several GeV, while the EGRET data below 1 GeV were well accounted for by 
models that describe the photon spectrum as entirely due to secondary interactions 
of charged cosmic rays. The excess was found with a very similar spectral shape 
for all sky directions, i.e. it does not only occur when observing, for example, the 
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Galactic Centre region. 

An excess of 7-rays from all directions is a signature that could be provided by 
dark matter annihilations taking place throughout the Milky Way. The shape of the 
energy spectrum of the excess photons is determined by the mass and the couplings 
of the dark matter particle and should indeed be independent of the observation 
direction. An analysis of the total EGRET data set has been performed to investi- 
gate possible dark matter contributions to the diffuse photon flux [24] . The analysis 
extracts the background shape from a model of cosmic ray propagation and 7-ray 
production in the Galaxy (GALPR0P[3]) with a particular sets of parameters, 
while the signal shape is based on expectations from supersymmetry [26] combined 
with the known fragmentation of quarks. The absolute normalizations of signal and 
background are allowed to vary, and a fit of the energy spectra is used to limit 
the mass of the dark matter particle to the 50-100 GeV range (Fig. 3). The spatial 
distribution of the detected photons is used to extract very detailed information 
about the dark matter distribution in the Galaxy. It is found that the EGRET 
data can be described with a contribution from annihilating sypersymmetric dark 
matter when the annihilation rate is enhanced by a 'boost factor' of 20 which could 
be possible due to clumping on small scales. 

The proper modelling of the sizable fraction of background photons created by 
interactions of charged cosmic rays is a key element when looking for dark matter 
induced photons. The GALPROP model combined with cosmic ray measurements 
(i.e. local (anti)proton and electron spectra) has been used to describe the EGRET 
data on diffuse photon emission without any dark matter contribution. This alter- 
native analysis [27] uses the diffuse emission itself to adjust the electron injection 
spectrum and also scales up the initial proton spectum using antiproton data. It ar- 
rives at a different GALPROP parameter set and somewhat relaxes the requirement 
on the agreement of locally measured and predicted electron and proton spectra. 
This is not implausible since details of the galactic structure are anyway not taken 
into account in the background model. Recently, it has been pointed out that the 
description of the EGRET excess with the help of dark matter as proposed by [24] 
would result in an antiproton flux higher than measured by balloon experiments [28] . 

5. Very High-Energy 7-rays from the Galactic Centre 

The Galactic Centre region is considered as a promising target for dark matter 
searches since an increased dark matter density is expected there. The dynami- 
cal centre of the Galaxy is dominated by Sagittarius A*, a putative black hole of 
w 3 X 10^ solar masses [29]. Sgr A* was first discovered as compact radio source and 
exhibits a very low bolometric luminosity. It was found to be variable on sub-hour 
timescales in the infrared band and in X-rays [30-33] suggesting that the observed 
radiation is created close to the event horizon of the black hole. Sgr A* is sur- 
rounded by a supernova remnant (Sgr A East) and a region with a high density 
of ionized hydrogen (Sgr A West). Recently, also an energetic pulsar candidate has 
been discovered close to Sgr A* [43]. 

In 2004, results of observations of the Galactic Centre in very high-energy 7-rays 
(> 100 GeV) have been published by several experiments [34-36] using the imaging 
atmospheric Cherenkov technique. The most significant and precise 7-ray signal was 
reported by the H.E.S.S. experiment above an energy threshold of 160 GeV [36] . The 
detected source appears point-like compared with the point spread function of the 
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instrument. Within errors, the source position is compatible with Sgr A* , but also 
with the supernova remnant Sgr A East and the energetic pulsar candidate. The 
energy spectrum measured by H.E.S.S. is well described by a power law with 
a spectral index F of 2.25 (symbols in Fig. 4), and extends well beyond to 20TeV. 
No indications for flux variability were found in the data on time scales down to 10 
minutes. 

The 7-ray signal from the Galactic Centre has been interpreted both as dark 
matter [37,44] and in terms of 'ordinary' astrophysical processes. When ascribing 
the total photon signal to dark matter annihilations the H.E.S.S. energy spectrum 
implies dark matter masses in the order of lOTeV (Fig. 4) regardless whether the 
dark matter candidate comes supersymmetric or Kaluza-Klein theories. The power- 
law shape of the spectrum does not fit to the cutoff that is in general expected from 
the annihilation spectra of dark matter particles. It has also been tried to search 
for a dark matter annihilation spectrum that is hidden under an astrophysical 
background with a power-law shape, but no significant dark matter component was 
found [44] . 

Astrophysical explanations of the Galactic Centre signal have been studied as 
well and could be viable alternatives. Accretion onto the black hole Sgr A* and as- 
sociated acceleration processes can produce high-energy protons and electrons that 
create secondary 7-rays due to synchrotron and curvature radiation or due to inter- 
actions with ambient matter and photon fields. The low luminosity of Sgr A* implies 
that it appears transparent for TeV 7-rays, and it has been demonstrated [39] that 
the H.E.S.S. data can be explained assuming populations of high-energy protons or 
electrons associated with Sgr A* . On the other hand, the lack of temporal variability 
makes an association of the TeV 7-ray signal with Sgr A* less likely. 

Other astrophysical interpretations focus on the supernova remnant Sgr A East [40] . 
The spectral index measured by H.E.S.S. is close to the value expected for par- 
ticle acceleration in the shock front of supernova explosions. The observed 7-rays 
would then be created by protons and/or electrons accelerated to energies of around 
100 TeV. This explanation appears not unmotivated since Sgr A East is a very 
powerful supernova remnant with an above-average explosion energy, and the pres- 
ence of a 100 TeV-particles has lately been demonstrated for other supernova rem- 
nants [41,42]. It has also been argued that the energetic young pulsar would be 
luminous enough to explain the TeV 7-radiation. 

6. Summary 

A number of signals recorded by balloon experiments, satellites and imaging 
Cherenkov telescopes have been interpreted as being due to dark matter. Dark 
matter explanations, often of the same signal, invoke various theoretical frame- 
works, and there seems to be no scenario that can explain more than experimental 
signature. In most cases, astrophysical sources of the observed signals offer a rea- 
sonable alternative. 
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Abstract 

Astrophysical bounds on the properties and abundances of primordial quark nuggets and 
cosmic ray strangelets are reviewed. New experiments to search for cosmic ray strangelets 
in lunar soil and from the International Space Station are described. Analogies with bary- 
onic and supersymmetric Q-balls are briefly mentioned, as are prospects for strangelets as 
ultra-high energy cosmic rays. 



1. Introduction 

Quark nuggets, nuclearites and strangelets are different names for lumps of a 
hypothetical phase of absolutely stable quark matter, so-called strange quark matter 
because of the admixture of slightly less than one-third strange quarks with the up 
and down quarks. Whether strange quark matter is absolutely stable is a question 
yet to be decided by experiment or astrophysical observation (see [1,2] for reviews), 
but if it is the case, then strange quark matter objects may exist with baryon 
numbers ranging from ordinary nuclei to a maximum of order 2 x 10^^ corresponding 
to gravitational instability of strange stars. 

Truly macroscopic c^uark matter lumps surviving from the cosmological quark- 
hadron phase transition are often referred to as quark nuggets, and if they hit the 
Earth they are sometimes dubbed nuclearites. Strangelets are smaller lumps (baryon 
number A < 10^) where the electron cloud neutralizing the slightly positive quark 
charge mainly resides outside the quark core. Strangelets are unlikely to survive 
from the early Universe, but may form as a result of strange star binary collisions 
and/or acceleration from the surface of pulsars. The nomenclature is not strictly 
defined, and in the following the word strangelet will be used as a general name 
except when discussing leftovers from the early Universe. 

Q-balls are non-topological solitons suggested from various origins in the early 
Universe and as such have nothing to do with strange quark matter (their origin 
and general properties were described by Kuscnko at this workshop). However for 
some classes of Q-balls (baryonic and supersymmetric Q-balls), the astrophysical 
bounds that can be derived on strangelets arc easily generalized to these creatures 
as well and are therefore of interest in the context of this Workshop on Extreme 
Physics with Neutrino Telescopes. 

In the following I shall discuss the (unlikely) survival of cosmological quark 
nuggets, the more optimistic prospects for cosmic rays strangelets from strange 
stars, and two new experimental efforts to search for them. Due to space limita- 
tions I will not go into details with the Q-ball analogies but refer the reader to 
Kusenko's contribution and to [3]. 
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2. Sources, sizes, and fluxes 

2.1. Primordial quark nuggets 

In a first order cosmological quark-hadron phase transition at T « 100 MeV, 
supercooling may lead to concentration of baryon number inside shrinking bubbles 
of quark phase, that may reach nuclear matter density and form quark nuggets. The 
baryon number inside the horizon during the cosmic quark-hadron phase transition 
(an upper limit for causal formation of quark nuggets) is Ahor ~ lO''^. Witten [4] 
predicted that typical nuggets would be somewhat smaller than this and argued that 
quark nuggets might explain the cosmological dark matter problem. Quark nuggets 
would decouple from the radiation bath very early in the history of the Universe and 
behave as cold dark matter in the context of galaxy formation. Today the nuggets 
would move with typical galactic halo velocities of a few hundred kilometers per 
second through the Milky Way. 

Later studies showed that the hot environment made cosmological nuggets un- 
stable against surface evaporation [5] and boiling [6] , effectively destroying nuggets 
with A below 10'^^~^^ depending on assumptions. Small traces of primordial nuggets 
with lower baryon numbers could also be left over from the destruction processes 
in the early Universe. Even such traces may in fact be "observed" using the astro- 
physical detectors discussed below. Let v = 250km s^^V250 and p = 10~^'*gcm~^p24 
be the typical speed and mass density of nuggets in the galactic halo. The speed 
is given by the depth of the gravitational potential of our galaxy, whereas p24 ~ 1 
corresponds to the density of dark matter. Then the number of nuggets hitting the 
Earth is 

jc- « 6.0 X 10^^-1/924^250 cm-2 sT-\ (1) 

Quark nuggets have a positive electrostatic surface potential (several MeV) of 
the quark phase because quarks are stronger bound than electrons, so during Big 
Bang nucleosynthesis (T < 1 MeV), nuggets absorb neutrons but not protons. This 
reduces the neutron-to-proton ratio, thereby lowering the production of ''He. The 
helium- production is very sensitive to the total amount of nugget- area present, and 
in order not to ruin the concordance with observations, one finds [7] that only 
nuggets with A > Abbn ~ ^^^^^nugfn ^^'^ allowed during nucleosynthesis. Here 
J7nug is the present-day nugget contribution to the cosmic density in units of the 
critical density, and /n < 1 is the penetrability of the nugget surface. 

In spite of carrying baryon number, primordial quark nuggets do not contribute 
to the usual nucleosynthesis limit on llbaryon- The baryon number is "hidden" in 
quark nuggets long before Big Bang nucleosynthesis begins, and the nuggets only 
influence nucleosynthesis via the neutron absorption just described. The same is 
true for baryon number carrying Q-balls. 

While primordial quark nuggets remain a possibility within the tight restrictions 
mentioned, the main problem with this scenario is the need for a first order quark- 
hadron phase transition which is currently not favored in lattice QCD studies at 
zero chemical potential. 

2.2. Strangelets from compact stars 

If strange quark matter is absolutely stable all compact stars are likely to be 
strange stars (see the following Section), and therefore the galactic coalescence 
rate estimated for neutron star binaries that inspiral due to loss of orbital energy 
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by emission of gravitational radiation, believed to be of order one collision in our 
Galaxy every 10,000 years [8], is really the rate of strange star collisions. Each event 
involves a phase of tidal disruption as the stars approach each other before the final 
collision. During this stage small fractions of the total mass may be released from 
the binary system in the form of strange quark matter. No realistic simulation of 
collisions involving two strange stars has been performed. Simulations of binary 
neutron star collisions, depending on orbital and other parameters, lead to the 
release of anywhere from 10~^ — 1O~^M0 {Mq is the solar mass), corresponding 
to a total mass release in the Galaxy of 10"^" — 1O~^M0 per year. The equation 
of state for strange quark matter is stiff, so strange star collisions probably lie 
in the low end of the mass release range. A conservative estimate of the galactic 
production rate of strangelets is lO~^''M0yr~^. 

Quark matter lumps released by tidal forces are macroscopic [9] , but subsequent 
collisions lead to fragmentation, and if the collision energy compensates for the sur- 
face energy involved in making smaller strangelets, a significant fraction of the mass 
released from binary strange star collisions might end up in the form of strangelets 
with A « 10^ — 10^ [9]. Incidentally, a similar range of strangelet masses is expected 
if the surface of strange stars consists of a layer with strangelets embedded in an 
electron gas rather than pure quark matter all the way to the surface [10]. 

Assuming that strangelets from binary collisions are accelerated and propagate 
like cosmic ray nuclei in our Galaxy, taking proper account of their small charge- 
to-mass ratio, as well as energy loss, spallation, escape from the Galaxy, etc., it 
was shown in [11] that the expected flux for color-flavor locked strangelets (charge 
Z = 0.3A2/3 [12]) near Earth is 

lO-^^-i-^W^s-^sr-i. (2) 

Most of these nuggets have rigidities (momentum divided by charge) of a few GV, 
but with a powerlaw tail at higher rigidity. Apart from the slightly different A- 
dependence this is some 12 orders of magnitude smaller than the flux estimate 
for dark matter nuggets, which is not unreasonable because the total strangelet 
mass originating from binary collisions over the age of the Galaxy is around 1 Mq, 
compared to 10^^ Mq of dark matter. 

Another possible cosmic ray strangelet source is extraction from the surface of 
pulsars and acceleration in the strong pulsar electric fields. A measurable flux is 
predicted in [13] within the scenario where the strange star surface consists of 
strangelets embedded in an electron gas [10]. Formation and acceleration in super- 
nova explosions has also been suggested [14]. 

3. Detection of cosmic ray strangelets (or why either all or no compact 
stars are strange) 

De Riijula and Glashow [15] argued that quark nuggets hitting the Earth would 
show up as unusual meteor-events, earth-quakes, etched tracks in old mica, in me- 
teorites and in cosmic-ray detectors. A negative search for tracks in ancient mica 
corresponded to a lower nugget flux limit of 8 x 10~^^ cm~^ s~^ sr~^ for nuggets 
with A > 1.4 X 10^^ (smaller nuggets being trapped in layers above the mica samples 
studied). 

Later investigations have improved these flux limits by a few orders of magnitude 
and extended them to lower A, though with higher flux limits (see other contribu- 
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tions to these proceedings for examples). This has excluded quark nuggets with 
3 X 10^ < ^ < 5 X 10^^ as dark matter, but a low flux from the Big Bang or from 
collision of strange stars cannot be ruled out. 

Neutron stars and their stellar progenitors may be thought of as alternative large 
surface area, long integration time detectors leading to much tighter flux limits 
[16], see also [17]. The presence of a single quark nugget in the interior of a neutron 
star is sufficient to initiate a transformation of the star into a strange star [4,18]. 
The time-scale for the transformation is short, between seconds and minutes, so 
observed pulsars would have been converted long ago if their stellar progenitors 
ever captured a quark nugget, or if neutron stars absorbed one after formation. 

To convert a neutron star into strange matter a quark nugget should not only 
hit a supernova progenitor but also be caught in the core [16]. A main sequence 
star is capable of capturing non-relativistic quark nuggets with baryon numbers 
below AgTOP ~ lO'^^, which are braked by inertia, i. e. they are slowed down by 
electrostatic scatterings after plowing through a column of mass similar to their 
own, and settle in the stellar core. Relativistic nuggets may be destroyed after 
collisions with nuclei in the star, but even a tiny fraction of a nugget surviving such 
an event and settling in the star is sufficient to convert the neutron star to a strange 
star, so the non-relativistic flux limits may still apply. 

For nuggets with A < ^stop the sensitivity of main sequence stars as detectors 
is given by the limit of one nugget hitting the surface of the supernova progenitor 
in its main sequence lifetime. Converted into a flux, JF, of nuggets hitting the Earth 
it corresponds to 

« 4 X 10"''2w|5oCm-2 s"i ST-\ (3) 

This is a factor of 10^"-10^° more sensitive than direct detection experiments! 

// it is possible to prove that some neutron stars are indeed neutron stars rather 
than strange stars, the sensitivity of the astrophysical detectors rules out quark 
nuggets as dark matter for A < 10'^^"'^*. And it questions the whole idea of sta- 
ble strange quark matter, since it is impossible to avoid polluting the interstellar 
medium with nuggets from strange star collisions or supernova explosions at fluxes 
many orders of magnitude above the limit measurable in this way. 

// on the other hand strange quark matter is stable, then all neutron stars are 
likely to be strange stars, again because some pollution can not be avoided. 

4. Experiments underway 

Several experiments have searched for strangelets in cosmic rays. While some 
interesting events have been found that are consistent with the predictions for 
strangelets, none of these have been claimed as real discoveries. Interpreted as flux 
limits rather than detections these results are consistent with the flux estimates 
given above. For discussions see [19,20]. 

If the interesting events were actual measurements, two experiments that are 
currently underway will reach sensitivities that would provide real statistics. 

AMS-02: The Alpha Magnetic Spectrometer (AMS) is a space-based particle 
physics experiment involving several hundred physicists from more than 50 insti- 
tutions in 16 countries, led by Nobel laureate Samuel Ting of MIT. A prototype 
(AMS-01) flew in June 1998 aboard the Space Shuttle Discovery [21], and AMS- 
02 is currently scheduled to fly to the International Space Station (ISS) in 2009. 
Once on the ISS AMS-02 will remain active for at least three years. AMS-02 will 
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provide data with unprecedented accuracy on cosmic ray electrons, positrons, pro- 
tons, nuclei, anti-nuclei and gammas in the GV-TV range and probe issues such as 
antimatter, dark matter, cosmic ray formation and propagation. In addition it will 
be uniquely suited to discover strangelets characterized by extreme rigidities for a 
given velocity compared to nuclei [19,20]. AMS-02 will have excellent charge reso- 
lution up to Z w 30, and should be able to probe a large mass range for strangelets. 
A reanalysis of data from the AMS-01 mission has given hints of some interesting 
events, such as one with Z = 2, A = 16 [22] and another with Z = 8, but with the 
larger AMS-02 detector running for 3 years or more, real statistics is achievable. 

LSSS: The Lunar Soil Strangelet Search (LSSS) is a search for Z = 8 strangelets 
using the tandem accelerator at the Wright Nuclear Structure Laboratory at Yale 
[20,23]. The experiment involves a dozen people from Yale, MIT, and Arhus, led by 
Jack Sandweiss of Yale. The experiment which is about to begin its real data taking 
phase, studies a sample of 15 grams of lunar soil from Apollo 11. It will reach a 
sensitivity of 10~^^ over a wide mass range, sufficient to provide detection according 
to Eq. 2 if strangelets have been trapped in the lunar surface layer, which has an 
effective cosmic ray exposure time of around 500 million years and an effective 
mixing depth due to micrometeorite impacts of only around one meter, in contrast 
to the deep geological and oceanic mixing on Earth. Combined with the fact that 
the Moon has no shielding magnetic field, this results in an expected strangelet 
concentration in lunar soil which is at least four orders of magnitude larger than 
the corresponding concentration on Earth. 

5. Strangelets as ultra-high energy cosmic rays 

The existence of cosmic rays with energies well beyond lO^^eV, with measured 
energies as high as 3 x lO^^eV, is one of the most interesting puzzles in cosmic 
ray physics [24]. It is almost impossible to find a mechanism to accelerate cosmic 
rays to these energies. Furthermore ultra-high energy cosmic rays lose energy in in- 
teractions with cosmic microwave background photons, and only cosmic rays from 
nearby (unidentified) sources would be able to reach us with the energies mea- 
sured. Strangelets circumvent both problems [25] , and therefore provide a possible 
mechanism for cosmic rays beyond the so-called GZK-cutoff. 

Acceleration: All astrophysical "accelerators" involve electromagnetic fields, 
and the maximal energy of a charged particle is proportional to its charge. The 
charge of massive strangelets has no upper bound in contrast to nuclei, so highly 
charged strangelets are capable of reaching energies much higher than those of 
cosmic ray protons or nuclei using the same "accelerator" [25]. 

The GZK-cutofT is a consequence of ultra-relativistic cosmic rays hitting a 2.7K 
background photon with a Lorentz-factor 7 large enough to boost the 7 x 10~^ cV 
photon to energies beyond the threshold of energy loss processes, such as photo- 
pion production or photo-disintegration. The threshold for such a process has a 
fixed energy, E'Thr, in the frame of the cosmic ray, e.g., i^xhr ~ lOMeV for photo- 
disintegration of a nucleus or a strangelet, corresponding to 7Thr = E^hi/ E2.7K ~ 
10^°, or a cosmic ray total energy 

-Exotai = 7Thr Amoc^ « 10^^ A eV. (4) 

Since strangelets can have much higher A- values than nuclei, this pushes the GZK- 
cutoff energy well beyond the current observational limits for ultra-high energy 
cosmic rays [25,26]. 
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6. Conclusion 

Lumps of strange quark matter (quark nuggets, nuclearites, strangelets) may 
form in a first-order cosmological quark- hadron phase transition (unlikely), or in 
processes related to compact stars (more likely). Flux estimates for lumps reaching 
our neighborhood of the Galaxy as cosmic rays are in a range that makes it realistic 
to either detect them in upcoming experiments like AMS-02 or LSSS, or place severe 
limits on the existence of stable strange quark matter. A similar line of reasoning 
applies to Q-balls. 
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Supersymmetric extensions of the Standard Model predict the existence of Q-balls with 
baryon and lepton numbers. Stable Q-balls can form at the end of inflation from the 
fragmentation of the Affleck-Dine condensate and can exist as dark matter. The best 
current limits come from Super-Kamiokande and MACRO. The search beyond these limits 
can be conducted using the future water Cherenkov detectors. 



1. Introduction 

Supersymmetry is a theoretically appealing possibility for physics beyond the 
Standard Model. It predicts the existence of new particles and extended objects, Q- 
balls [1,2]. Both the Q-balls and the lightest supersymmetric particle are candidates 
for dark matter. Here we will briefly review the former possibility. 

2. Q-balls from Supersymmetry 

In a class of theories with interacting scalar fields cf) that carry some conserved 
global charge, the ground state is a Q-ball [3,4], a lump of coherent scalar condensate 
that can be described scmiclassically as a non-topological soliton of the form 



Q-balls exist whenever the scalar potential satisfies certain conditions that were first 
derived for a single scalar degree of freedom [3] with some abclian global charge and 
were later generalized to a theory of many scalar fields with different charges [1]. 
Non-abclian global symmetries [5] and abelian local symmetries [6] can also yield 
Q-balls. 

For a simple example, let us consider a field theory with a scalar potential U (if) 
that has a global minimum U{0) = at ip = Q. Let U{cp) have an unbroken global^ 
U(l) symmetry at the origin, tp ~ 0. And let the scalar field cp have a unit charge 
with respect to this U(l). 

The charge of some field configuration ip{x,t) is 



^ Q-balls associated with a local symmetry have been constructed [6]. An important qualitative 
difference is that, in the case of a local symmetry, there is an upper limit on the charge of a stable 
Q-ball. 
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Since a trivial configuration ip{x) = has zero charge, the solution that minimizes 
the energy, 



(3) 



and has a given charge Q > 0, must differ from zero in some (finite) domain. 
This is a Q-ball. It is a time-dependent solution, more precisely it has a time- 
dependent phase. However, all physical quantities are time- independent. Of course, 
we have not proved that such a "lump" is finite, or that it has a lesser energy than 
the collection of free particles with the same charge; neither is true for a general 
potential. A finite-size Q-ball is a minimum of energy and is stable with respect to 
decay into free iy9-particles if 

U{lp) /ip^ = min, for = LpQ> Q. (4) 

One can show that the equations of motion for a Q-ball in 3+1 dimensions are 
equivalent to those for the bounce associated with tunneling in 3 Euclidean dimen- 
sions in an effective potential Uuj{ip) = U{(f) — (l/2)ti;^(/3^, where oj is such that it 
extremizes [8] 

Ec. = ^3(^)+wQ. (5) 

Here 5*3(0;) is the three-dimensional Euclidean action of the bounce in the potential 
U^{ip). The Q-ball solution has the form (1), where (p{x) is the bounce. 

The analogy with tunneling clarifies the meaning of condition (4), which simply 
requires that there exist a value of ui, for which tJ^{tp) is negative for some value of 

= </5o 7^ separated from the false vacuum by a barrier. This condition ensures 
the existence of a bounce. (Clearly, the bounce does not exist if Uuj{v) > for all 
ip because there is nowhere to tunnel.) 

In the true vacuum, there is a minimal value too, so that only for to > loq, Ut^{ip) 
is somewhere negative. If one considers a Q-ball in a metastable false vacuum, then 
ujQ = 0. The mass of the ip particle is the upper bound on lu in cither case. Large 
values of to correspond to small charges [8]. As Q —^ 00, lo —^ loq. In this case, the 
effective potential U^{lp) has two nearly-degenerate minima; and one can apply the 
thin- wall approximation to calculate the Q-ball energy [3]. For smaller charges, the 
thin- wall approximation breaks down, and one has to resort to other methods [8]. 

The above discussion can be generalized to the case of several fields, with 
different charges, qt [1]. Then the Q-ball is a solution of the form 

Vk{x,t) = e''"^^'pk{x), (6) 

where (p{x) is again a three-dimensional bounce associated with tunneling in the 
potential 

U^{p) = U{p) - (7) 

k 

As before, the value of w is found by minimizing E^^ in equation (5). The bounce, 
and, therefore, the Q-ball, exists if 



for|'^ol'>0- (8) 
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The soliton mass can be calculated by extremizing E„ in equation (5). If lygp 
defined by equation (8) is finite, then the mass of a soliton M{Q) is proportional 
to the first power of Q: 

M{Q)^fiQ, if|(polVoo- (9) 

In particular, if Q — > oo, ji ^ ^ (thin- wall limit) [3,4]. For smaller values of Q, 
ji was computed in [8]. In any case, jl is less than the mass of the particle by 
definition (8). 

However, if the scalar potential grows slower than the second power of (j), then 
|</5oP = oo, and the Q-ball never reaches the thin-wall regime, even if Q is large. 
The value of cj) inside the soliton extends as far as the gradient terms allow, and the 
mass of a Q-ball is proportional to , p < 1. In particular, if the scalar potential 
has a flat plateau U{(j)) ^ m at large <j>, then the mass of a Q-ball is [12] 

M(Q) ~ mg3/4. (10) 

This is the case for the stable baryonic Q-balls in the MSSM discussed below. 

It turns out that all phenomenologically viable supersymmctric extensions of 
the Standard Model predict the existence of non-topological solitons [1] associated 
with the conservation of baryon and lepton number. If the physics beyond the 
standard model reveals some additional global symmetries, this will further enrich 
the spectrum of Q-balls [7]. The MSSM admits a large number of different Q- 
balls, characterized by (i) the quantum numbers of the fields that form a spatially- 
inhomogeneous ground state and (ii) the net global charge of this state. 

First, there is a class of Q-balls associated with the tri-linear interactions that 
are inevitably present in the MSSM [1]. The masses of such Q-balls grow linearly 
with their global charge, which can be an arbitrary integer number [8]. Baryonic 
and leptonic Q-balls of this variety are, in general, unstable with respect to their 
decay into fermions. However, they could form in the early universe through the 
accretion of global charge [9,10] or, possibly, in a first-order phase transition [11]. 

The second class [12] of solitons comprises the Q-balls whose VEVs are aligned 
with some flat directions of the MSSM. The scalar field inside such a Q-ball is a 
gauge-singlet [13] combination of squarks and sleptons with a non-zero baryon or 
lepton number. The potential along a flat direction is lifted by some soft supersy- 
mmetry-breaking terms that originate in a "hidden sector" of the theory at some 
scale Ag and are communicated to the observable sector by some interaction with 
a coupling g, so that gA ^ 100 GeV. Depending on the strength of the mediating 
interaction, the scale A^, can be as low as a few TeV (as in the case of gauge- 
mediated SUSY breaking), or it can be some intermediate scale if the mediating 
interaction is weaker (for instance, g ~ /"'^pia^ck ^^'^ ~ ^^^'^ GeV in the 
case of gravity-mediated SUSY breaking). For the lack of a definitive scenario, one 
can regard A^ as a free parameter. Below A^ the mass terms are generated for all 
the scalar degrees of freedom, including those that parameterize the flat direction. 
At the energy scales larger than A^, the mass terms turn off and the potential is 
"flat" up to some logarithmic corrections. If the Q-ball VEV extends beyond A^, 
the mass of a soliton [12,14] is no longer proportional to its global charge Q, but 
rather to Q^^^. A hybrid of the two types is yet another possibility [15]. 

This allows for the existence of some entirely stable Q-balls with a large baryon 
number B (B-balls). Indeed, if the mass of a B-baU is Af^ - (1 TeV) x B^/^, then 
the energy per baryon number {M^/B) ~ (1 TeV) x B^^^'^ is less than 1 GeV for 
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B > 10^^. Such large B-balls cannot dissociate into protons and neutrons and are 
entirely stable thanks to the conservation of energy and the baryon number. If they 
were produced in the early universe, they would exist at present as a form of dark 
matter [14]. 



3. Fragmentation of Affleck Dine Condensate into Q-balls 

Several mechanisms could lead to formation of B-balls and L-balls in the early 
universe. First, they can be produced in the course of a phase transition [11]. Sec- 
ond, thermal fluctuations of a baryonic and leptonic charge can, under some con- 
ditions, form a Q-ball. Finally, a process of a gradual charge accretion, similar to 
nucleosynthesis, can take place [9,10,16]. However, it seems that the only process 
that can lead to a copious production of very large, and, hence, stable, B-balls, is 
fragmentation of the Affleck-Dine condensate [14]. 

At the end of inflation, the scalar fields of the MSSM develop some large ex- 
pectation values along the flat directions, some of which have a non-zero baryon 
number [17]. Initially, the scalar condensate has the form given in eq. (1) with 
4>{x) = const over the length scales greater than a horizon size. One can think 
of it as a universe filled with Q-matter. The relaxation of this condensate to the 
potential minimum is the basis of the Affleck-Dine (AD) scenario for baryogenesis. 

It was often assumed that the condensate remains spatially homogeneous from 
the time of formation until its decay into the matter baryons. This assumption 
is, in general, incorrect. In fact, the initially homogeneous condensate can become 
unstable [14] and break up into Q-balls whose size is determined by the potential 
and the rate of expansion of the Universe. B-balls with 12 < log^Q B < 30 can form 
naturally from the breakdown of the AD condensate. These are entirely stable if 
the flat direction is "sufficiently flat" , that is if the potential grows slower than 0^ 
on the scales or the order of 0(0). The evolution of the primordial condensate can 
be summarized as follows: 

Affleck-Dine condensate 



baryonic Q-balls 



unstable stable 
(decay) 



related _ ^ ^ . . 

^ Dark Matter 



This process has been analyzed analytically [14,26] in the linear approximation. 
Recently, some impressive numerical simulations of Q-ball formation have been 
performed [27]; they confirm that the fragmentation of the condensate into Q-balls 
occurs in some Affleck-Dine models. The global charges of Q-balls that form this 
way are model dependent. The subsequent collisions [14,28] can further modify the 
distribution of soliton sizes. 
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In supcrsymmctric extensions of the Standard Model, Q-ball formation occurs 
along flat directions of a certain type, which appear to be generic in the MSSM [29]. 

4. SUSY Q-balls as Dark Matter 

Conceivably, the cold dark matter in the Universe can be made up entirely of 
SUSY Q-balls. Since the baryonic matter and the dark matter share the same origin 
in this scenario, their contributions to the mass density of the Universe are related. 
Most of dark-matter scenarios offer no explanation as to why the observations 
find ^oARK ^ within an order of magnitude. This fact is extremely difficult 
to explain in models that invoke a dark-matter candidate whose present-day abun- 
dance is determined by the process of freeze-out, independent of baryogenesis. If one 
doesn't want to accept this equality as fortuitous, one is forced to hypothesize some 
ad hoc symmetries [30] that could relate the two quantities. In the MSSM with AD 
baryogenesis, the amounts of dark-matter Q-balls and the ordinary matter baryons 
are related [14]; one predicts [18] ^oark ^ lOfis for B-balls with B ^ 10^''. How- 
ever, the size of Q-balls depends on the supersymmetry breaking terms that lift the 
flat direction. The required size is in the middle of the range of Q-ball sizes that can 
form in the Affleck-Dine scenario [14,26,27]. Diffusion effects may force the Q-balls 
sizes to be somewhat smaller, B ~ 10^^ — lO^'* [19]. 

The value B ^ 10^^ is well within the present experimental limits on the baryon 
number of an average relic B-ball, under the assumption that all or most of cold 
dark matter is made up of Q-balls. On their passage through matter, the electrically 
neutral baryonic SUSY Q-balls can cause a proton decay, while the electrically 
charged B-balls produce massive ionization. Although the condensate inside a Q- 
ball is electrically neutral [13], it may pick up some electric charge through its 
interaction with matter [20]. Regardless of its ability to retain electric charge, the 
Q-ball would produce a straight track in a detector and would release the energy 
of, roughly, 10 GeV/mm. The present limits [20,31-33,35] constrain the baryon 
number of a relic dark- matter B-ball to be greater than 10^^. Future experiments 
are expected to improve these limits. It would take a detector with the area of 
several square kilometers to cover the entire interesting range B ~ 10^^. ..10"^". 

5. Interactions with matter, experimental and astrophysical bounds 

A Dirac fermion scattering off a Q-ball can convert into an antifermion, as long 
as the fermion number is spontaneously broken inside the Q-ball by the scalar 
condensate [22]. The baryonic Q-balls can, therefore, interact with matter nuclei 
converting them into pions whose decay can produce a signal in various detectors, 
including Super-Kamiokande [20,32-34]. 

Dark-matter superballs pass through the ordinary stars and planets with a neg- 
ligible change in their velocity. However, Q-balls can stop in the neutron stars and 
accumulate there [21]. As soon as the first Q-ball is captured by a neutron star, 
it sinks to the center and begins to absorb the baryons into the condensate. High 
baryon density inside a neutron star makes this absorption very efficient, and the 
B-ball grows at the rate that increases with time due to the gradual increase in 
the surface area. After some time, the additional dark-matter Q-balls that fall onto 
the neutron star make only a negligible contribution to the growth of the central 
Q-ball [21]. So, the fate of the neutron star is sealed when it captures the first 
superball. 
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Baryonic Q-balls can destroy neutron stars by stimulated nucleon decay in nuclear 
matter [21-23]. Neutron stars are stable in some range of masses. In particular, there 
is a minimal mass (about 0.18 solar mass), below which the force of gravity is not 
strong enough to prevent the neutrons from decaying into protons and electrons. 
While the star is being consumed by a superball, its mass gradually decreases, 
reaching the critical value eventually. Neutron stars are known to exist for as long 
as 10 Gyr, which sets a bound on the rate of absorption by SUSY Q-balls. 

The astrophysical limits on SUSY Q-balls depend on the nature of the flat di- 
rection and the non-renormalizable operators that "lift" it for a sufhciently large 
VEV. The generic lifting terms can be written in the form 

F("'"'(</')nfti„g«A_M4(^^j (^l^j , (11) 

where M is the characteristic high scale, presumably, of the order of the Planck 
scale. The possible lifting terms [24] may preserve the baryon number (m = 0), or 
they may cause a violation of the baryon number for a large scalar VEV (m ^ 0). 
If the leading lifting terms in eq. (11) have to = 0, the baryon number is conserved 
even for very large values of the VEV. The flat directions of this kind generate Q- 
balls that can grow rapidly inside a neutron star, and stringent constraints exists on 
such relic Q-balls [23]. In contrast, if the flat direction is lifted by terms with m ^ 
0, the corresponding Q-balls cannot grow beyond certain size, and the limits from 
neutron stars or white dwarfs do not apply [23]. In the MSSM, the flat directions 
which are not constrained by the stability of neutron stars include QLe, QLd, Lude, 
QLde, QLud, QLude, in the notation of Ref. [24]. 

6. Conclusion 

Supersymmetric models of physics beyond the weak scale offer two plausible 
candidates for cold dark matter: the lightest supersymmetric particle and a stable 
non-topological soliton, or Q-ball, carrying some baryonic charge. 

SUSY Q-balls make an appealing dark-matter candidate because their formation 
is a natural outcome of the Affleck-Dine baryogenesis. The basic assumptions are 
supersymmetry and inflation. The search beyond the current limits can be con- 
ducted using future water Cherenkov detectors. 

This work was supported in part by the DOE grant DE-FG03-91ER40662 and 
the NASA ATP grants NAG 5-10842 and NAG 5-13399. 
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Abstract 

A search for Q-balls has been carried out in Super-Kamiokande II with 541.7 days of live 
time. A neutral Q-ball can interact with nuclei to produce pions, generating a signal of 
successive contained pion events along a track. No candidate for successive contained event 
groups has been found, and upper limits on the flux of such Q-balls are obtained. 



1. Introduction 

In the framework of some supersymmetric models, stable non-topological solitons 
called Q-balls [1] can be produced in the early universe and contribute to the dark 
matter [2-5] . If Q-balls are formed in the early universe and have survived until the 
present time, they would be part of the dark matter located in the halo of Galaxy. 
Q-balls are solitons with theoretically appealing roles in some dark matter scenarios 
and in explanations of the baryon asymmetry [6] . 

In this model [3] , the Q-ball mass and radius are given by 

Mq = ^^MsQK Rq = -^A/5"'Q^ (1) 

where Ms is the energy scale of SUSY breaking and Q is the Q-ball's baryon 
number. 

Q-balls are classified into two groups according to the properties of their interac- 
tions with matter: supersymmetric electrically charged solitons and supersymmetric 
electrically neutral solitons. Scintillator detectors may be more suitable for the de- 
tection of charged Q-balls with large energy loss in matter. Thus, we will confine 
our attention here to neutral Q-balls. 

When a neutral Q-ball collides with a nuclcon it absorbs its baryonic charge and 
induces the dissociation of the nuclcon into free quarks. In this process, about 1 GeV 
energy is released by the emission of typically two or three pions [7,8]. 

The cross section of the interaction between neutral Q-balls and the matter is 
roughly estimated to be the geometrical size of the Q-ball, 

a = ttRq^ (2) 

For 100 mb cross section, the energy loss of Q-ball is estimated to be about 60 MeV/cm. 
The mass of a Q-ball with this cross section is about lO^^McV/c^. Therefore the 
effect of the energy loss in the Earth is expected to be negligible compared to the 
kinetic energy of a Q-ball. 
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2. Super-Kamiokande detector 

Supcr-Kamiokandc is a cylindrical 50kton water Chcrenkov detector located at 
a depth of 2,700 m water equivalent. The water tank is optically separated into two 
regions; the inner detector (ID) instrumented with inward facing 20 inch diame- 
ter photomutiplier tubes (PMT) and the outer detector (OD) instrumented with 
outward facing 8 inch PMTs. Super-Kamiokande II (SK-II) is the second phase of 
Super-Kamiokande experiment and started in December 2002 with 5,182 PMTs in 
the ID and 1,885 PMTs in the OD. See [12] for more details on the detector. The 
inner PMTs are instrumented with acrylic and FRP covers to avoid a chain reaction 
implosion. The OD is used to veto entering cosmic ray muons and to tag exiting 
charged particles. 

In this paper we report the results from a search for neutral Q-balls using 
541.71ive days of data for the SK-II period. A Q-ball interacting in the detec- 
tor will produce several pions which emit Cherenkov light. Thus, we can observe 
these pion events along the Q-ball track in the detector. 

3. Simulation of Q-ball events 

A Monte Carlo simulation is developed to estimate the detection efhciency of 
Q-ball events. The direction of a Q-ball is randomly generated, and the distance 
between the center of the detector and a Q-ball track is assumed to be uniformly 
distributed within 50 m. The Q-ball velocity is generated assuming a Maxwellian 
distribution with a mean value of 270 km/sec — a typical galactic velocity. Therefore, 
the velocity of Q-balls observed on the earth is determined by superposing the 
velocity of solar system 220 km/sec on the Q-ball velocity in the halo [9]. The pions 
produced in the interactions with a Q-ball and a nucleon can be determined using 
KNO scaling [10]. The total energy of the released pions is equal to the nucleon 
mass. The charges of the pions are randomly generated within the constraint of the 
charge conservation. The kinematics of generated pions are determined by assuming 
the Q-ball and nucleus scatter into the Q-ball and the several pions. 

4. Data Reduction 

In this work we look for a Q-ball signature consisting of at least two successive 
events in which 1 GeV of energy is released in the form of pions in the ID. A Q-ball 
with (7=10 mb could interact about 18 times within 100 /^scc, since the interaction 
length of a Q-ball with P~1Q~^ is 30 m. The details of the event reconstruction in 
Super-Kamiokande can be found in [11]. 

The first selection criteria are: 

1-1) The number of total photoelectrons (p.e.'s) in the ID should be larger than 
300 (the top of Fig. 1) and the maximum number of p.e.'s in the ID within sliding 
300nsec time window should be less than 10,000 (the bottom of Fig. 1), which 
corresponds to about 60MeV and 2 GeV energy for electrons, respectively. These 
cuts reject electron events from decays of cosmic ray muons, high energy through- 
going muon events and electrical noise events. 

1-2) At least one delayed event should exist within a 100/xsec window after the 
event selected by 1-1). This cut rejects single events. 

All events selected by 1-1) and 1-2) are considered to be an event group. The 
event selected by criterion 1-1) is regarded as the first event in an event group. 
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Fig. 1. The charge distributions for the data 
and Q-ball Monte Carlo events. The soUd black 
line indicates the data and the dotted red 
line indicates the Monte Carlo events with 
a-=10mb. 



Fig. 2. The top figure shows the distribution 
of the number of hits in the largest OD cluster 
for the data and Q-ball Monte Carlo events. 
The bottom figure shows the distribution of the 
goodness of the vertex fit. The lines are same 
as in Fig. 1. 



After the first selection, there are stiU background event groups which are caused 
by cosmic ray muons, electrons from a decay of cosmic ray muons and electrical 
noise caused by large muon pulses. Background event groups still survive the above 
criteria with an approximate rate of 0.06 /sec. In the second reduction, we require 
at least two events which are entirely contained in the ID. 

2-1) The number of hit PMTs in the largest OD hit cluster should be less than 
16, in order to remove primarily cosmic ray muon events as shown in the top of 
Fig. 2. 

2-2) The number of p.e.'s in the ID within SOOnsec time window should be larger 
than 500 and the goodness of vertex fit should be larger than 0.6 (the bottom of 
Fig. 2) to better reject events associated with electrical noise. 

2-3) In an event group two or more events should satisfy both 2-1) and 2-2). 

Two event groups have been found, both containing two successive events. How- 
ever, they are accidental coincidence events consisting of a low energy event and a 
cosmic ray muon within the 1.3 fisec trigger gate. The cosmic ray muon is separately 
recorded over two events since the cosmic ray muon hits at the end of the trigger 
gate. These event groups are not rejected in the former reduction stages because of 
the absence of hit OD PMTs and the large number of total p.e.'s in the ID due to 
the muons. 

We design the following selection criteria to reject the accidental coincidence 
events without losing any Q-ball Monte Carlo events. 

3) If an event group contains two events, for the first event, the ratio of the 
number of hit ID PMTs in a lOOnsec time window at the end of the trigger gate 
to the total number of hit ID PMTs should be less than 50% (see in Fig. 3). 

No group of events has survived after applying all of the selection criteria. The 
number of background event groups due to the accidental coincidence of 2 or more 
atmospheric neutrino events within 100/Ltsec during the 541.7day exposure is esti- 
mated to be 8.6x10"^. 
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Fig. 3. The distribution of the ratio of the num- 
ber of hit ID PMTs in a 100 nsec time window 
at the end of the trigger gate to the total hit ID 
PMTs for the data and the Monte Carlo events. 



Fig. 4. 90% C.L. upper flux limits on the neu- 
tral Q-ball flux as a function of the Q-ball-nu- 
eleon interaction cross section. The circles show 
our results and the dashed line in the region 
below 1 mb cross section is obtained assuming 
that the flux limits should be in inverse propor- 
tion to the square of cross section. 



5. Results 

Since there is no candidate event group, the upper Hmits on the Q-ball flux for 
various mean interaction lengths can be calculated from the following formula, 

Flux < i90%C.L.), (3) 

where Nq = 2.3 for no observed event group, T is the detector live time (541.7days), 
and S^fffl, the effective aperture, which is defined as: 

Seffi^ = cTrr^ X in [m^ sr], (4) 

where e is the detection efficiency estimated by Monte Carlo events and r is the 
impact parameter, which is defined to be 50 m in this analysis. 

In another theory, Q-balls could convert matter into antimatter on their sur- 
faces [8]. In this case, up to 2GeV energy can be released. However, the present 
selection criteria is almost equally efficient for Q-balls from this theory. We also 
find that varying the multiplicity of the generated pions from two to four does not 
change the detection efficiency. The velocity of the neutral Q-ball has uncertainty 
of as much as ±10% [9]. The efficiency is not affected by this uncertainty. Table 1 
summarizes the detection efficiencies and the obtained flux limits. 

Fig. 4 shows the 90% C.L. upper limits on the neutral Q-ball flux as a function 
of the Q-ball cross section. The upper right hatched region is excluded because the 
Q-ball mass density exceeds the Galactic mass density [13]. 

The hatched region is previously excluded by the experiments which searched 
for monopole-catalyzed nucleon decays. The triangle points are the results from 
Kamiokande [14] and the square point is obtained by MACRO [15]. The present 
results, shown in circles, are the most stringent limits on neutral Q-ball flux for 
Q-ball cross section below 200 mb. 
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Table 1 

Summary of flux upper limits on neutral Q-balls for several Q-ball cross-sections. 



Cross section (mb) 


Detection efficiency (%) 


Flux upper limits (cm ^sec ^sr ^] 


0.02 


0.01 


5.0x10-13 


0.05 


0.03 


2.0x10-13 


0.1 


0.07 


7.7x10-11 


0.2 


0.6 


9.1x10-15 


0.5 


2.1 


2.4x10-15 


1 


4.5 


1.1x10-15 


10 


13.6 


3.7x10-1^ 


100 


14.5 


3.4x10-1'^ 


200 


14.4 


3.5xlO-i'5 



In summary, wc report on a search for groups of successive contained events 
at Supcr-Kamiokandc II. We have found no evidence for the successive contained 
event groups in 541.7days of SK-II data. New upper hmits at 90 % C.L. for neutral 
Q-ball flux are obtained. These limits are the most stringent bounds on Q-ball flux 
for Q-ball cross section below 200 mb. 
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Abstract 

Neutrino telescopes are sensitive to a variety of hypothetical super-heavy exotic particles. 
We review the status of current searches for magnetic monopoles, nuclearites, and Q-balls 
using data taken with the AMANDA-II detector. 



1. Introduction 

The existence of magnetic monopoles is mandatory in a large class of Grand Uni- 
fied Theories (GUTs). Monopoles are supposed to be copiously produced as topo- 
logical defects in symmetry breaking phase transitions in the early universe, but 
their density will have been strongly diluted by inflation. GUTs predict monopole 
masses to range from 10® - 10^^ GeV, depending on the symmetry group and unifi- 
cation scale of the underlying theory [1] . The monopole magnetic charge will be an 
integer multiple of the Dirac charge gn ~ e/ {2a), where e is the electric elementary 
charge and a = 1/137 is the fine structure constant. Since magnetic monopoles 
are topologically stable, they should still be present in today's universe and can 
be searched for in cosmic radiation. Once created, monopoles can efficiently be ac- 
celerated in large scale magnetic fields or by gravitation. Monopoles with masses 
below ~ 10^'' GeV can have been accelerated to relativistic velocities [2] and could 
be detected in neutrino telescopes through their direct Cherenkov emissions. In 
some GUTs, monopoles can act as "catalysts" in interactions that violate baryon 
number conservation [3]. In a neutrino telescope, the signature of these "catalyz- 
ing" monopoles would be a series of closely spaced light bursts from nucleon decay 
products produced along the monopole trajectory. Other massive particles have 
also been hypothesized to exist in cosmic radiation: Nuclearites (nuggets of strange 
dark matter) and Q-balls (supersymmetric coherent states of squarks, sleptons and 
Higgs fields, predicted by supersymmetric generalizations of the standard model). 
Electrically neutral Q-balls will have the same experimental signature in a neutrino 
telescope as catalyzing monopoles, although the underlying process [4] is different 
from nucleon decay catalysis. Nuclearites and charged Q-balls might also be de- 
tectable, as, traveling through matter, they would generate a thermal shock wave 
which emits blackbody radiation at visible wavelengths [6,7]. 

2. The AMANDA-II Neutrino Telescope 

AMANDA-II is a neutrino telescope located at a depth between 1500 and 2000 m 
under the ice at the geographic South Pole. A cylindrical volume of roughly 200 m 
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diameter of the Polar ice was instrumented with a total of 677 optical modules 
(OMs), consisting of a photomultiplier tube (PMT) and supporting electronics en- 
closed in a transparent pressure sphere. The OMs were deployed on 19 vertical 
strings, which are arranged in three concentric circles (see figure 1). 

The detector was built in three stages, each in- 
corporating step-wise improvements to the PMT sig- 
nal transmission techniques. The inner ten strings of 
the detector are read out electrically via coaxial or 
twisted-pair cables, while the outermost strings use 
optical fiber transmission. The detector is operated 
with a variety of different triggers, two of which are 
relevant to the analyses presented here: First, the 24- 
fold multiplicity trigger requiring a minimum of 24 
OMs hit within a fixed coincidence window of 2.5 /is, 
and second, a so-called correlation trigger, requiring Fig. 1. Arrangement of the 19 
n OMs to be hit in any group of m adjacent OMs strings of the AMANDA-II de- 
on the same string. For each triggered event, PMT t^^*"'' ™ horizontal plane, 
pulse data is recorded over a time window of ~ 33 ^s. For each OM, up to eight 
subsequent pulses (or hits) can be recorded. The vast majority of triggers are due 
to down-going atmospheric muons, yielding an average event rate of roughly 90 Hz. 




3. Search for Relativistic Magnetic Monopoles 

The number of Cherenkov photons iV^ emitted per path length dx and photon 
wavelength dX radiated from a relativistic magnetic monopole carrying one Dirac 
charge passing through matter with index of refraction n is [8] 

dxdX A2 V e / ^ /32n2 J ' ^ > 

where /3 is the velocity of the monopole as a fraction of the velocity of light. The 
Cherenkov light intensity in ice radiated from the monopole is enhanced by a factor 
{9d ■ n/e)'^ = 8300 compared to the intensity radiated from a particle with electric 
charge e and the same velocity. Thus, in a neutrino telescope, a relativistic magnetic 
monopole will stand out as an extremely bright event relative to the background of 
atmospheric muons. 

We are presently searching for the signal of relativistic magnetic monopoles in 
data taken with AMANDA-II during the year 2000. This data set corresponds to 
194 days of effective livetime. Only events that fulfill the 24-fold multiplicity trigger 
are included. The analysis is "blind" , meaning that the entire data selection chain 
is optimized exclusively on Monte Carlo simulations, and only a small subset of the 
experimental data is used to verify the detector simulation. 

We have simulated the detector response to relativistic magnetic monopoles 
carrying one unit Dirac charge passing the detector's sensitive volume with four 
different velocities (3 — 0.76,0.8,0.9, and (3=1. The background of down-going 
atmospheric muon bundles was simulated with the air-shower simulation package 
CDRSIKA [9]. In several data filtering steps we reject the bulk of low-energy atmo- 
spheric muons. We select events in which a large number of photons are detected. 
Observables like the number of hit OMs, the total number of PMT pulses recorded, 
and the pulse amplitudes are measures for the light deposition in the detector. We 
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Fig. 2. Left: Zenith angle distribution from a likelihood track reconstruction. The region cos(©) < 
corresponds to events in which the particle direction is reconstructed as up-going. Right: Number 
of hit fiber OMs for the "up-going" event sample {cos(0) < 0). This observable serves as final 
cut parameter. In both figures, the full histograms represent the simulated atmospheric muon 
background. Black markers represent 20% of experimental data taken during the year 2000, and 
the four open histograms represent the simulated signal of relativistic magnetic monopoles with 
four different velocities. 



use these observables either as one dimensional cut parameters or as input vari- 
ables to a discriminant analysis [10]. By selecting events with large light deposition 
alone, we can identify relativistic magnetic monopoles among the background of 
down-going atmospheric muons, even if they were arriving from above the horizon. 
In this document, however, we present the search of up-going monopoles only. The 
mass range over which this search is sensitive is limited, since only monopoles with 
masses larger than 10^^ GeV can penetrate the entire Earth and still be relativistic 
upon reaching the detector [11]. 

We reconstruct the particle track using a likelihood method [12] and consider only 
events for which the zenith angle of the track is bigger than 90°, corresponding to 
a particle entering from below the horizon. Some of the down-going atmospheric 
muon bundles are mis-reconstructed as up-going particles, posing a background to 
the search for up-going monopoles (see figure 2). We reject this remaining back- 
ground with a final cut on the number of hit OMs with fiber readout. We use only 
fiber OMs for this last cut. Located at the outermost strings, these OMs define 
the detector's surface area and are consequently most crucial to the acceptance of 
extremely bright signals. Also, their response to large amounts of light is simulated 
most accurately by the detector simulation (which is essential in a "blind" analy- 
sis). We optimize the final cut such that we achieve the optimum "sensitivity". In 
this context, "sensitivity" means the 90% C.L. fiux upper limit that we expect to 
obtain for a given background prediction, if no true signal were present (see [13], 
and references therein). The final cut parameter is shown in figure 2. The signal 
acceptance, and hence the sensitivity of this analysis, depends on the monopole 
velocity. We expect the analysis of 194 days of experimental data to yield sensitiv- 
ities between - 5 x 10"^^ (/? = 1) and - 1 x 10"^^ cm-^g-isr"! {(3 = 0.76). This 
sensitivity range does not yet include the effect of systematic uncertainties. As of 
this writing, our sensitivity is comparable to the best upper limit on the flux of 
relativistic magnetic monopoles published by the BAIKAL Collaboration [14]. 
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Fig. 3. Left: a background event with a triggering and a non-triggering muon. Right: a simulated 
signal event from a particle traveling at speed /3 = 0.01. 



4. Search for Subrelativistic Particles 

In the search for sub-relativistic particles, two hypothesized mechanisms for light 
production are exploited. Magnetic monopoles have been suggested to catalyze 
nucleon decay with a cross section which is large enough to make them detectable 
but uncertain with several orders of magnitude [3]. Here only protons arc considered, 
which have several possible catalyzed decay channels. The channel p e+7r° creates 
an electromagnetic shower with energy close to the proton mass. Other channels 
lose some of their shower energy to neutrinos. 

Neutral Q-balls dissociate baryons in a fundamentally different way, but the 
experimental signature is quite similar [4]. The cross section is their geometric 
size. By limitations given in [5], it ranges from ~ 10~^^cm^ and many orders of 
magnitude upwards. 

The luminosity of thermal shock waves from nuclearites and charged Q-balls. as 
given by [6,7], is also determined by their geometric size. Their luminosity would 
exceed that of neutral Q-balls by several orders of magnitude. 

All particles are simulated with isotropic directions and with speed /? = 10~^. 
In the simulations, the luminosity is expressed as the mean distance A between 
two electromagnetic showers, whose energy is the proton mass. The only value 
for A used so far is 2 cm, and only hydrogen decays in the decay channel c'^tt'' 
with a branching ratio of 0.9 arc considered [15]. The corresponding cross section 
of monopolc catalyzed decay is 9 • 10~^'*cm^, which is at the upper edge of what 
appears to be allowed by the theoretical uncertainties. For neutral Q-balls. oxygen 
nucleon decay is considered, making A correspond to a cross section of 9 • 10"^^cm^. 
For nuclearites and charged Q-balls. the chosen A corresponds to a much lower 
luminosity than that given by [6,7]. 

Simulations show that the correlation trigger is substantially more sensitive to 
this type of signal than the multiplicity trigger. In each triggered event, hits are 
collected during a time window of 33/iS. Background relativistic muons emit light 
during ~ 3/iS, whereas slow particles emit during a large fraction of the time win- 
dow. A comparison is shown in fig. 3. The left picture shows a background event 
with a triggering muon at time 19/xs, and an accidental early non-triggering muon 
at 9fis. The right picture shows a simulated signal event. The signal separation 
from background is based on hits at times when no light from triggering muons is 
expected, the early and late hits outside the interval 16 — 24/^s. 

A period of 113 days in 2001 when a constant correlation trigger definition was 
used, is considered here. The background properties and a preliminary expected 
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sensitivity is determined using 20% of the data. A first filter reduces the data by 
99%, requiring a total of at least 14 early and late hits. The second filtering uses 
additional trigger cleaning and three additional cuts based on early and late hits. 

The events passing the filter have an 
exponential distribution in the number 
of early hits. It is shown in fig. 4, along 
with a fit. Since about 80% of any sig- 
nal events would be expected above 20 
(as seen from simulations), but none 
were found, fig. 4 must be almost signal 
free. Thus, the fit parameters are suit- 
able for background estimation. Using 
the same sensitivity optimization scheme 
as in the relativistic case, the optimal 
final cut for the 80% sample requires 
> 24 early hits. The corresponding expected sensitivity is 5 • 10~^^cm~^s~^sr~^. 
This updates the oral presentation, where the correlation trigger had not been 
included in the simulations and data filtering. 




10 15 20 

Number of early hits 



Fig. 4. An exponential fit parametrizes the back- 
ground distribution in the number of early hits. 



5. Discussion and Outlook 

The AMANDA neutrino telescope is an excellent instrument to search for several 
postulated super heavy exotic particles. In this document, we present first studies 
of the sensitivity of AMANDA to magnetic monopoles, Q-balls and nuclearites. 
The given sensitivities are still preliminary. Specifically, systematic uncertainties 
are not yet included. So far, we have used relatively small sub-sets of the available 
AMANDA data in order to outline our analysis strategies. The sensitivity of the 
two analyses will improve substantially with more data. 
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Abstract 

We present the results of a search for relativistic magnetic monopoles (/9 > 0.8) with the 
Baikal underwater Cherenkov telescope NT200, based on data taken in 1998-2002. The 
upper limit on the flux of monopoles with /3 > 0.8 from the lower hemisphere is the most 
stringent at the present time. 



1. Introduction 

Fast magnetic monopoles with Dirac charge g=68.5e are attractive objects to 
search for with deep underwater neutrino telescopes. The intensity of Cherenkov 
radiation from a relativistic monopole in water is 8300 times higher than that of a 
muon. The hght flux from a monopole is equivalent to the light flux from a muon 
with an energy of about 10 PeV. An optical module (OM) of the Baikal experiment 
could detect such a bright object from distances up to 100 m. 

In 1974 Polyakov [1] and Hooft [2] discovered a monopole solution of the SOS 
Georgi-Glashow model. In a wide class of elementary particle models the mass of 
monopoles may range between 10^ — lO'^* GeV. Propagating through the Universe, 
monopoles could be accelerated by magnetic flelds up to energies of 10^^ — 10^^ GeV 
[3], [4]. Hence monopoles with mass below lO^** GeV are expected to be relativistic. 
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Fig. 1. Distributions of experimental data (solid line), MC atmospheric muons (dashed line) and 
MC from monopolc (dotted line) for cor^z (a), 9 (b), Rrec (c), Ni^a (d). 

The monopole energy loss is about 10 for a Lorentz factor < 10^, but 

significantly rises with energy, by a factor of 1000 for a Lorentz factor 10^ [5]. 
Therefore, light monopoles with mass less than 10^ GeV could not cross the Earth. 
In order to suppress the huge background we search for monopoles only from the 
lower hemisphere and assume a mass range 10^—10^^ GeV. In this paper we present 
the result of a search for rclativistic monopoles with the Baikal neutrino telescope 
NT200, based on data taken in the years 1998 - 2002 (April 1998 - February 2003). 

2. Detector and search strategy 

The neutrino telescope NT200 [6] takes data since April 1998 and consists of 
192 OMs in 8 strings. The OMs arc arranged in pairs, and each pair defines a 
channel. A trigger is formed by the requirement of N^n > 4 fired channels within 
500 ns. For such events, the amplitude and time of all fired channels within a 
time window of 2000ns arc digitized and sent to the shore. The space-time pattern 
of the light recorded from a monopole depends on the water characteristics. The 
maximal absorption length of Baikal deep water is Lobs(480?T.m) = 22 -f- 24 m and 
varies slightly during years. Scattering in Baikal water is characterized by a strongly 
anisotropic function f{6) with a mean cosine cos{9) = 0.85 0.9 and a scattering 
length Ls = 30-f70m. Both scattering length and f{0) vary considerably during the 
year. Our Monte-Carlo simulations indicate that the uncertainty of the scattering 
length from 15 to 30 m leads to an uncertainty of the effective area of nearly 20% 
for fast monopole registration. 

The most natural way to search for such bright objects as fast monopoles is 
to select events with a large number of fired channels, Nhu > 30. To reduce the 
background from atmospheric muons we select only up-going particles. For such 
particles the times of the hit channels increase with rising z-coordinatcs from bot- 
tom to top of the detector. A cut on the value of the z-time correlation, corxz, has 
been applied to suppress downward-moving particles: 
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corrz = — > W 

where ti and Zi are time and z-coordinate of a fired channel, T and z are mean values 
for times and z-coordinates of the event, at and cr^ the root mean square errors for 
time and z-coordinates. The first cut (CUTl) was Nhu > 30 and cottz > 0. These 
requirements leave 1.5 • 10~^ of those events that satisfy the trigger Nhu > 6 on > 
3 strings. CUTl reduces the effective area for monopoles with /3 = 1 by about a 
factor of 2. 

The main background for the fast monopole signature is provided by muon bun- 
dles, high energy muons and cascades from muons. The simulation chain of atmo- 
spheric muons starts with cosmic ray air showers generated with CORSIKA [7]. We 
use the QGSJETl [8] model and a primary composition according to [9]. The MUM 
[10] program is used for muon propagation through the water. To obtain the de- 
tector response to muons passing through the detector volume, we simulate all 
muon energy loss processes, taking into account the characteristics of the OMs and 
the measurement system. In fig. 1 we compare experimental data and simulated 
atmospheric muon events with respect to all parameters which have been used for 
background rejection. The expected distributions for fast (/? = 1) monopoles are 
presented also. One can see that the simulation describes the experimental data 
quite well, even on the level of rare events. 

Within 1003 days live time (between April, 1998 and February, 2003), 21240 
events satisfy CUTl . All events have been reconstructed as single muon tracks. 
For further background suppression we use additional cuts, which essentially reject 
background events from atmospheric muons and at the same time only slightly 
reduce the area for registration of fast monopoles. 

The final cuts were : 

CUTl: N^t > 30 and corrz > 

CUT2: Nhit > 35 and corrz > 0, events should be reconstructable 

CUTS: CUT2 and determined from reconstruction less than 3 

CUT4: CUTS and > 100°; 0: reconstructed zenith angle 

CUT5: CUT4 and Rrec > 10 25 m; Rrec- distance from the center of NT200 

CUT6: CUT5 and corrz > 0.25 ~ 0.65 . 

NT200 took data in various configurations, due to different numbers of oper- 
ating channels for different time periods. The range of cuts on corrz and Rrec 
correspond to different configurations of NT200. No events from the experimental 
sample pass CUT 6. Passing rates versus cut level for experimental and simulated 
events are shown on fig. 2. The figure also shows the dependence of the effective 
area for monopoles as a function of cut level. 



3. Limit on the flux of relativistic monopoles 

The acceptance, A^ff, for monopoles with p = 1,0.9,0.8 has been calculated 
for each configuration of NT200 separately, depending on the number of operat- 
ing channels and concrete values for rejection cuts. From the non-observation of 
candidate events in NT200 and the earlier stages NT36 and NT96 [11], an upper 
limit on the flux of relativistic (/3 > 0.8) monopoles at the 90% confidence level is 
obtained. The cumulative acceptances A^ff ■ T as well as the 90%C.L. upper limits 
are presented in Table 1. These results are an update of an earlier analysis [12], 
mainly differing by a more realistic accounting for light scattering in Baikal water. 
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Fig. 2. Passing rates versus cut level for ex- 
perimental data (triangles), MC for atmo- 
spheric muons (asterix), MC for monopoles 
(squares). 



Fig. 3. Upper limit (90% C.L.) on the 
monopole flux obtained in this work 
(BAIKAL), compared to other experiments. 



In fig. 3 the 90% C.L. upper limit obtained with the Baikal neutrino telescope for 
an isotropic flux of fast monopoles is compared to the limits from the experiments 
Ohya [13] , MACRO [14] and AMANDA [15]. The search for fast magnetic monopole 
with the Baikal underwater detector NT200 results in the strongest experimental 
bound on the flux of monopoles with /? > 0.8. 



4. Outlook 

In order to increase the sensitivity for the searches of bright objects, NT200 was 
significantly upgraded. In April 2005, NT200+ was put into operation. It consists of 
NT200 and 3 external strings spaced at 100 m distance from the center of NT200 
[16]. The main advantage of NT200-I- is the possibility to identify high energy 
cascades which are the main source of background for the monopole search. It 
permits to reject background by softer cuts on Nhu compared to the NT200 analysis. 
We expect an increase of the effective area for fast magnetic monopoles in NT200+ 
of about 1.5 times compared to NT200. 

The construction NT200+ can be considered as a first step in the construction 
of a neutrino telescope with an effective volume of about 1 km'^, (Gigaton Volume 
Detector, GVD). The design of the new telescope GVD has started recently [17]. 
According to preliminary considerations, the telescope GVD will consist of 90-100 
strings. 12-16 OM's will be placed on each string and the instrumented length 
of strings will be 300-350 m. Following our pilot simulations we hope to reach a 
sensitivity to fast monopoles on the level of 5 • 10^^^ cm~^ si~^ after one 
year of data taking. 



5. Conclusion 

The neutrino telescope NT200 in Lake Baikal is taking data since April 1998. 
From the non-observation of candidate events in NT200 and the early stages NT36 
and NT96, an upper limit on the flux of fast monopoles is obtained. With 0.46-10""'^^ 
cm^^- sec^^- sr^^ for (3 = 1 (90% C.L.) it is the most stringent limit at present. 
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Table 1 

Ae// * T and 90% C.L. upper limits on the flux of fast monopoles. 





/3 = 1 


P = 0.9 


(3 = 0.8 


A^ff ■ T NT200 
cm^- sec- sr 


4.84- 10^^ 


3.48- 10" 


1.23-10" 


Ae// ■ T NT36+NT96 
cm^- sec- sr 


0.37 ■ 10" 


0.25 • 10" 


0.094 - 10" 


90 % C.L. upper flux limit 
cm~^- sec"-' - sr"-*^ 


0.46 - 10"" 


0.65 ■ 10"" 


1.83 - 10"" 



In April 2005, the upgraded array NT200+ was put into operation. It will have a 
sensitivity to fast monopoles of about 1.5 times that of NT200. Our future plans 
are connected to the Gigaton Volume Detector. Its sensitivity to fast monopoles 
will be about ~ 5 • 10^^^ cm"^ sr"^ sec~^ for one year of data taking. 
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Abstract 

The ANTARES collaboration is building an undersea neutrino telescope at 2400 m depth 
in the Mediterranean Sea. The experiment aims to detect high-energy cosmic neutrinos 
using a 3D array of 900 photomultipliers (PMTs) arranged in 12 strings. The advantages of 
neutrinos as astrophysical and cosmic messengers are that they open a new window to ob- 
serve known astrophysical objects as well as to look for new Physics, such as dark matter. 
In many supersymmetric models, the favourite dark matter candidate is the lightest neu- 
tralino whose annihilation in the core of massive celestial objects can lead to the emission 
of neutrinos in the subsequent decay chains. The expected performance of ANTARES is 
discussed. 



1. Introduction 

Rapid advances and improvements in observational cosmology are leading to the 
beginning of a precision era in the sense that many of the key cosmological parame- 
ters are being determined with greater and greater accuracy. A series of observations 
have provided sound foundations for our current picture of the universe as mainly 
consisting of non-baryonic, non-luminous and non-absorbing matter, commonly re- 
ferred to as Dark Matter, immersed in some form of energy featuring a negative 
pressure, named Dark Energy. Weakly Interacting Massive Particles or WIMPs can 
explain the growth of structures and are the natural candidate for dark matter 
particles. The original WIMP candidate, and the only one known to exist, is the 
neutrino. However, only a tiny fraction of about 0.3% of the total cosmic density 
can be attributed to neutrinos. Moreover the hot dark matter scenario, though not 
completely ruled out, seems less and less favourable since the simulations of the 
structure formation appear to be inconsistent with observations. Unlike the neu- 
trino, all the other candidates for cold WIMPs are hypothetical and most of them 
are associated with a new symmetry of nature known as supersymmetry (SUSY). 
A hypothetical particle x which is a good candidate for Dark Matter if it satisfies 
a number of constraints coming from current observational results: 

- The cosmological abundance of x's must be about il^ ~ — — 0.23; 

- X niust not be a baryon; 

- X must be electrically and colour neutral, otherwise clectromagnctically and 
strongly interacting x particles would become bound by normal matter forming 
heavy isotopes: very sound upper limits exists for abundances of heavy nuclear 
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Fig. 1. Time differences between hits in a pair 
of OMs in storeys at different distances from 
the Hght source. The top distribution refers to 
OMs at 150 m from the emitter: the effect of 
light scattering is visible. The bottom plot cor- 
responds to a 70 m distant storey: the distri- 
bution is Gaussian and has a width of 0.7 ns. 
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Fig. 2. Reconstructions of two events: height 
(Z) along the string vs. the measured hit times 
wrt an arbitrary zero. Crosses correspond to all 
the hits within a I fis time window, black dots 
to hits in coincidence while those used in the 
event fit are surrounded by a red square. 



isotopes which strongly support, over a wide mass range (1 GeV< <1 TeV), 
the hypothesis of at most a weakly interacting particle. 

2. The ANTARES detector 

The ANTARES collaboration is building a neutrino telescope in the Mediter- 
ranean Sea, at a depth of 2.4 km, located roughly 40 km off the southern coast of 
France at a latitude of 42°50 [1,2]. Neutrinos are detected when they interact with a 
nucleon TV via a charged current {viN — > IX) weak interaction. The detector oper- 
ates by detecting the intensity and arrival time of the Cherenkov radiation emitted 
by charged particles produced in the neutrino interaction. The detector will consist 
of 12 autonomous mooring lines called "strings" . Each of them is equipped with 75 
10-inch Hamamatsu PMTs, arranged in triplets (storeys), oriented at 45 degrees 
with respect to the downward vertical. A series of LED and laser beacons, dis- 
tributed across the storeys and at the bottom of the line, provide a tool for timing 
calibrations (Fig. 1) while a system of acoustic transponders and receivers allows to 
perform precise position measurements (to within a precision of less than 20 cm) . 
The first two lines of the detector have already been deployed and are currently 
being operated, while the complete 12-line detector will be operational by the end 
of 2007. As an example, two reconstructed muon tracks are shown in Fig. 2. 

3. Neutralino detection at Neutrino Telescopes 

The WIMPs in the galactic halo will pass through massive bodies in the galaxy 
and can lose energy by scattering off nuclei so as to become gravitationally bound. 
Over time, WIMPs concentrate near the centres of these bodies and annihilate pro- 
ducing Standard Model particles. The exact rate depends on the time the WIMPs 
have had to accumulate and the annihilation cross section. The products of these 
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annihilations will possibly decay and produce neutrinos which will be able to es- 
cape from the centre of these bodies and would potentially be visible as a neutrino 
flux at the surface of the Earth. The advantage of this detection method, which is 
common to the case of gamma rays too (but different from antimatter searches), 
is that neutrinos do not interact in outer space, and therefore the direction from 
which they arrive points at the location where they were produced. A high-energy 
neutrino signal in the direction of the centre of the Sun or of the Earth is therefore 
an excellent experimental signature, which may stand up against what is the main 
limitation of the technique itself, namely the neutrino background generated by 
cosmic-ray interactions in the Earth's atmosphere. 

Quantitatively, the neutrino flux from neutralino annihilations depends on the 
one hand on the particle physics setup, i.e. on the details of the decay chain of a 
neutralino of a given mass and composition; on the other hand, a crucial role is 
played by the capture versus annihilation balance in the core of the celestial bodies 
and on the physics of the propagation of the relevant SM decay products. The 
differential neutrino flux is given by 

dE^ 47rL>2 ^dE^ ^ ' 

where Ta is the annihilation rate, D is the distance of the detector from the source, / 
is the neutralino pair annihilation final state and are the branching ratios into 

the final state, each giving rise to the energy distribution of neutrinos -jg^ The 
differential flux of neutrino-induced muons reaching the detector is then provided 

by 

- y... ''''' Jo ^""Je^ '""''"^ dE'^ 

where from right to left we can identify 

- : the differential neutrino flux given by eq. 1 ; 

- — ^^^] the production cross section of a /i of energy E'^ from a Vfj^ of energy 

- P (^Efj_, E'^; Ry. the probability for a ^ of initial energy E'^ to end up with an 
energy E^^ after traversing a path-length R; 

- R is the /Li range in the material surrounding the detector. 

At this point, the only missing ingredient for the computation is the annihilation 
rate Ta which, in turn, depends on the details of the neutralino interactions. Some 
assumptions and models must then be made which can affect the final result of the 
computation. Moreover the Sun and the Earth compositions are quite different be- 
cause the Sun is mainly consisting of hydrogen while the Earth is mostly composed 
of nuclei with zero spin. This in turn implies that the flux from the two sources are 
respectively dominated by spin-dependent and spin-independent processes. 



4. ANTARES sensitivity to Neutralinos 

The sensitivity of ANTARES to neutrinos originating from Dark Matter anni- 
hilation in the context of mSUGRA [3] has been extensively studied through Mon- 
teCarlo simulations [4]. One of the results is shown in Fig. 3. The effective volume 
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Fig. 3. ANTARES sensitivity to a muon flux 
from ncutralinos in the Sun. 
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Fig. 4. Effective volume for isotropic upgoing 
neutrinos 



(which is the detector sensitive volume inchiding reconstruction efRciency) is calcu- 
lated as a function of ncutralino energy and optimised with respect to the angular 
size cut (Fig. 4). For any assumed neutralino mass, the expected background inside 
the cone is calculated and the 90% confidence level limit is set using the Feldman 
and Cousins technique [5] to yield a neutrino flux limit. This is in turn converted to 
a muon flux limit using the muon yield per neutrino provided by the DarkSUSY [6] 
program. Superimposed on Fig. 3 are a number of points which correspond to the 
theoretical predictions within the mSUGRA framework. The fall-off in sensitivity 
at low energies arises from using a reconstruction package optimised for higher 
energies: dedicated studies to improve low energy reconstruction efficiencies and to 
evaluate the sensitivity to neutrino fluxes from neutralino annihilation in the Earth 
are under way. 
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Abstract 

The IceCube neutrino observatory, currently under construction at the South Pole, will be 
used to search for particles and interactions beyond the standard model. These searches 
include both direct and indirect searches for dark matter, searches for non-standard neu- 
trino interactions, searches for supersymmetric particles, and searches for both relativistic 
and non-relativistic monopoles. The expected sensitivity of IceCube to exotic physics will 
be discussed. 



1. Introduction 

The IceCube neutrino observatory will be the largest volume instrumented de- 
tector ever constructed. Nearly one cubic kilometer of South Pole ice will be instru- 
mented with highly sensitive photomultiplier tubes and associated readout elec- 
tronics. The large size is driven by the desire to detect high-energy (TeV-PcV) 
extraterrestrial neutrinos. A detector of this size can also open up new possibilities 
for exotic physics searches. 

High energy neutrinos have many properties that make them ideal probes for 
physics beyond the standard model. The small standard model neutrino interac- 
tion cross-section means that any new types of interactions will have pronounced 
effects on the observed neutrino spectra. At this conference, we heard talks on 
several new types of neutrino interactions, including TeV scale gravity [1], micro 
black- holes [2], and supersymmetric stau production [3]. The observation of these 
signatures requires the existence of a high-energy extra-terrestrial neutrino flux. 
Once IceCube establishes the existence and makes quantitative measurements of a 
high energy extra-terrestrial neutrino flux we will use this lepton beam to probe 
physics beyond the reach of accelerators. 

IceCube will also be a unique instrument for dark matter searches. There is a long 
history of indirect dark matter searches in neutrino telescopes utilizing neutrinos 
produced in neutralino annihilations in the center of the Earth or Sun (for a review 
of such searches in AMANDA and IceCube. see the contribution by Hubert [4]). 
There are also dark matter candidates where IceCube has sensitivity for direct 
detection. They are generally very massive and therefore have extremely low fluxes 
due to the dark matter bound, are expected to have virial velocities {(3 w 10"'^), 
and emit photons via non-Cherenkov mechanisms. Three potential candidates are: 
Supersymmetric Q-balls. These are coherent states of squarks and Higgs fields 
that may be produced during the decay of the proposed Affleck-Dine conden- 
sate [5]. Q-balls carry large baryon number but would evade the non-baryonic 
constraint on dark matter from Big Bang Nucleosynthesis . The observation of 
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these objects would validate the Affleck-Dine mechanism as the origin of the 
baryon-antibaryon asymmetry in the universe. These states arc only stable for 
masses > 10^^ GeV. They produce light via the catalysis of nucleon decay (which 
leads to pion production and finally photons). 

Heavy Strangelets. These objects are states with nearly equal numbers of up. 
down, and strange quarks [6]. As they become very massive, they have sizes 
of order atomic dimensions or above. Objects this size passing through ice will 
trigger a thermal shock and produce light via black-body radiation. 

Massive Magnetic Monopoles. The limits on massive monopoles with masses 
> 10^^ GeV can be improved by IceCube. These objects also register in the 
detector via the catalysis of nucleon decay. 

A search for slowly moving massive particles using AMANDA is in progress [7] . 
As with other Cerenkov neutrino telescopes, IceCube will be highly sensitive to 

relativistic monopoles. A fast monopole with unit Dirac charge will produce ^^8300 

times as much Cerenkov radiation as an identical particle with unit electric charge. 

2. Detector 

The possibility of using the South Pole Ice as a neutrino detector was first demon- 
strated by AMANDA [8]. The South Pole has several obvious and non-obvious 
advantages as a location for high energy neutrino detection. The ice thickness at 
the South Pole is ~2800 m, which allows the active instrumentation to be shielded 
from the vast majority of cosmic-ray muons. At depths below ~1400m the trapped 
air bubbles in the ice convert to an air hydrate crystal form with a substantially 
increased scattering length. Below 1500m, the effective scattering length varies be- 
tween '--jl5-45m for wavelengths between 300 and 600nm except near the dust peak 
at ~2100m[9]. For wavelengths less than 400 nm the absorption length averages 
-110m. 

Although the optical properties of water (effective scattering lengths —200 m and 
average absorption length —50m) are superior to ice. South Pole Ice has several 
technical advantages for operation of a neutrino telescope. For the instrumented 
depths (1450m to 2450m), the temperature of the ice varies from -40 C at the highest 
point to -15 C at the lowest point. At these low temperatures, the dark noise rates 
in the photomultiplier tubes are reduced (averaging about 450 Hz for the IceCube 
DOMs). The lack of biolumincscence allows the optical modules to deployed in 
isolation rather than in pairs. While it is difflcult to deploy the optical modules 
in ice, once the IceCube strings are installed they are mechanically stable and do 
not require constant alignment. Several AMANDA strings have been operated for 
nearly 10 years with minor optical module failures. 

IceCube is both larger and technologically superior to AMANDA. The instru- 
mented volume of AMANDA is —0.015 km'^, while IceCube will eventually approach 
1 km'^ once all 70-80 strings are installed (as of this conference, IceCube has 9 de- 
ployed strings). IceCube, however, is more sparsely instrumented with a horizontal 
spacing between strings of 125m compared to 55-75 m for AMANDA. Vertically, the 
modules are spaced every 17m in IceCube versus 10-20m in AMANDA. In addition 
to the in ice optical modules, a surface array (IceTop) based on the same optical 
module technology will be installed. The acceptance of the surface/deep ice joint 
array will approach — 0.3km^sr and will allow cosmic ray composition studies from 
below the knee up to 1 EeV. 
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Fig. 1. Plot of the arrival time of the first photoelectron in the DOM 17m above an LED flasher. 

Both AMANDA and IceCube utilize large (10" in the case of IceCube) photomul- 
tiplier tubes encased in spherical pressure spheres as the active optical elements. 
In the case of IceCube, however, each digital optical module (DOM) acts as an 
autonomous data collection unit. Each DOM has its own high voltage power sup- 
ply and waveform digitization electronics. The digitized waveforms can be buffered 
in the DOMs, and are communicated digitally at a rate of up to 1 MBit/s over 
the cable. In the case of AMANDA, HV was provided from the surface, and the 
photomultiplier signals were sent to the surface analog communications (with cor- 
responding signal degradation). 

For the purposes of exotic physics studies, the in-ice signal digitization will give a 
improvement in the sensitivity of the detector. The waveform digitization is done by 
three analog transient waveform digitizers (ATWDs) that run at three different gain 
settings. This leads to a large increase in dynamic range and will reduce saturation 
for the brightest events (e.g. relativistic monopoles). In addition to the ATWDs 
which record 128 waveform samplings with a sampling rate between 200 and 700 
MHz, each DOM is equipped with a slow 40 MHz digitizer that records the waveform 
for up to 6.4/is. This ability to record multiple photoelectrons over an extended 
period will allow for improved background rejection. 

In the case of slowly moving massive particle searches, the IceCube hardware 
opens up new possibilities as well. The DOMs operate essentially dead-time free. 
Particles moving at 10^'^c will spend a few milliseconds crossing the detector and 
the DOMs are able to record signals during this entire period. The IceCube trigger 
is entirely software based and topological triggers for slow particles are under devel- 
opment. The minimum velocity to which IceCube is sensitive is largely determined 
by the dark noise rates, and these are substantially lower than similar water based 
Cerenkov detectors. 

3. Detector Performance and Sensitivity Estimates 

As of 2006, 9 IceCube strings have been deployed with 540 DOMs. A total of 16 
IceTop stations have been constructed with 64 DOMs. Extensive testing has been 
done on the deployed hardware to verify the performance of the IceCube DOMs and 
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Fig. 2. Simulation of a mass =1 PcV/c^ monopole with 7 = 10 passing through the 2006 nine-string 
IceCube detector. 

data acquisition system[10]. The detector has been used to reconstruct down-going 
muons as well as up-going atmospheric neutrinos. Additionally, LEDs in each DOM 
have been used to test the performance of the array. Figure 1 shows the arrival time 
distribution of the first photon recorded in the DOM 17m above a flashing LED. 
This shows that the timing resolution of the deployed DOMs is better than 2ns. 

First simulations are underway to determine the sensitivity of the IceCube array 
to exotic physics signals. The simulations use a detailed model of South Pole Ice, 
simulations of the photon propagation through the pressure sphere and glass, de- 
tailed modeling of the photomultiplier tube response and electronics, and a sim- 
ulation of the detector trigger. Figure 2 shows a simulated event consisting of a 
M=l PeV/c^, 7 = 10 monopole passing through the 2006 IceCube configuration. 
As monopole events are very bright, the large cross-sectional area of IceCube will 
allow even the present detector stage with 9 strings to perform monopole searches 
with increased sensitivity over existing limits. It is expected that analysis of the 
2006 IceCube data for will have a flux sensitivity of ~10~^^ cm~^s~^sr^^. This will 
represent a factor of 3-4 improvement of the current best limits from Baikal[ll]. The 
full IceCube detector operated for several years will have a sensitivity to relativistic 
monopolcs of ~10^® — 10^^^ cm^^s^^sr^^. 

IceCube will have a similar flux sensitivity to slowly moving massive particles. 
Figure 3 shows the current limits as well as future IceCube sensitivity to massive 
strangelets, neutral Q-balls, and GUT monopolcs that catalyze nucleon decay. In 
this case the flux limit is shown versus the mass of the slowly moving massive 
particle. The diagonal black line is the local dark matter limit [p — 0.3 GeV/cm'^) 
assuming a velocity (3 — The black line horizontal line (IceCube24) is for 

the expected 21-23 string 2007 IceCube configuration, and the blue line is for three 
years of full IceCube. The strongest current limits on monopolcs arc from searches 
for nuclear tracks in ancient mica [12]. IceCube should exceed that sensitivity level 
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Fig. 3. Current flux limits on /3 = 10"'^ massive particles and predicted flux sensitivity for 2007 
24 string IceCube configuration and full IceCube configuration. 

and will not be subject to the same theoretical uncertainties. 
4. Conclusions 

The IceCube neutrino observatory will be the most sensitive detector for many 
of the exotic physics topics discussed at this conference. Extraterrestrial neutrino 
beams will probe particles and interactions beyond the standard model. IceCube 
will conduct both indirect and direct searches for dark matter. The initial experience 
with the IceCube hardware indicates that a rich exotic physics program is feasible. 
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Abstract 

The NESTOR collaboration advances towards a large, deep-sea neutrino telescope in the 
Mediterranean Sea. We highlight the achievements of operating the NESTOR neutrino 
telescope at a depth of about 4000m and outline the next phases of NESTOR for a 
neutrino telescope with large sensitive area for high energy astrophysics. 



1. Introduction 

In attempting to detect intra- and extra-galactic High Energy Neutrinos on Earth 
one is faced with a minuscule signal compared to the overwhelming backgrounds 
due to cosmic rays and atmospheric neutrinos. Furthermore, to observe as much 
as possible of the given neutrino fluxes from extraterrestrial sources, requires vast 
sensitive detection areas, in the order of km^. 

The Neutrino Extended Submarine Telescope with Oceanographic Research, NES- 
TOR, is a deep underwater neutrino detector designed to detect neutrino induced 
muons and showers of charged particles through their Cherenkov radiation produced 
in sea water [1,2]. After a description of the NESTOR site and detector, we high- 
light the main results of operations in deep sea, at about 4000m depth, [3-5] and 
report on the NESTOR plans for future cubic kilometer size Neutrino telescopes. 

2. NESTOR Site and Detector 

Located off the South- West of the Peloponnese, in Greece, near the Ionian Sea 
with deepest waters, 5200m, in the Mediterranean Sea, the NESTOR project, up 
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Fig. 1. A schematic drawing of a NESTOR 
tower made of vertically placed stars with Op- 
tical Modules at their apex, the readout elec- 
tronics in the Titanium sphere and calibration 
modules in-between. 



Fig. 2. Pulse height distribution of one PMT 
using data from 4000m depth (main plot) with 
three spectra, the PMT dark counts, and single 
and double photoelectron, compared to calibra- 
tion data taken in the lab (insert plot). 



to now makes use of a large horizontal plateau 8km x 9km at a mean depth of 
4000?Ti [1] about 7.5 nautical miles from the shore. Several studies of the envi- 
ronmental properties, [6-8], show a water transmission length of 55m ± lOm at a 
wavelength of A = 460nm and minimal underwater currents, few cm/s [9]. 

The NESTOR detector, is composed of rigid hexagonal stars which can be stacked 
vertically to form a tower (see figure 1). The six-arm star (or floor) is made of 
Titanium latticed tubes with a diameter of typically 32m. At each end of each 
arm two Optical Modules are mounted, each with a photomultiplier tube (PMT), 
with one photocathode facing down- wards and the other up- wards, providing Air 
coverage. 

The NESTOR Optical Module [10] is made of a Benthos glass sphere housing a 
15-inch Hamamatsu PMT. Some important characteristics of these PMTs are large 
photocathode area, high quantum-efficiency (20%) and small transit time jitter 
(FWHM w 5.5n.s) at the single photoelectron level. Optical coupling between the 
glass housing and the PMT is improved by glycerine and shielding from the Earth's 
magnetic field is achieved with a mu-metal cage. 

At the centre of a star a Titanium sphere (Im diameter) hosts the electronics 
for voltage control of all PMTs, and the front-end readout electronics developed by 
the Lawrence Berkley National Laboratory [11], for signal acquisition, triggering, 
digitization and data transmission to shore. The electronics can apply a threshold 
to each PMT signal and a selectable majority trigger, allowing event selection by 
requiring a time coincidence for any subset of all PMTs on a star. When a trigger 
condition is satisfied the analogue signals of all PMTs on a star are digitized and 
transmitted to shore, to be monitored in the control room and recorded for offline 
analysis. 

Calibration modules, above and below a star, use bright blue LEDs with con- 
trolled light pulse sequences for in-situ calibration and monitoring of the PMTs. 
The relative timing between the PMTs was measured to be about 0.5ns. 
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Fig. 3. Total trigger rates (solid line), bi- 
olumincscence contribution to total trig- 
ger rates (open circle), experimental trigger 
rates from periods without bioluminescence 
(closed circle). Bioluminescence occurs only for 
1.1% it 0.1% of the active experimental time. 



Fig. 4. Vertical atmospheric muon intensity of 
various experiments. NESTOR measurements 
pointed at by arrow. Insert plot, ditto in dou- 
ble-log scale. 



The junction box mounted on the sea-bottom Anchor Unit connects the telescope 
to the shore through an electro-optical cable. 

3. Results from the deep-sea 

The NESTOR collaboration has deployed and operated environmental sensors 
attached to the Anchor Unit at a depth of 4000m in the Mediterranean Sea (project 
LAERTIS; Laboratory in the Abyss of Europe with Real-Time data Transfer to 
shore for Interdisciplinary Studies) in January 2002 [3] and a telescope module at 
3800m depth in March 2003 [4]. The event samples collected have been used to 
compare the data and Monte Carlo predictions of the cosmic ray muon flux [5]. 
The simulation tools and methods, [12] and the reconstruction, [13] are described 
in detail elsewhere; here we highlight the main results. 

The star deployed at 3800m had a diameter of 12?7i and 12 OMs, 6 up-/down- 
looking OM pairs at the end of each arm. The baseline rate of a PMT in the 
deep sea was measured to be of the order of 50kHz per PMT and constant with 
time and independent of trigger criteria. A typical PMT pulse height distribution 
recorded in deep-sea (fig. 2) is well described by the thermionic PMT noise and 
single and double photoclcctron distributions as measured in calibration runs in 
the laboratory [4]. 

At periods, high PMT rate bursts were observed above the baseline of 50kHz for 
a prolonged duration of 1 — 10s. These are attributed to bioluminuous activity of 
marine organisms. In requiring a majority trigger, the effect of the bio-activity on 
the active experimental time has been estimated, see figure 3, and found to be in 
the order of 1% at the NESTOR site, in agreement with several previous findings 
using autonomous free drop measurements [14]. 

The atmospheric muon flux was found to be in good agreement with MC pre- 
dictions, based on the atmospheric muon model of Okada [15]. Combined with 
previous NESTOR measurements [1], the flux of atmospheric muons arriving at 
the detector depth per unit solid angle dfi, time dt and area dS can be parame- 
terized as dN / {(KldtdS) ~ Iqcos^IO) and the vertical atmospheric muon intensity 
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Fig. 5. The NuBE array showing 4 strings sur- 
rounding the NESTOR tower. Also shown is 
a downward going muon with its Cherenkov 
cone. 




Fig. 6. Concentric arrangement of NESTOR 
towers with variable density detector geometry 
can cover energy ranges from GeV to PeV of 
neutrino induced particles. 



at 3800m depth was found to be, /q = 9.0 x 10"^ ± 0.7 x W-^{stat) ± 0.4 x 
10^^ {syst)cm~^ sr~^ , with a = 4.7 ± 0.5{stat) ± 0.2{syst). Figure 4 shows the 
vertical atmospheric muon intensity as measured by NESTOR along measurements 
of other experiments, [5,13]. 



4. Neutrino Burst Experiment 

Bursts of very-high energy neutrinos (E> arc expected to accompany 

Gamma-Ray Bursts (GRBs) originating in relativistic fireballs. The relativistic fire- 
ball model has been used by Waxman and Bahcall [16] to predict that about 10—100 
neutrino-induced muons per year could be detected in coincidence with GRBs in a 
Ifcm^-area detector. The particles produced by those neutrino interactions have dis- 
tinct signatures; muons producing showers (catastrophic bremsstrahlung) of many 
km long range; electrons or neutral-current events (cascades) producing short range 
but very bright showers. 

NuBE [17], proposes to build a high-energy neutrino detector to detect neutrinos 
of energy E > lOTeV and temporal coincidences (in the order of 100s) between neu- 
trino events and gamma-ray bursts as observed by satellites. The proposed experi- 
ment can be accomplished by deploying a NESTOR tower of four stars surrounded 
by four autonomous strings, as shown in figure 5. A NuBE string will consist of two 
vertically separated photon-detector nodes, with each node made of two clusters 
of Optical Modules. Eight Optical Modules will form a cluster, arranged such as 
to provide omnidirectional coverage. NuBE strings at distances of 300to from the 
NESTOR tower will yield an effective are of > Ikm^ omnidirectional in sensitivity 
for neutrino energies E^, > lOTel^. For GRB durations of about 10 — 100s, an ab- 
solute timing accuracy of a few seconds per year is achievable, with relative timing 
accuracy of 0.3/is in the strings and about 5ns in the tower. While the primary 
identification with GRBs lies in time correlations, the angular resolution capabil- 
ity of FWHM w 14°, see figure 6 and 9 in [5], provides further verification of the 
correlation with GRBs or other astrophysical sources. 
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5. The KM3NeT 

A deep-sea fcm'^-scalc observatory for high energy neutrino astronomy and asso- 
ciated platform for deep-sea science study, cahed KM3NeT [18], is in progress in 
Europe with 37 participating institutions. A fcm'^-scale detector in the Northern 
hemisphere will be complementary to the IceCubc detector [19] at the South Pole. 

Figure 6 depicts a birds eye view of a possible knt^ size detector [20], made of 13 
NESTOR towers of 100m diameter, each made of a stack of 12 stars. 

The central tower would be surrounded by other towers on radii in the order of 
600m from each other. Variable density detector arrays can cover different energy 
ranges from GeV to PeV of neutrino induced particles. The expertise of designing, 
deploying and operating the NESTOR stars and strings as proposed for the NuBE 
aids in the development of the KM3NeT project. 

6. Conclusions 

The NESTOR collaboration has operated neutrino detector modules at a depth 
of about 4000m in the Mediterranean Sea. The good agreement with the published 
measurements of the vertical muon intensity in deep-sea demonstrate the precise 
understanding of the detector response, calibration, simulation and efficiencies. The 
Neutrino Burst Experiment is planned to measure time correlations of v's and 7 
rays of GRBs and to contribute to the development of a knt^ sized neutrino telescope 
in the Mediterranean Sea. 
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Abstract 

Neutrinos are the best candidates to test the extreme Universe and ideas beyond the 
Standard Model of particle Physics. Once produced, neutrinos do not suffer any kind of 
attenuation by intervening radiation fields like the Cosmic Microwave Background and 
are not affected by magnetic fields. In this sense neutrinos are useful messengers from the 
far and young Universe. In the present paper we will discuss a particular class of sources 
of Ultra High Energy Cosmic Rays introduced to explain the possible excess of events 
with energy larger than the Graisen-Zatsepin-Kuzmin cut-off. These sources, collectively 
called top-down, share a common feature: UHE particles are produced in the decay or 
annihilation of superheavy, exotic, particles. As we will review in the present paper, the 
largest fraction of Ultra High Energy particles produced in the top-down scenario are 
neutrinos. The study of these radiation offers us a unique opportunity to test the exotic 
mechanisms of the top-down scenario. 



1. Introduction 

The Ultra High Energy (UHE) neutrino detection is one the most important step 
forward in the Cosmic Ray (CR) Physics. The discovery of neutrinos with energy 
larger than 10^'' eV will start the neutrino astronomy, enabhng the observation of 
the most distant and powerful sources in the Universe. The existence of neutrinos 
with such high energy is intimately related to the observation of Ultra High Energy 
Cosmic Rays (UHECR). 

Soon after the discovery of the Cosmic Microwave Background (CMB) radiation 
it was shown that the flux of UHECR, with an energy larger than 10"'^^ eV, should 
be characterized by a sharp steepening at energies ~ 10^° eV, due to the absorption 
processes on the CMB radiation. This effect is the well known Graisen-Zatsepin- 
Kuzmin (GZK) cut-off [1]. After a few decades of observations the detection of 
the GZK steepening is still one of the major open problems in UHECR physics 
and the experimental data are not conclusive. The 11 Akeno Grand Air Shower 
Array (AG AS A) events with energy larger than 10^" eV [2] contradict the expected 
suppression of the UHECR spectrum. On the other hand the HiRes data seem to be 
consistent with the GZK cut-off picture [2] . If the UHECR primaries arc protons 
and if they propagate rectilincarly, as the claimed correlation with BL-Lacs at 
energy 4 — 8 x 10^^ eV implies [3], than their sources must be seen in the direction 
of the highest energies events with energies up to 2 — 3 x 10^*^ eV detected by HiRes, 
Fly's Eye and AGASA [2]. At these energies the proton attenuation length is only 
about 20 — 30 Mpc and no counterparts in any frequency band was observed in 
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the direction of these UHECR events. This is a strong indication that CR particles 
with energies larger than 10^" eV may have a different origin from those with lower 
energies. 

Models of origin of UHECR fall into two categories, top-down and bottom-up. 
In the bottom- up scenario UHECR originates from cosmic accelerators. In these 
accelerators particles of relatively low energy are brought to UHE through multiple 
interactions at the source. The most promising accelerators known in nature are 
based on the diffusive shock acceleration mechanism, in which the particle acceler- 
ation is realized through multiple interactions with a shock front. This mechanism 
works fairly well in Super Nova Remnants (SNR) that are believed to be the re- 
sponsible for the acceleration of Galactic Cosmic Rays of energy E < 10^^ eV (for a 
review see [4]). At the highest energies different bottom- up scenario have been pro- 
posed, among them, the most promising, are those in which acceleration is realized 
through the interaction with a relativistic shock front in Active Galactic Nuclei 
(AGN) and Gamma Ray Bursts (GRB). In the framework of bottom-up models 
the observed UHECR flux should show the predicted GZK steepening and UHE 
neutrinos are produced by the interaction of UHECR with different backgrounds, 
at the source and during their journey to us. These are the so-called cosmogenic 
neutrinos (first proposed in [5]) that we will not discuss here. 

The presence of an excess of events as claimed by AGASA inspired the intro- 
duction of several exotic models for the production of UHECR. These models, 
collectively called top-down, explain the excess of AGASA and give also a clear 
explanation for the lacking of any counterpart of the highest energy events. Many 
different ideas have been proposed among top-down models: strongly interacting 
neutrinos [6] and new light hadrons [7] as unabsorbed signal carriers, Z-bursts [8], 
Lorentz-invariance violation [9], Topological Defects (TD) (see [10] for a review), 
and Superheavy Dark Matter (SHDM) (see [11] for a review). 

In the present paper we will concentrate our attention on the two last models 
that show common features: UHE particles are produced in the decay of superheavy 
particles, that we shall call collectively X particles, with a typical mass of the order 
of the Grand Unified energy scale Mqut — 10^"* eV. In the case of TD the X particle 
once produced, by the internal dynamics of the defect or through the interaction 
of different defects, immediately decays. While in the case of SHDM the X particle 
itself is long-lived contributing to the Dark Matter of the universe. 

From the point of view of elementary particle physics the X particle decay process 
proceed in a way similar to e^e~ annihilation into hadrons: two or more off-mass- 
shell quarks and gluons are produced and they initiate a QCD cascade. Finally the 
partons are hadronized at the confinement radius. Most of the hadrons in the final 
state are pions and thus the typical prediction of all these models is the dominance 
of neutrinos and photons at the highest energies E > 5 x 10^^ eV. It is important 
to stress here that these models predict neutrino fluxes most likely within reach of 
the first generation neutrino telescopes such as AMANDA, and certainly detectable 
by future kilometer-scale neutrino observatories [12]. 

2. Hadrons spectrum in X decay 

The first step to determine the neutrino flux produced in the decay of X par- 
ticles is the determination of the hadron spectrum. Moreover, this evaluation is 
particularly important because it represents a direct signature of the production 
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mechanism that, in principle, can be detected experimentally. As discussed in the 
introduction, the mass of the decaying particle, AIx, that represents the total CMS 
energy ^/s, is in the range 10^^ - 10^^ GeV. 

The existing QCD Monte Carlo (MC) codes become numerically unstable at 
much smaller energies, e.g., at ^/s ~ 10'' GeV and the computing time increases 
rapidly going to larger energies. In this section we will briefly review the main results 
obtained, in the computation of the top-down spectrum of UHE particles, using two 
different computational techniques: one based on a new MC scheme [13,14] and the 
other based on the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equations [13,16]. In both cases SUSY is included in the computation. 

Monte Carlo simulations are the most physical approach for high energy calcula- 
tions which allow to incorporate many important physical features as the presence 
of SUSY partons in the cascade and coherent branching. The perturbative part of 
a QCD Monte Carlo simulation is quite standard with the inclusion of SUSY. For 
the non-perturbative hadronization part an original phcnomenological approach is 
used in Ref. [13]. The fragmentation of a parton i into an hadron h is expressed 
through perturbative fragmentation hmction of partons Dj {x, Mx), that represents 
the probability of fragmentation of a parton i into a parton j with momentum frac- 
tion X = 2p/Mx, convoluted with the hadronization functions ,fj'{x,Qo) at scale 
Qq, that is imderstood as the fragmentation function of the parton i into the hadron 
h at the hadronization scale Qq ~ 1.4 GeV [13]. To obtain the fragmentation func- 
tions of hadrons one has: 

Dlix,Mx)^Y. t-Di{-,Mx)fl'iz,Qo) (1) 

where the hadronization functions do not depend on the scale Mx ■ This important 
property of hadronization functions allows the determination of //'(x, Qo) from the 
available LEP data, D'^{x, Mx) at the scale Mx = Mz- 

The fragmentation functions D^{x^Mx) at a high scale Mx can be calculated 
also evolving them from a low scale, e.g. Mx = Mz, where they are known experi- 
mentally or with great accuracy using the MC scheme. This evolution is described 
by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [15] which 
can be written as 

dtD^ = E ^P'^^ (^) ® D'jix/z, t) , (2) 

3 

wheret = ln(s/so) is the scale, Cg) denotes the convolution /(g)^ = dx/xf{x)g{x/ z), 
and Pij is the splitting function which describes the emission of parton j by parton 
i. Apart from the experimentally rather well determined quark fragmentation func- 
tion D'^{x, Mz), also the gluon fragmentation function Dg{x, Mz) is needed for 
the evolution of Eq. (2). The gluon FF can be taken either from MC simulations 
or from fits to experimental data, in particular to the longitudinal polarized e+e" 
annihilation cross-section and three-jet events. The first application of the DGLAP 
method for the calculation of hadron spectra from decaying superheavy particles 
has been made in Refs. [16]. The most detailed calculations have been performed 
by Barbot and Drees [16], where more than 30 different particles were allowed to 
be cascading and the mass spectrum of the SUSY particles was taken into account. 
The results obtained with the two different techniques discussed above agree fairly 
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well [13]. The accuracy in the hadron spectrum calculations has reached such a 
level that one can consider the spectral shape as a signature of the model. The 
predicted hadron spectrum is approximately cx dE/E^'^ in the region of x relevant 
for UHECR observations. 

3. Spectra of Neutrinos, Photons and Nucleons 

The spectra of neutrinos and photons produced by the decay of superheavy parti- 
cles are of practical interest in high energy astrophysics and can be computed from 
the decay of charged pions [13]. The FFs for charged pions and protons+antiprotons 
can be determined, following [13], from the FFs of hadrons Dh simply introducing 
the ratios SNix) and e^{x) as: Dn{x) ~ £iq(x)Dh(x) and Dt^{x) = eTr(x)Dh(x). 
The spectra of pions and nucleons at large Mx have approximately the same shape 
as the hadron spectra, and one can use in this case = 0.73 ± 0.03 and En = 
0.12 ±0.02 [13]. 

An interesting feature of the up-dated calculations performed in [13] and by 
Barbot and Drees in [16] is the ratio of photons to nucleons. 7/iV. At x ~ 1 x 10""^ 
this ratio is characterized by a value of 2 - 3 only [13]. This result has an important 
impact for SHDM and topological defect models because the fraction of nucleons 
in the primary radiation increases. However, in both models photons dominate (i.e. 
their fraction becomes > 50%) at > (7 - 8) x 10^^ eV. 

Let us now concentrate our attention on UHECR from superheavy dark matter 
(SHDM) [18] and topological defects (TD) [19]. The comparison of the UHECR 
spectrum obtained with the AGASA data, will provide us with the correct neutrino 
flux normalization. 

Production of SHDM particles naturally occurs in a time-varying gravitational 
field of the expanding universe at the post-inflationary stage. The relic density of 
these particles is mainly determined (at fixed reheating temperature and inflaton 
mass) by their mass Mx- The range of practical interest is (3 — 10) x lO^'^ GeV, 
at larger masses the SHDM is a subdominant component of the DM. SHDM is 

,halo Jialo 

accumulated in the Galactic halo with the overdensity 5 — ^-^^r = „ — 1 where 

« 0.3 GeV/cm3, p„ = 1.88 x IQ-'^^h^ g/cm^ and ncDuh^ = 0.135 [20]. With 
these numbers, i5 « 2.1 x 10^. Because of this large local overdensity, UHECRs from 
SHDM have no GZK cutoff. 

dumpiness of SHDM in the halo can provide the observed small-angle clustering. 
The ratio rx = ^x(to/Tx) of relic abundance Qx and lifetime tx of the X particle 
is fixed by the observed UHECR flux as rx ^ 10^^^. In the most interesting case 
of gravitational production of X particles, their present abundance is determined 
by their mass Mx and the reheating temperature Tfj. Choosing a specific particle 
physics model one can fix also the life-time of the X particle. There exist many 
models in which SH particles can be quasi-stable with lifetime tx ^ 10^*^ yr. 
The measurement of the UHECR fiux, and thereby of rx, selects from the three- 
dimensional parameter space {Mx-, Ifl, r^) a two-dimensional subspace compatible 
with the SHDM hypothesis. 

In Figure 1 (left panel) we have performed a fit to the AGASA data using the 
photon flux from the SHDM model and the proton flux from uniformly distributed 
astrophysical sources. For the latter we have used the non-evolutionary model of 
[21]. The photon flux is normalized to provide the best fit to the AGASA data at 
E>Ax 10^9 eV. 
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Fig. 1. [Right Panel] Comparison of SHDM prediction with the AGASA data. The calculated 
spectrum of SHDM photons is shown by the label 7 and by the label pextr the spectrum of 
extragalactic protons (uniformly distributed astrophysical sources). The sum of these two spectra 
is shown by the thick black curve. The red thick line is the SHDM neutrino flux. [Left Panel] 
Diffuse spectra from necklaces. The red thick curve shows neutrino flux, the black thick curve is 
the sum of protons and photons fluxes produced by necklaces (labeled thin black lines). 



One can see from the fit of Figure 1 (left panel), that the SHDM model can explain 
only the excess of AGASA events at > 1 x 10^° eV: depending on the SHDM 
spectrum normalization and the details of the calculations for the extragalactic 
protons, the flux from SHDM decays becomes dominant only above (6—8) x 10^^ eV. 

Topological Defects (for a review see [22] and reference therein) can naturally 
produce UHE particles. The following TD have been discussed as potential sources 
of UHE particles: superconducting strings, ordinary strings, monopolonium (bound 
monopolc-antimonopolc pair), monopolonia (monopolc-antimonopole pairs connec- 
ted by a string), networks of monopoles connected by strings, vortons and necklaces 
(see Ref. [22] for a review and references). Monopolonia and vortons are clustering 
in the Galactic halo and their observational signatures for UHECR are identical to 
SHDM. However the friction of monopolonia in cosmic plasma results in monopolo- 
nium lifetime much shorter than the age of the universe. Of all other TD which are 
not clustering in the Galactic halo, the most favorable for UHECR arc necklaces. 

Necklaces are hybrid TD produced in the symmetry breaking pattern G ^ H x 
U{1) — > H X Z2. At the first symmetry breaking monopoles are produced, at the 
second one each (anti-) monopole get attached to two strings. This system resembles 
ordinary cosmic strings with monopoles playing the role of beads. Necklaces exist 
as the long strings and loops. The symmetry breaking scales of the two phase 
transitions, rjm and rjs, are the main parameters of the necklaces. They determine 
the monopole mass, m Anrim/e, and the mass of the string per unit length fi ~ 
27r77^. The evolution of necklaces is governed by the ratio r ~ m/iid, where d is the 
average separation of a monopole and antimonopole along the string. As it is argued 
in Ref. [23], necklaces evolve towards configuration with r ^ 1. Monopoles and 
antimonopoles trapped in the necklaces inevitably annihilate in the end, producing 

heavy Higgs and gauge bosons {X particles) and then hadrons. The rate of X 

2 

particles production in the universe can be estimated as [23] fix ~ tHdx ' ^■'^^'^^ ^ 
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is the cosmological time. 

The photons and electrons from pion decays initiate e-m cascades and the cascade 
energy density can be calculated as Wcas = 5/^'"^^ Jo° W (i+z)^ = fA^^F, where 
z is the redshift and /tt ^ 1 is the fraction of the total energy release transferred to 
the cascade. The parameters of the necklace model for UHECR are restricted by 
the EGRET observations [24] of the diffuse gamma-ray flux. This flux is produced 
by UHE electrons and photons from necklaces due to e-m cascades initiated in colli- 
sions with CMB photons. In the range of the EGRET observations, 10^ - 10^ MeV, 
the predicted spectrum is (x with a = 2 [25]. The EGRET observations de- 
termined the spectral index as a ~ 2.10 ± 0.03 and the energy density of radiation 
as Wobs ~ 4 X 10~^ eV/cm'^. The cascade limit consists in the bound Wcas < ^^obs- 
According to the recent calculations, the Galactic contribution of gamma rays to 
the EGRET observations is larger than estimated earlier, and the extragalactic 
gamma-ray spectrum is not described by a power-law with a ~ 2.1. In this case, 
the limit on the cascade radiation with a = 2 is more restrictive and is given by 
oJcas < 2 X 10~^eV/cm'^; we shall use this limit in further estimates. Using Wcas with 
fn ~ ^ and to = 13.7 Gyr [20] we obtain from the limit on the cascade radiation 
r^H < 8.9 X 10^7 GeV^. 

The important and unique feature of necklaces is their small separation D, which 
ensures an high density. The distance D is given hy D ^ r^^^^to [23]; since r^/i 
is limited by e-m cascade radiation we can obtain a lower limit on the separation 

between necklaces as D - ( j/j^^ )^^^^o > 10(^/10^ GeV^)!/^ kpc, this small 
distance is a unique property of necklaces allowing the unabsorbed arrival of parti- 
cles with the highest energies. The fluxes of UHECR from necklaces are shown in 
Figure 1 (right panel). We used in the calculations r^/z = 4.7 x 10^^ GeV^ which 
corresponds to Wcas = 1-1 x 10~^ eV/cm"^, i.e. twice less than allowed by the bound 
on Wcas- The mass of the X particles produced by monopole-antimonopole annihi- 
lations is taken as Mx = 1 x lO^"' GeV. From Figure 1 (right panel) one can see 
that the necklaces model for UHECR can explain only the highest energy part of 
the spectrum, with the AGASA excess somewhat above the prediction. Thus UHE 
particles from necklaces can serve only as an additional component in the observed 
UHECR flux. This result has a particular impact on the possible UHE neutrino 
detection. In fact, the necklaces model is only under constrained by the available 
UHECR data, in this context only a clear UHE neutrino observation with a typical 
spectrum as in figure 1 (right panel) can confirm (or falsify) the model. 
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Abstract 

When high-energy cosmic rays impinge on a dense dielectric medium, radio waves are 
produced through the Askaryan effect. At wavelengths comparable to the length of the 
shower produced by an Ultra-Ifigh Energy cosmic ray or neutrino, radio signals are an 
extremely efficient way to detect these particles. This approach offers, for the first time, the 
realistic possibility of measuring UHE neutrino fluxes below the Waxman-Bahcall limit. 



1. Introduction 

Askaryan predicted as early as 1962 [1] that particle showers in dense media 
produce coherent pulses of microwave Cerenkov radiation. Recently this prediction 
was confirmed in experiments at accelerators [2] and extensive calculations have 
been performed on the development of showers in dense media to yield quantita- 
tive predictions for this effect [3]. The pulses emitted when UHE particles strike 
the lunar regolith are detectable at Earth with radio telescopes, an idea first pro- 
posed by Dagkesamanskii and Zheleznyk [5]. Several experiments have since been 
performed [6,7] to find evidence for UHE neutrinos. All of these experiments have 
looked for this coherent radiation near the frequency where the intensity of the 
emitted radio waves is expected to reach its maximum. Since the typical lateral 
size of a shower is of the order of 10 cm the peak frequency is of the order of 3 GHz. 

We have proposed [8] a different strategy to look for the radio waves at con- 
siderably lower frequencies where the wavelength of the radiation is comparable 
in magnitude to typical longitudinal size of showers. The lower intensity of the 
emitted radiation, which implies a loss in detection efficiency, is compensated by 
the increase in detection efficiency due to the near isotropic emission of coherent 
radiation. The net effect is an increased sensitivity by several orders of magnitude, 
for the detection of UHE cosmic rays and neutrinos [8] . 

2. Model for Radio Emission 

There exist two rather different mechanisms for radio emission from showers trig- 
gered by UHE cosmic rays or neutrinos. One is the emission of radio waves from a 
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shower in the terrestial atmosphere. Here the primary mechanism is the synchrotron 
acceleration of the electrons and positrons in the shower due to the geomagnetic 
field, called geosynchrotron radiation, which has recently been confirmed with new 
digital radio techniques [9]. The second mechanism, known as the Askaryan ef- 
fect [1], applies to showers in dense media, such as ice, salt, and lunar regolith, 
where the front end of the shower has a surplus of electrons. Since this cloud of 
negative charge is moving with a velocity which exceeds the velocity of light in the 
medium, Cerenkov radiation is emitted. For a wavelength in the radio-frequency 
range, coherence builds up and the intensity of the emitted radiation reaches a 
maximum. 

The intensity of radio emission from a hadronic shower, with energy Es, in the 
lunar regolith, in a bandwidth Av at a frequency v and an angle 9, can be param- 
eterized as [8] 

3.86 X 104 e-^' At/ / sing E^ rfmoon Y 

t{U,v,t.s} 100 MHz Uin0elO2O eVdz/o(l + Wi^o)^-^^)/ 

where Z = (cos 6* - l/n){-j^=^ (^)' " ^'^ ^ distance to the 

observer, and dmoon ~ 3.844 x 10® m is the average Earth-Moon distance. The angle 
at which the intensity of the radiation reaches a maximum, the Cerenkov angle, is 
related to the index of refraction (n) of the medium, cos 0c = ^/n. 

The spreading of the radiated intensity around the Cerenkov angle, Ac, is, on 
the basis of general physical arguments, inversely proportional to the shower length 
and the frequency of the emitted radiation. Based on the results given in Ref. [10] 
it can be parameterized as 



1 



where L{Es) is the shower length which depends on the energy. In Fig. 1 the spread- 
ing angle is compared at two frequencies with different phenomenological descrip- 
tions and with different analytical models. At the higher frequency all descriptions 
give the same result while at the lower frequency the prediction of 1 is consistent 
with the analytical models. 

In our simulations we have taken into account the attenuation of radio waves in 
the Lunar regolith. As a mean value for the attenuation length for the radiated 
power we have taken \r = (9/i/[GIIz]) m [11]. 

In Fig. 2 we compare the detection limits for UHE neutrinos at different frequen- 
cies with the results obtained from the GLUE experiment [7]. One sees that with 
decreasing frequency one loses sensitivity at lower-energies. This follows directly 
from 1 since with decreasing frequency the maximum signal strength decreases. If 
the energy is more than a factor 4 above this threshold value the detection limit 
improves rapidly with decreasing radio-frequency until one reaches a frequency of 
100 MHz where one obtains the optimum sensitivity. Decreasing the frequency even 
lower provides no gain since the detection limit has already reached the optimum, 
given by the heavy black line. 

Our result at 2.2 GHz lie close to that of the GLUE experiment. The dashed- 
dotted curve in Fig. 2 shows the results of a calculation where wc have reproduced 
the simulation for the GLUE experiment, i.e. including (in a somewhat simplified 
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Fig. 1. The angular spread around the 
Cerenkov angle. The left (right) hand displays 
the results for 2.2 GHz (lOOMHz) respectively. 
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Fig. 2. Flux limits (assuming a null observa- 
tion) for UHE neutrinos as can be determined 
in a 100 hour observation. 



manner) the effects of averaging over lunar-surface slopes of 10°. This result lies 
close to the published limits (the open squares in Fig. 2). 



3. WSRT &: LOFAR predictions 

The Westcrbork Synthesis Radio Telescope (WSRT) consists of fourteen 25 m 
parabolic dishes located on an east-west baseline extending over 2.7 km. The low 
frequency band which concerns us here covers 115-170 MHz [12]. Each WSRT 
element has two receivers with orthogonal dipoles enabling measurement of all 
four Stokes parameters. In tied-array mode the system noise at low frequencies is 
Fnoise =600 Jy. To observe radio bursts of short duration, the new pulsar back- 
end (PuMa II) is used. It can provide dual-polarization baseband sampling of eight 
20 MHz bands, enabling a maximum time resolution of approximately the inverse 
of the bandwidth. In the configuration which we propose to use, four frequency 
bands will observe the same part of the moon with the remaining four a different 
section. In total, coverage of about 50% of the lunar disk can be achieved. 

An even more powerful telescope to study radio flashes from the moon will be the 
LOFAR array [13]. With a collecting area of about 0.05 km^ in the core (which can 
cover the full moon with an array of beams), LOFAR will have a sensitivity about 
25 times better than that of the WSRT. LOFAR will operate in the frequency band 
from 115-240 MHz where it will have a sensitivity of about Fnoise =20 Jy. 

The probability for finding during the complete observation time a signal simul- 
taneously in four frequency bands with a power of 20x Fnoise per frequency band, 
where Fnoise is the noise level per polarization direction, as a random fluctuation of 
the background noise is less than 0.001 (equal to Scr significance). For this reason 
we have assumed in the calculations a detection threshold of 25x Fnoise for both 
the WRST and the LOFAR telescopes. 

A simulation for LOFAR, taking = 120 MHz, bandwidth of = 20 MHz, 
and an observation time of 30 days is shown in Fig. 3 for neutrinos. The results 
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be determined with WSRT and LOFAR obser- 
vations (see text) are compared with various 
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Fig. 4. Flux limits on UHE cosmic rays as can 
be determined in a 1 year observation with the 
LOFAR antenna system compared to that with 
SKA in 2 frequency bands. 



arc compared with the limits that can be obtained from a presently proposed ob- 
servation for 500 hours at the WSRT observatory assuming a 50% Moon coverage. 
The predicted spectrum of GZK neutrinos is taken from Ref. [14]. The prediction of 
top-down (TD) neutrinos is based on the decay of a topological defect with a mass 
of 10^'* eV. The Waxman-Bahcall (WB) limit [15], based on theoretical arguments, 
is an upper limit on the neutrino flux which is consistent with the data on the fluxes 
of UHE cosmic rays. The limits are also compared with those from the RICE [16], 
GLUE [7], ANITA [17], and FORTE [4] experiments which are calculated as model 
independent limits similar to our limits. The limit from the RICE [16] experiment 
has been calculated assuming different power-law spectra for the neutrinos. In gen- 
eral such a limit lies below the model-independent limit [4] . 



4. Outlook 

Observations at the WSRT are in progress. 200 hours observation time has been 
assigned to the experiment with a promise for an additional 300 hours. In the first 
observations we have shown that the subtraction of RFI is well under control. So-far 
we have not seen evidence for a signal from the moon and we are in the process of 
setting improved limits. 

For the longer term future observations are planned with the Lunar Orbiting 
Radio Detector (LORD) and with the Square Kilometer Array (SKA) telescope. 
The LORD project is a collaboration between Russia (LPI, Lavochkin Association, 
MSU, JINR) and Sweden (ISP) [18] which plans to launch a satelhte (around 2009) 
in an orbit around the Moon. For the detection of the radio pulse the satellite will 
carry an antenna system which is sensitive to a wide band, 100-1000 MHz. The 
SKA project [19] is a collaboration between institutions in 18 countries. The SKA 
is an aperture synthesis radio telescope where the location of the site will be decided 
in 2007. The system will operate in a very wide frequency band, from 100 MHz to 
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25 GHz. The detection sensitivity of SKA for UHE particles is even greater than 
that of LORD. To show the possibihties with these new antennas, the predicted 
sensitivities are given for observations in the low, 100-300 MHZ, (LFB) and the 
medium, 300-500 MHz (MFB) frequency bands with SKA. The exciting possibility 
of a real measurement of the GZK-neutrino flux with several hundreds of count will 
become available. 
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Probing the variation of relic neutrino masses with 
extremely high-energy cosmic neutrinos 



Deutsches Elektronen- Synchrotron DESY, Notkestrajie 85, 22607 Hamburg, Germany 



We analyze the prospects for testing the cosmic neutrino background and its interpre- 
tation as source of Neutrino Dark Energy with the neutrino telescope LOFAR. 



1. Introduction 

Certainly, one of the most challenging questions in modern Cosmology and Par- 
ticle physics is, what is the nature of Dark Energy? Strong observational evi- 
dence hints at the existence of this smooth, exotic energy component which drives 
the apparent accelerated expansion of the universe. Recently, Fardon, Nelson and 
Weiner [1] have shown that Big Bang relic neutrinos - the analog of the cosmic 
microwave photons - are promoted to a natural Dark Energy candidate if they 
interact through a new non-Standard Model force. Due to this scalar force, the ho- 
mogeneously distributed relic neutrinos can form a negative pressure fluid and thus 
exhibit just the right properties to act as Dark Energy. As a further consequence 
of this new interaction, the neutrino mass becomes a function of neutrino energy 
density which decreases as the universe expands. Thus, intriguingly, the neutrino 
mass is not a constant but is promoted to a dynamical quantity. After discussing the 
details of this so-called Mass Varying Neutrino (MaVaN) scenario, I will consider an 
astrophysical possibility of testing it with extremely high energy-cosmic neutrinos 
(EHECi/). A more detailed discussion of the results may be found in Rcf. [2]. 

2. Mass Varying Neutrinos 

I will concentrate on a concrete realization of the non-Standard Model neutrino 
interaction as preferred in the literature [1,3] which implements the seesaw mecha- 
nism [4] . Generically, varying mass particle scenarios exhibit a so-called dark sector 
which only consists of Standard Model singlets. In the considered case of MaVaNs it 
contains a light scalar field, the acceleron A, which has an associated fundamental 
potential Vo{A). The acceleron interacts with a second field of the dark sector, a 
right-handed neutrino TV, through a Yukawa coupling, kANN, and thus generates 
its mass Mn{A) = kA. The dark force mediated by the acceleron is transmitted 
to the active neutrino sector via the seesaw mechanism which in addition provides 
a natural explanation for the smallness of the left-handed neutrino mass TOi/. At 
scales well below 100 GcV the Lagrangian contains a Majorana mass term [1,2], 
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for the active left-handed neutrino, where the Dirac type mass m^) originates from 
electroweak symmetry breaking and Mjv 3> has been assumed. Consequently, 
the neutrino mass mi,{A) is light and, since it is generated by the value of the 
acceleron, neutrinos interact through a new force mediated by A. 

The coupling leads to a complex interplay between the acceleron and the neutri- 
nos which links their dynamics. Since the neutrino energy density is a function of the 
neutrino mass m^(y^), it becomes an indirect function of the value of the acceleron, 
Piyi'rn^{A), z). As a direct consequence, it stabilizes the acceleron by contributing 
to its effective potential, 



Vcs{A, z) = p^{m^{A), z) + Va{A), where 



Pu{mi,{A),z) : 



ev + 1 



(2) 
(3) 



with z denoting the redshift and T^^iz) = riy(,(l + z) the neutrino temperature with 
TjyQ ^ 1.69 X 10^** eV. Since cosmic expansion causes a dilution of the neutrino 
energy density, also Y^s evolves with time. Assuming the curvature scale of Vcft and 
thus the mass of A to be much larger than the Hubble scale, dV^^jdA? = m}/^ ^ 
77^, the adiabatic solution to the equations of motions apply [1]. As a result, the 
acceleron instantaneously tracks the minimum of its effective potential, the total 
energy density of the coupled system, and thus varies on cosmological time scales. 
Consequently, the neutrino mass my(^A) which is generated from the value of the 
acceleron, is not a constant but promoted to a dynamical quantity. As long as 
drriu/ dA does not vanish, its time variation is determined by 



dVcji{A) dmy{A) ( dpi,{m^,z) 



OA 



dA 



dVo{m^) 



drrii, 



0. (4) 



Fig. la shows the evolution of both the effective potential V;,ft(.4, z) as well as its 
minimum due to changes in Pu{z). 

In the following I will address a possible shortcoming of MaVaN models and how 
it can be avoided. Since MaVaNs attract each other through the force mediated by 
the acceleron, they are possibly subject to a phenomenon similar to gravitational 
instabilities of cold dark matter (CDM) [5]. As long as the neutrinos are still rel- 
ativistic, their random motions can prevent a collapse. However, when they turn 
non-relativistic the system can become unstable leading to the possible formation 
of so-called 'neutrino nuggets' [5]. If the neutrinos really clump, they would redshift 
like CDM with w ^ ^ —\ and thus cease to act as Dark Energy. However, firstly, 
certain constraints on the scalar potential Vb and on the function m^{A), can lead to 
stable MaVaN models even in the highly non-relativistic regime [6] . Secondly, neu- 
trino oscillation experiments allow one neutrino to be still relativistic today, which 
could be responsible for cosmic acceleration until the present time [5,3]. In Ref. [3] 
the latter case emerges naturally after a straightforward super-symmetrization of 
the standard MaVaN scenario. However, the modified model relies on a slightly 
different acceleration mechanism known from Hybrid inflation [7] to be considered 
in the following. A scalar field keeps a second scalar field in a metastable minimum 
due to its large value and the energy density stored in the false minimum can drive 
cosmic acceleration. In the MaVaN hybrid model the first scalar field can be identi- 
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Fig. 1. a) Evolution of the effective potential V^f({rni,,z) and the neutrino {z) due 

to changes in the neutrino energy density pi/{z)\ b) Exact m^^ (^) (solid), approximated by 
•mu^(z) cx: (1 + z)"^ and m,j^{z) oc (1 + 2)-^/^ (dashed and dotted) in the low and high red- 
shift regime, respectively, assuming m^Q j = 10~^ eV and a normal neutrino mass hierarchy such 
that neutrino oscillation experiments fix m^g ^ = 8-3 x 10~^ eV and mug 3 = 5.17 X 10^^ cV today. 

fied with the acceleron, which is driven to larger values due to the stabilizing effect 
of the fermionic neutrino background. It keeps the scalar partner N of the light- 
est neutrino, naturally present in a super-symmetric theory, in a false minimum 
until it reaches a critical value y^crit ■ The combined scalar potential Vq {A, N) ap- 
pears as dark energy with w ^ —1 and can therefore drive accelerated expansion. 
This hybrid model provides a microscopic origin for a quadratic scalar potential 
Vo(^) oc and thus according to eq. 4 fixes the neutrino mass evolution. After 
generalizing eq. 4 to include three neutrino species, the evolution of the neutrino 
masses in the low rcdshift regime is found to be well approximated by a simple 
power law [2], 

m^.{z) ~ m^.,o{^ + , where m^.fi — 771^.(0) and i = 1,2,3, (5) 

as shown in fig. lb, where the lightest neutrino is assumed to be still rclativistic 
today, TTii^j = 10~^ eV{T^g. Finally, we are in a position to do MaVaN phenomenol- 
ogy. In the following I will discuss a possible astrophysical test [2] for the neutrino 
mass variation involving extremely high-energy cosmic neutrinos. 

3. Extremely high-energy cosmic neutrinos 

We arc living in exciting times for extremely high-energy cosmic neutrinos [8,2]. 
Existing and planned observatories cover an energy range of 10^ GcV < i?o < 
10^^ GeV and promise appreciable event samples (cf. [8,2] and references therein). 
Thus it seems timely to consider the diagnostic potential of EHECi^'s for astro- 
physical processes. A particular example will be discussed in the following. If 
the energy of an EHEC;^ coincides with the resonance energy, Ej"^^ = = 

4.2 X 10^2 GcV, i 1,2,3, of the process vv Z, an EHECi^ can anni- 

hilate with a relic anti-neutrino and vice versa into a Z Boson [9] with mass Mz- 
This exceptional loss of transparency of the cosmic neutrino background (Ci^B) 
with respect to cosmic neutrinos is expected to lead to seizable absorption dips in 
the diffuse EHEC;/ fluxes to be detected at earth (cf. fig. 2). Independent of the 
nature of neutrino masses, their resolution would constitute the most direct evi- 
dence of the Ci'B so far. In addition, since the annihilation process is sensitive to 
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Fig. 2. a), b) (first two columns): Flavor survival probability Pi^ = Pv^, a = e,fi,T, without and 
with inclusion of thermal background eff'ects, respectively, integrated back to Zs = 5; c) (third 
column) Normalized sum Pv = ^ ^ Pi^a, including thermal effects for Zs = 5, = 20 and Zs = 50 
from top to bottom. All panels assume a normal neutrino mass hierarchy with m^g ^ = 10~^ eV 
and varying (solid) and constant (dotted) neutrino masses as a function of their energy Eq at 
earth. 



the neutrino mass and thus to its possible time variation, it could serve as a test 
for Neutrino Dark Energy (MaVaNs). 

In fig. 2 the survival probability of an EHEC;^ both for varying and constant 
neutrino masses is plotted as a function of its energy Eq as measured at earth. 
It encodes the physical information on possible annihilation processes (for details 
see [2]) and takes values between and 1. In order to disentangle the different in- 
fluences on the EHECz/ survival probability, let us first assume the relic neutrinos 
to be at rest (cf. fig. 2a). In this case the absorption features for constant neutrino 
masses are only subject to the effect of cosmic redshift which causes an EHECi' 
emitted at its source at Zg with energy E to arrive at earth with the red-shifted 
energy Eq = E / {1 + z^) < E. As a result, the survival probability is reduced in the 
interval + Zg) < E^ < E'q'^I and the absorption minima arc thus redshift 

distorted. In case of a possible mass variation, however, the survival probability 
exhibits sharp spikes at the resonance energies E^^f (cf. fig. 2a). This can be un- 
derstood by recalling the mass behavior in the low redshift regime stated in eq. 5. 
The neutrino masses m,^. (z) introduce a z dependence into the resonance energies, 
£;™«(z) = M|(l + z)/{2'm^„ J = £^5^1(1 + z) which is compensated by the effect of 
cosmic redshift. Accordingly, all annihilations occurring at < z < Zg contribute 
to an absorption dip at £'[°*'(z)/(l + z) = -Eq"^. As a second step let us take into 
account that the relic neutrinos are moving targets with a Fermi-Dirac momentum 
distribution. In fig. 2b the corresponding Fermi-smearing both for constant and 
varying neutrino masses results in a thermal broadening (and thus merging) of the 
dips produced by the mass eigcnstatcs (cf. fig. 2a). 

The plot of the fiavor summed survival probability in fig. 2c shows that the dip 
depth increases with Zs independent of the nature of neutrino masses. However, 
generically, the MaVaN features are clearly shifted to higher energies and the min- 
ima deeper than for neutrinos with constant mass. 
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Fig. 3. EgF with F = ^ _F^^ + ^ F,^^ for a flux saturating today's observational bound for 
varying (solid) and constant (dashed) neutrino masses and for Zb = 20, assuming a normal neutrino 
mass hierarchy with m,^^ ^ = 10~^ eV. a) Together with projected sensitivity of LOFAR expressed 
in terms of the diffuse neutrino flux per flavor, corresponding to one event per energy decade and 
indicated duration; b) Together with error bars indicating the corresponding statistical accuracy. 

4. Outlook 

The fig. 3 shows the expected EHEC;^ flux from an astrophysical source at a red- 
shift Zs — 20 both for constant was well as for varying neutrino masses normalized 
to saturate today's observational bound. In fig. 3a it is plotted together with the 
projected sensitivity of the neutrino telescope LOFAR [10] which corresponds to 
maximally 3500 events detected per energy decade and indicated duration. This 
translates into very small statistical error bars as included in fig. 3b. a blow up of 
fig. 3a. Accordingly, independent of the nature of neutrino masses, the detection of 
absorption dips with LOFAR and thus the most direct evidence for the existence 
of the Gi'B so far can be expected within the next decade. Furthermore, if LOFAR 
achieves a decent energy resolution, the variation of neutrino masses and thus the 
interpretation of the Ci^B as source of Dark Energy could be tested. 
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Abstract 

An exciting prediction of higher dimensional, TeV-scale quantum-gravity models is the 
copious production of microscopic black holes in the interactions of high energy particles. 
Such black holes could be produced at the LHC, as well as in the interaction of high energy 
cosmic neutrinos. In this short review we discuss the current bounds obtained from the 
non- observation of neutrino-induced black holes in neutrino detectors as well as future 
perspectives for their detection. 



1. Introduction 

The possibility that black holes will be produced in the collision of two light 
particles at center-of-mass (cm) energies above the Planck scale is an exciting and 
realistic perspective in the context of theories with 5 ~ D — 4 > 1 large compact [1] 
or warped [2] extra dimensions and a low fundamental Planck scale Mjy > 1 TeV 
characterizing quantum gravity. In these theories one might expect the copious 
production of black holes in high energy collisions at cm energies above Md [3]. 
Correspondingly, the Large Hadron Collider (LHC) expected to have a first pilot 
run in late 2007, may turn into a factory of black holes at which their production 
and evaporation may be studied in detail [4,5] (see [6] for a recent review and a 
more complete list of references). 

However, even before the LHC starts operating, important measurements can 
still be made with high-energy neutrino telescopes. In fact, among the strongest 
constraints today are obtained from a (zero-)measurement with neutrino telescopes! 
Black hole production and subsequent decay in the scattering of ultra-high energy 
cosmic neutrinos can be detected by large scale neutrino-detectors [7,8], such as 
IceCube [9] or RICE [10]. Alternatively, black holes produced in the scattering on 
nucleons in the atmosphere may initiate quasi-horizontal air showers far above the 
Standard Model rate [11-14]. 

The current bounds and the perspective for the future are topic of this short re- 
view. In addition we discuss some experimental signatures of black hole production 
which will allow to discriminate against standard model background. 

2. Black hole creation and decay 

The original cross-section calculations for black hole production follow simple, 
semi-classical reasoning. It is assumed that at trans-Planckian parton-parton cm 
energies squared, s 3> M'^, and at impact parameters smaller than the Schwarzschild 
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Fig. 1. Cross section c^J^, Eq. (2), for black hole production in neutrino-nucleon scattering, for 



Miy = 1 TeV, = 5 TeV, and 5 = 2, 4, 6 extra dimensions (solid lines, from bottom to top). 

Also shown is the Standard Model (SM) charged current (CC) neutrino-nucleon cross section 
(dashed-dotted line). Plot taken from [13] 



radius rs of a (4 + (5)-dinicnsional black hole with mass A/bh 
will form. The cross section for black hole production is then: 

with the Schwarzschild radius, 



a/1, a black hole 
(1) 



1 




Frequently one introduces the ratio x = MbijMD and demands that a; ^ 1 to 
ensure that the semi-classical cross-section estimation is justified. 

The contribution of black hole production to the neutrino-nucleon cross section 
is then obtained as the convolution 



The sum extends over all partons in the nucleon, with parton distribution func- 
tions fi{x,fj.) and factorization scale /i. The results are not very sensitive on the 
parton distribution function or the exact value of the factorization scale [13]. The 
cross-sections are displayed in Fig. 1 for several different extra dimensions and 
Mbh/MD > 5. 

There are many ongoing efforts to make the somewhat simplistic geometric cross- 
section calculation more accurate. These include the incorporation of gray body 
[17], form factor [18] or impact parameter effects [19]. Their discussion is beyond 
the scope of this review, however, with a few exceptions they do not alter the 
cross-sections by more than a factor of a few (see [6] for further reading.) 

Once a TeV-mass black hole is created, it will decay primarily through Hawking 
radiation on extremely short time-scales ^ 10^^'' s, thereby producing a large num- 
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Fig. 2. Predictions of the cosmogenic neutrino flux (the sum over all flavors). Long-dashed (long— 
dashed-dottcd) line: Flux from Rcf. [20] for cosmological evolution parameters 
(777 = 4, 2:max = 4). Solid (dotted) line: Flux from Ref. [21], assuming a maximum energy of 
Eniax = 3 ■ lO'^'^''^^' eV for the ultra-high energy cosmic rays. Shaded: Theoretical upper limit 
of the neutrino flux from "hidden" astrophysical sources that are non-transparent to ultra-high 
energy nucleons [24]. The theoretical upper limit is by now surpassed by experimentally upper 
bounds from AMANDA [25] . Plot taken from [13] 

bcr of 0(20) hard quanta [7,4,5]. Since tlic temperature associated with the 
black hole is so large, T ~ Mbh/-^V, the full spectrum of standard model particles 
will be generated in a rather flavor democratic way. 

3. The Neutrino Beam 

Clearly, an essential ingredient to probe the TeV-gravity scenarios is a neutrino 
beam with high enough neutrino energies (~ 10^ GeV). So far no such flux has 
been detected. However, there are a number of model predictions for astrophysical 
neutrino productions. JVIore or less guaranteed, and therefore comprising a reason- 
able lower bound, are the so-called cosmogenic neutrinos which are produced when 
ultra- high energy cosmic protons inelastically scatter off the cosmic microwave back- 
ground radiation in processes such as p7 — > A — > n7r+, where the produced pion 
subsequently decays. Estimates of these fluxes can be found in Refs. [20-23], some 
of which are shown in Fig. 2. While these model predictions might be viewed as a 
lower bound, a theoretical upper bound [24] is also included in Fig. 2. 

4. Current Constraints 

Several dedicated neutrino telescopes and air shower detectors look for high- 
energy neutrinos with energies above 10 PeV. Concerning black hole production it 
was shown [11,13,14] that already now sensible constraints can be obtained from the 
non-observation of horizontal showers by the Fly's Eye collaboration [15] and the 
Akeno Giant Air Shower Array (AGASA) collaboration [16]. These constraints turn 
out to be competitive with other currently available constraints on TeV-scale gravity 
which are mainly based on interactions associated with Kaluza-Klein gravitons, 
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Fig. 3. Lower bounds on Mo as a function of the ratio of the minimum black hole formation 
threshold, MJ^", to M^,. The upper curve for each flux model is the lower bound when the geo- 
metrical cross-section of Eq. 1 is used for the black hole cross section. The lower curve corresponds 
to the lower bound once impact parameter effects arc included. ESS, KKSS and PJ refer to Refs. 
[22], [21] and [23], respectively. Plot taken from [26]. 

according to which a fundamental Planck scale as low as Mu =0(1) TeV is still 
allowed for S > 5 [14]. Constraints from one year operation of the Pierre- Auger- 
Observatory should allow to constrain Mo up to ~ 2 TeV [11,13,14] assuming the 
geometrical cross-section to hold down to M^^"^ ~ Md and a cosmogenic neutrino 
flux according to the more conservative model of [21]. 

Dedicated neutrino telescopes such as IceCube [9] or RICE [10] have other means 
to search for ultra high-energy neutrino interactions. At the energies relevant for 
black hole production the Earth is opaque to neutrinos, hence one only expects 
neutrinos from above the horizon. Due to the high energy, the background of 
atmospheric muons becomes insignificant and one can search for neutrino-induced 
through-going muons as well as for neutrino-induces showers in or near the detec- 
tor. Details of the signature are discussed in the next section. It was estimated that 
with a km'^ sized neutrino detector constraints of the order of a Mo ^ 2 TeV can be 
obtained, both from contained events as well as through-going muon events [7,8]. 
Recently the RICE collaboration has reported the results from a search for very 
high energy neutrinos using a detector of 20 radio antennas distributed within an 
instrumented volume of just ^ 0.01 km'' [10]. However, due to the km-long attenu- 
ation length of the radio signal in cold polar ice, the effective volume is significantly 
larger than the instrumented volume. Already at 1 EeV energies, the volume for 
hadronic cascades exceeds one cubic kilometer. The large volume, combined with an 
analyzed live-time of nearly two years and the fact that no neutrino-induced event 
candidate was found [10] was used to put constrain on the fundamental Planck 
mass Md [26]. These constraints arc shown in Fig. 3 for three different cosmogenic 
neutrino flux models and two different calculations of the black hole production 
cross-section. Assuming the geometrical cross-section calculation, A/bh/A/£) > 1 
and a flux from [21] they obtain a bound on Me, of 3 TeV. The bound weakens 
to about 1.5 TeV, once the cross-section calculation includes effects due to the im- 
pact parameter [19]. Note that these results have only a small dependence on the 
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number of extra dimensions. If the model predictions for cosmogenic neutrino are 
to be believed, these bounds on the fundamental mass scale are currently the most 
restrictive available for (5 > 5. 

5. The Future 

Several future projects are currently planned which will have orders of magnitude 
larger sensitivities for neutrinos of ultra-high energies. Radio detectors such as 
ANITA [27] or SALSA [28], or a hybrid extension of IceCube with radio and acoustic 
sensors [29] will have effective volumes of 0(100 km^), and will therefore be capable 
to detect dozens events due to standard model interactions of cosmogenic neutrinos. 

Individual black hole events will consist of a large hadronic component, a some- 
what smaller electro-magnetic component, as well as muons and taus produced in 
the event. The sensitivity to probe non-standard model interactions will depend 
on observable signatures to discriminate against the standard model background. 
While it is a challenge for sparse high energy neutrino detectors to identify the rel- 
ative contributions, some promising ideas exist nevertheless. The LPM effect will 
lead to a significant stochasticly variable elongation of electro-magnetic showers 
above ~ 10^ GeV energies, while it has less of an influence for hadronic cascades. If 
simultaneously intensity and longitudinal profile can be measured, one might dis- 
criminate the two contributions. The contribution of taus can be identified through 
its decay — a second shower displaced relative to the primary vertex [30]. The pos- 
sibility to identify "soft" muons from a black hole event was studied in [31]. In a 
large hybrid detector which combines optical with radio and acoustic sensitivity, all 
the above signatures could be used for identifying black hole events. One can also 
test the statistical properties of an observed event distribution. Clearly, a sudden 
rise in the energy spectrum will be a first indication of additional non-standard 
cross-sections. An essentially model independent signature of additional interac- 
tion channels would be an observed suppression of the angular distribution of high 
energy neutrino events towards the horizon [32]. 

Summarizing it can be said, that while current bounds from neutrino telescopes 
on Mu arc already probing an important mass-regime, detectors currently under 
construction have the sensitivity to improve the bounds significantly, or better yet, 
discover the first micro black holes. 
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Abstract 

In this talk I review the potential of Icecube for revealing physics beyond the standard 
model in the oscillation of atmospheric neutrinos [1]. 



1. Introduction 

With its high statistics data [2] Supcr-Kamiokandc (SK) cstabhshed beyond 
doubt that the observed deficit in the fi-\ike atmospheric events is due to oscillations, 
a result supported by the K2K and MINOS long-baseline (LBL) experiments [3,4]. 

Mass oscillations are not the only possible mechanism for atmospheric (ATM) 
i/f^ — > Vr flavour transitions. These can be also generated by a variety of nonstan- 
dard physics characterized by the presence of an unconventional v interaction that 
mixes neutrino flavours [5] . Examples include violations of the equivalence principle 
(VEP) [6,7], non-standard neutrino interactions with matter [8], neutrino couplings 
to space-time torsion fields [9], violations of Lorentz invariance (VLI) [10] and of 
CPT symmetry [11,12]. In contrast to the E energy dependence of the conventional 
oscillation length, new physics can produce neutrino oscillations with wavelengths 
that are constant or decrease with energy. [13,14]. At present these scenarios can- 
not explain the data[15] and a combined analysis of the ATM and LBL data can 
be performed to constraint them even as subdominant oscillation effects [16]. 

IceCube, with energy reach in the 0.1 ~ lO^TeV range for ATM neutrinos, is 
the ideal experiment to search for new physics. For most of this energy interval 
standard Arn^ oscillations are suppressed and therefore the observation of an an- 
gular distortion of the ATM neutrino flux or its energy dependence provide a clear 
signature for the presence of new physics mixing neutrino flavours [1] . 

2. Propagation in Matter of High Energy Oscillating Neutrinos 

We concentrate on v^j,-Vr flavour mixing mechanisms for which the propagation 
of v^s {+) and P's (— ) is governed by the following Hamiltonian [12]: 

Am^ is the mass-squared difference between the two neutrino mass eigenstates, 
accounts for a possible relative sign of the new physics (NP) effects between v^s and 
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v^s and A(5„ parametrizes the size of the NP terms. By we denote the possible 
non- vanishing relative phases. 

/ cos6' sin6'\ _ / cos^„ sin^„e=^*''"\ , , 

'~\-s\ne cose)' ~ lv-sin^„e™. cose„ )' ^' 

If NP strength is constant along the neutrino trajectory the oscillation probabil- 
ities take the form [12]: 

P.,^., = 1 - P.,^.. - 1 - sin^ 29 sin^ {^^^ U 
sin^ 28 = -ij (sin^ 261 + RI sin^ 2^™ + 2i?„ sin 261 sin 2e„c?7„) , 
n = Vl + P^ + 2i?„ (cos 26* cos 2^„ + sin 26* sin 2^„c7/„) , i?„ = cr, 



2 Am2 ' 

where, for simplicity, we have assumed scenarios with one NP source characterized 
by a unique A(5„ . c?7„ = cos ?7„ 

Eq. 1 describes, for example, flavour mixing due to new tensor-like interactions 
for which n — 1 leading to a contribution to the oscillation wavelength inversely 
proportional to the neutrino energy. This is the case for j/^'s and i^t-'s of different 
masses in the presence of violation of the equivalence principle due to non- universal 
coupling of the neutrinos, 71 7^ 72 to the gravitational potential [6], so A61 = 
2|0|(7i - 72) = 2|(/)|A7. lyi and 1^2 are related to Vfj^ and v-r by a rotation ^1 — £,vep- 

For constant potential 0, this mechanism is phenomenologically equivalent to the 
breakdown of Lorentz invariance resulting from different asymptotic values of the 
velocity of the neutrinos, ci C2, ASi = (ci — C2) = 5c/c, with i/i and 1^2 being 
related to and by a rotation ^1 = (^yu [10]. 

For vector-like interactions, n = 0, the oscillation wavelength is energy- independent. 
This may arise, for instance, from a non-universal coupling of the neutrinos, fci 7^ k2 
so A5q = Q{ki — ^2) (j^i and V2 is related to the i^^ and i>r by a rotation ^0 ~ S.q), to 
a space-time torsion field Q [9]. Violation of CPT resulting from Lorentz-violating 
effects such as the operator, i^t^'^^^fj.i'^, also leads to an energy independent contri- 
bution to the oscillation wavelength [11,12] which is a function of the eigenvalues of 
the Lorentz violating CPT-odd operator, 6^, ASq = bi —62, and the rotation angle, 
Co = CfZPT, between the corresponding eigenstates Vi and the flavour states Va- 

For most of the neutrino energies considered here, Am^ oscillations are suppressed 
and the NP effect is directly observed. Thus the results will be independent of the 
phase rjn and we can chose the NP parameters in the range A(5„ > 0, < ^„ < 7r/4. 

The Hamiltonian of Eq. (1) describes the coherent evolution of the I'^-i^r ensem- 
ble for any neutrino energy. High-energy neutrinos propagating in the Earth can 
also interact inelastically with the Earth matter either by charged current (CC) and 
neutral current (NC) and as a consequence the neutrino flux is attenuated. This 
attenuation is qualitatively and quantitatively different for j^t-'s and v^'s. j/^'s are 
absorbed by CC interactions while j^t-'s are regenerated because they produce a r 
that decays into another tau neutrino before losing energy [17]. As a consequence, 
for each Vr lost in CC interactions, another i>r appears (degraded in energy) from 
the T decay and the Earth never becomes opaque to i/'^s. Furthermore, a secondary 
flux of P^'s is also generated in the Icptonic decay r —^ nVpiVr [18]. 
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Attenuation and regeneration effects of incoherent neutrino fluxes can be consis- 
tently described by a set of coupled partial integro-differential cascade equations 
(see for example [19] and references therein) or by a Monte Carlo simulation of the 
neutrino propagation in matter [17,18,26]. For astrophysical r^'s, because of the long 
distance traveled from the source, the oscillations average out and at the Earth the 
neutrinos can be treated as an incoherent superposition of mass eigenstates. 

For ATM neutrinos this is not the case because oscillation, attenuation, and 
regeneration effects occur simultaneously when the neutrino beam travels across 
the Earth's matter. For conventional neutrino oscillations this fact can be ignored 
because the neutrino energies covered by current experiments are low enough for 
attenuation and regeneration effects to be negligible. But for non-standard scenario 
oscillations, future experiments probe high-energy neutrinos for which the attenu- 
ation and regeneration effects have to be accounted for simultaneously. 

In order to do so it is convenient to use the density matrix formalism to describe 
neutrino flavour oscillations. The evolution of the neutrino ensemble is determined 
by the Liouville equation for the density matrix p{t) ~ v{t) ® v{t)^ 

f = -^[H,p], (3) 

where H is given by Eq. (1). The survival probability is given by P^^ (i) = Tr[n^^_ 
where 11^^ — v^j,®i'^j, is the state projector, and with initial condition p(0) = H^^^ . 
An equivalent equation can be written for the v density matrix. 

In this formalism attenuation effects due to CC and NC interactions can be 
introduced by relaxing the condition Tr(/9) = 1. In this case 

'j^^-mE),,(EM-^^^^{^,,^(E.t)) , (4) 

whore |Ag„(E,tr" s[A5^(B,f)]-i+|ANc(E.ir'.[ASc(B,'r' -nrW^ScW, 

and [Anc(£'5 t)] ^ — nxix) aj^c{E) {nxix) is the number density of nucleons at the 
point X ^ ct). 

Vt regeneration and neutrino energy degradation can be accounted for by coupling 
these equations to the shower equations for the r flux, Fr {Er , t) (we denote by F 
the differential fluxes dcf) / {dE d cos 6)) . For convenience we define the neutrino Hux 
density matrix F^{E,x) = F^^(E,xq)p{E,x = ct) where F^^(E,xo) is the initial 
neutrino flux: 

^Ii^iE^ ^ -,[H, FAE., x)] - Ij, , {n„, P.(i?., x)} 

dx 2Xf^^{E^,x) 

1 F (E' ,) ^^Nc(i^-,i^.) ^^, 



\nc{EI,x) ' dE, 

1 .p^^E..x)'-^^^^^^^dE.^. 



^Scc(^T,a;) dE, 



dFr{Er,t) _ 1 



■Fr{Er,x) 



' -Tr[U.FAE.,t)]^^^^^^dE.. (6) 



Er ^Cc(-^'''*) dEr 
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X^^^{ErT x) ~ ^rCTr- Tj- is thc T lifetime and 7^ = Er/rrir is its gamma factor. 

dN^c(El.E„) ^ 1 dauc(El,E^) , dNcc{E^,E^) ^ 1 dal,c{E^,E^) , 
dEu — crNc(-E^) rf-B„ '^"'^ — o-J;c(-E^>-) "i-E^ 

easily computed. The r decay distributions d-^di^i^^'^") ^^fj dNi^^{E^,E^) 
found in Refs. [19,27]. 

The third term in Eq. (5) represents the neutrino regeneration by NC interactions 
and the fourth term represents the contribution from Vr regeneration, Vr ^ ^ 
u-r, describing the energy degradation in the process. The secondary flux from 9t 
regeneration, 9^ ^ ^ fj,'^ i/^, is described by the last term where we denote 
by over-bar the energies and fluxes of the r"*". Br^ = 0.18 is the branching ratio for 
this decay. In Eq. (6) the first term gives the loss of taus due to decay and the last 
term gives the r generation due to CC v-r interactions. In writing these equations we 
have neglected the tau energy loss, which is only relevant at much higher energies. 

An equivalent set of equations can be written for the P flux density matrix and 
for the r+ flux. Both sets are coupled due to the secondary v flux term. 

We solve this set of ten coupled evolution equations that describe propagation 
through the Earth numerically using the matter density profile of the Preliminary 
Reference Earth Model and obtain the neutrino fluxes in the vicinity of the detector 
%^=Tr[F^iE,L = 2Rcose)n^] . 

In Fig. 1 we illustrate the interplay between the different terms in Eqs. (5) and (6). 
The figure covers the example of VLI-induced oscillations with Sc/c = 10"^'' and 
maximal £,vii mixing. The upper panels show the final Vfj, and i^t fluxes for vertically 
upgoing neutrinos after traveling the full length of the Earth for the initial condi- 
tions d^{v^,)o/dEi, = d^{v^)o/dE^ cx E^^ and d^{vr)o/dEi, = d^{Dr)o/dE^ = 0. 

The figure illustrates that the attenuation in the Earth suppresses the neutrino 
fluxes at higher energies. The effect of the attenuation in the absence of oscillations 
is given by the dotted thin line in the left panel. Even in the presence of oscillations 
this effect can be well described by an overall exponential suppression [27,24] both 
for i/f^'s and the oscillated i^r's. In other words, we closely reproduce the curve for 
"oscillation + attenuation" simply by multiplying the initial flux by the oscillation 
probability and an exponential damping factor: 

= - 2i?cos0) cM~X{e){aMEHa-ccim , 

dEd cos o dEd cos & 

(7) 

where X{9) is thc column density of the Earth. 

The main effect of energy degradation by NC interactions (thc third term in 
Eq. (5)) that is not accounted for in the approximation of Eq.(7) is the increase 
of the flux in the oscillation minima (the flux docs not vanish in the minimum) 
because higher energy neutrinos end up with lower energy as a consequence of the 
NC interactions. The difference between the dash-dotted line and the dashed line is 
due to the interplay between the i/^- regeneration effect (fourth term in Eq. (5)) and 
the flavour oscillations. As a consequence of the first effect, we see in the right upper 
panel that the fiux is enhanced because of the regeneration of higher energy i^^ 's, 
Vt{E) t~ ^ Vt-{E' < E), that originated from the oscillation of higher energies 
Vfj.^s. In turn this excess of i^r's produces an excess of z/^'s after oscillation which 
is seen as the difference between the dashed curve and the dash-dotted curve in 
the left upper panel. Finally the secondary effect of Vt regeneration (last term in 
Eq. (5)), Vt{E) ^ t+ ^ a*"*" v^{E' < E), results into the larger z^^ flux (seen in 
the left upper panel as the difference between the dashed and the thick full lines). 
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Fig. 1. Vertically upgoing neutrinos after traveling the full length of the Earth taking into account 
the effects due to VLl oscillations, attenuation in the Earth, Ut regeneration and secondary V-r 
regeneration (see text for details). 

This, in turn, gives an enhancement in the i/r flux after osciUations as seen in the 
right upper paneL 

The lower panels show the final fluxes for an atmospheric-like energy spectrum 
d^{iy^)o/dE^ = d^{Uf,)o/dE^ cx E''^ and d^{ur)o/dE^ = d^{vr)o/dE^ = 0. In 
this case regeneration effects result in the degradation of the neutrino energy and 
the more steeply falling the neutrino energy spectrum, the smaller the contribution 
to the total flux. As a result the final fiuxes can be relatively well described by the 
approximation in Eq.(7). 



3. Example of Physics Reach: VLI-induced Oscillations 

The expected number of induced events at IceCube can be obtained by a 
semianalytical calculation as: 

/I /"OO /"OO /"OO /"OO /"OO 

dcose / / dl / dEl^ / dEl / dE, (8) 

dstt'cos e diff'erential muon neutrino neutrino flux in the vicinity of the de- 

tector after evolution in the Earth matter obtained as described in the previous 
section. We use the neutrino fluxes from Honda [20] extrapolated to match at 
higher energies the fluxes from Volkova [21]. At high energy prompt neutrinos from 
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charm decay are important and it is evaluated for two different models of charm 
production: the recombination quark parton model (RQPM) developed by Bugaev 
et al [22] and the model of Thunman ct al (TIG) [23] that predicts a smaller rate. 
'^dE^ -E") is the differential CC interaction cross section producing a muon of 

energy E'j[. T is the exposure time of the detector. Equivalently, muon events arise 
from i>n interactions that are evaluated by an equation similar to Eq.(9). 

After production with energy E'^, the muon ranges out in the rock and in the 
ice surrounding the detector and looses energy. We denote by F (E'^ , E^'^ , I) the 
function that describes the energy spectrum of the muons arriving at the detector. 
We compute the function F{E^^,E^^,l) by propagating the muons to the detector 
taking into account energy losses due to ionization, bremsstrahlung. e+e~ pair 
production and nuclear interactions according to Ref. [24]. 

The details of the detector are encoded in the effective area ^g^j for which we 
make a phenomenological parametrization to simulate the response of the IceCube 
detector after events that are not neutrinos have been rejected ( referred to as "level 
2" cuts in Ref. [25]). The explicit form of cn be found in Ref.[l]. 

Together with z/^-induced muon events, oscillations also generate /x events from 
the CC interactions of the Vr flux which reaches the detector producing a r that 
subsequently decays as r ^ ^iVf^Vr and produces a in the detector: 



dcosO 



dlL 



dl 



dE 



fin 



dE^. 



oo 



uifir 



dEr / dE^ (9) 



-Ex 



where — "^""dE" ' '^^^ found in Ref. [27]. Equivalently we compute the number 
of j^T-induced muon events. 

Neutrino oscillations introduced by NP effects result in an energy dependent 
distortion of the zenith angle distribution of ATM muon events. We quantify this 
effect in IceCube by evaluating the expected angular and Ej^"'^ distributions in the 
detector using Eqs. (9) and (10) in conjunction with the ffuxcs obtained after 
evolution in the Earth for different sets of NP oscillation parameters. 

For illustration we concentrate on oscillations resulting from VLI that lead to 
an oscillation wavelength inversely proportional to the neutrino energy. The results 
can be directly applied to oscillations due to VEP. We show in Fig. 2 the zenith 
angle distributions for muon induced events for different values of the VLI parameter 
Sc/c and maximal mixing (^yu = 7r/4 for different threshold energy E"^™ > ^threshold 
normalized to the expectations for pure Am^ oscillations. The full lines include both 
the I'^-induccd events (Eq.(9)) and i^T-induced events (Eq.(lO)) while the last ones 
are not included in the dashed curves. We see that for a given value of 5c/c there is 
a range of energy for which the angular distortion is maximal. Above that energy, 
the oscillations average out and result in a constant suppression of the number of 
events. Inclusion of the i/7-induced events leads to an overall increase of the event 
rate but slightly reduces the angular distortion. 

In order to quantify the energy-dependent angular distortion we define the vertical- 
to-horizontal double ratio 
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Fig. 2. Upper panels: Zenith angle distributions for muon induced events for different val- 
ues of tlic VLI parameter <5c/c and maximal mixing = n/A for different threshold energy 
E^"^ > ^threshold normalized to the expectations for pure Am^ oscillations . Lower panels: 
The predicted horizontal-to- vertical double ratio in Eq. (10) for different values of Sc/c. The data 
points in the figure show the expected statistical error corresponding to the observation of no NP 
effects in 10 years of IceCube. 



N;^"{eI'"-\ -0.6 < cosg < -0.2) 
- Ror f 7Vr"''''(-E,{'"'',-0.6 < cose* < -0.2) 

NJ^°~'''^'{El'"'\-l < cos9 < -0.6) 

where by we denote integration in an energy bin of width 0.2 logio(-BJ^'"'') 

using that IceCube measures energy to 20% in logj^Q E for muons. 

In what foUows we wiU use the double ratio in Eq. (10) as the observable to 
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Fig. 3. Sensitivity limits in the (5c/c, gvli '^^ 90i 95, 99 and 3 a CL. The hatched area in the upper 
right corner is the present 3cr bound from the analysis of SK data in Rcf. [16]. 

determine the sensitivity of IceCubc to NP-induced oscillations. We have chosen 
a double ratio to eliminate uncertainties associated with the overall normalization 
of the ATM fluxes at high energies. It is worth noticing that using this observable 
relies on the fact that the zenith angular dependence of the effective area is well 
understood. 

In Fig. 2 we plot the expected value of this ratio for different values of Sc/c. As 
mentioned above, IceCube measures energy to 20% in log^p E for muons. Accord- 
ingly, we have divided the data in 16 i?^'" bins: 15 bins between 10^ and 10'^ GeV 
and one containing all events above 10^ GeV. In the figure the full lines include 
both the i/^-induced events (Eq.(9)) and ;/T--induced events (Eq.(lO)) while the last 
ones are not included in the dashed curves. As described above, the net result of 
including the z^r-induced events is a slight decrease of the maximum expected value 
of the double ratio. The data points in the figure show the expected statistical error 
corresponding to the observation of no NP effects in 10 years of IceCube. In order 
to estimate the expected sensitivity we assume that no NP effect is observed and 
define a simple function including only the statistical errors. 

We show in Fig. 3 the sensitivity limits in the [(5c/c, ^vii]-plane at 90, 95, 99 and 
3 a CL obtained from the condition x^((5c/c, ^vii) < Xmaa: (^L, 2dof). We show in 
the figure the results obtained using the RQPM model and the TIG model. The 
difference is about 50% in the strongest bound on Sc/c. The figure illustrates the 
improvement on the present bounds by more than two orders of magnitude even 
within the context of this very conservative analysis. The loss of sensitivity at 
large Sc/c is a consequence of the use of a double ratio as an observable. Such an 
observable is insensitive to NP effects if Sc/c is large enough for the oscillations to 
be always averaged leading only to an overall suppression. 

When data becomes available a more realistic analysis is likely to lead to a further 
improvement of the sensitivity. 
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Abstract 

Several phenomenological models of physics beyond the Standard Model predict flavor 
mixing in the neutrino sector in addition to conventional mass-induced oscillations. In 
particular, violation of Lorentz invariance (VLI) results in neutrino oscillation effects 
parametrized by the maximal attainable velocity difference 5c/ c. We report on a study of 
the sensitivity of the AMANDA-II detector to such effects using distortions in the spec- 
trum of high-energy atmospheric neutrinos. For maximal mixing and six years of simulated 
data, the preliminary sensitivity of AMANDA-II to VLI of this type is 5c/c < 2.1 x 10"^'' 
at the 90% confidence level. 



1. Introduction 

Flavor oscillations in the neutrino sector provide an interesting method to test 
phenomenological models of physics beyond the Standard Model. While mass- 
induced oscillations of atmospheric neutrinos are on firm experimental footing [1-3] , 
subdominant effects may yet be present. In particular, violation of Lorentz invari- 
ance (VLI) can result in oscillations at high energies and can distort the atmospheric 
neutrino spectrum. 

The AMANDA-II detector, a subdetcctor of the IceCube experiment, is an array 
of 677 optical modules buried in the ice at the geographic South Pole which detects 
the Cerenkov radiation from charged particles produced in neutrino interactions 
with matter [4]. In particular, muons produced in charged-current and inter- 
actions deposit light in the detector with a track-like topology, allowing us to use 
directional reconstruction to reject the large background of down-going atmospheric 
muon events. After suitable quality selection criteria arc applied, AMANDA-II ac- 
cumulates atmospheric neutrino candidates above 50 GeV at a rate of 4 per day 
[5]. While conventional oscillations arc suppressed at these energies, VLI effects can 
be detected or constrained by their influence on the zenith angle distribution and 
energy-correlated observables. 

2. Phenomenology 

Various new physics scenarios can result in neutrino flavor mixing beyond con- 
ventional oscillations. We focus here on oscillations induced by differing maximally 
attainable velocities (MAVs) in the neutrino sector. MAV cigenstates can be dis- 
tinct from flavor cigenstates, resulting in oscillations characterized by the MAV 
difference (5c/c = (ci — C2)/c. 
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Conventional and VLI oscillations can be combined in a two-family scenario, with 
the following survival muon neutrino survival probability as a function of energy E 
and baseline L (in energy units) [6-8]: 

P.,^.. = 1 - sin' 26 sin^ f 7^j , (1) 



where 



sin^ 29 = ^ (sin^ 26* + sin^ 2^ + 2i? sin 2B sin 2^ cos x]) , (2) 



and 



cos 2^^ cos 2^ + sin 20 sin 2^ cos 77) , (3) 

Standard oscillations are characterized by the mass-squared difference Am^ and 
mixing angle 0, while VLI oscillation parameters include the velocity difference (Sc/c, 
the mixing angle ^, and the phase 77. If we take both conventional and VLI mixing 
to be maximal {0 = ^ = 7r/4) and set cos 77 = 1, this reduces to the following: 

/ • n 1 -2 /Am^i dcLE\ 
^V^-i.^ (maximal) = 1 - sm I + — — ) ■ (5) 

Note the different energy dependence of the two effects. For atmospheric neu- 
trinos, the zenith angle functions as a surrogate for the baseline L, allowing path 
lengths up to the diameter of the Earth. Figure 1 shows the survival probabil- 
ity as a function of neutrino energy and zenith angle for the maximal case, as in 
equation (5). 



3. Analysis Methodology 

First, to obtain a clean sample of atmospheric neutrinos, we must separate these 
from the large background of atmospheric muons. Selecting events with a recon- 
structed zenith angle below the horizon allows rejection of many such events, but we 
must generally apply further quality criteria to eliminate mis-reconstructed muons. 
For this study, we have used the selection criteria from the 2000-03 AMANDA-II 
point source search [5] and examine only zenith angles > 100°. 

Next, our goal is to measure or constrain the energy-dependent angular distor- 
tions caused by VLI effects. While AMANDA-II has an angular resolution of a few 
degrees [9], reconstruction of the neutrino energy is more difficult and fundamen- 
tally limited by the stochastic losses of the muon. Instead, we use a well-simulated 
energy-correlated observable, the number of triggered optical modules (Nch)- 

Now, to determine values of the parameters 9i of our hypothesis (in the simplest 
one-dimensional case, just 6c/c) that are allowed or excluded at some confidence 
level, we follow the likelihood prescription described by Feldman and Cousins [10]: 

- For each point in the parameter space 9i, we sample many times from the parent 
Monte Carlo distributions of the observable(s) (MC "experiments"). 

- For each MC experiment, we calculate the log likelihood ratio 
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Fig. 1. Atmospheric survival probability as function of neutrino energy and zenith angle. 
Conventional oscillations are present at low energies, wrhile high-energy oscillations are due to VLI 
(maximal mixing, (5c/c = lO"^'^). 

A/: = -21nU + 21nU,be.t , (6) 

where is the Poisson probabihty that the MC experiment is derived from the 
parent distribution at 9i (other hkehhood formulations are possible). 

- For each point 6i, wc find the value /S.Ccrit at which, say, 90% of MC experiments 
have a lower A>C. 

- Finally, we compare the A£ of the data (or in our case, a simulated data set ge- 
nerated under the null hypothesis) with the critical surface l^Ccrit , Sund regions of 
the parameter space at which AC > ACcrit are excluded at that confidence level. 
For a one-dimensional parameter space, this can likely be interpreted an upper 
limit, and one can calculate a median sensitivity by iterating over a number of 
simulated data sets. 

As noted in [10], the likelihood formulation has a number of desirable features 
compared to a standard approach, the most significant being proper coverage. 

4. Sensitivity of AMANDA-II 

We have performed a Monte Carlo study using six years of simulated AMANDA- 
II data: an integrated exposure of 1200 days, approximately 5100 events below the 
horizon under the null hypothesis (conventional oscillations only). For this initial 
study, we have tested only the Nch distribution across a one-dimensional param- 
eter space, varying the VLI strength Sc/c. To anticipate the impact of the inclu- 
sion of systematic errors in the future, we have left free the normalization of the 
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log(5c/c) 

Fig. 2. Likelihood ratio for VLI effects using the shape of the Af^-h distribution, for values of the 
parameter (5c/c. The critical curves for various confidence levels are shown, along with A£ for a 
simulated six- year data set. Values oi Sc/c to the right of the point of intersection with the critical 
curve are excluded. 

atmospheric neutrino flux (i.e. treating it as a nuisance parameter). Wc have not 
included the zenith angle distribution in this analysis, as we have not yet accounted 
for systematic uncertainties in the shape of the spectrum. The curves of ACcrit for 
the 90%, 95%, and 99% confidence levels are shown in Figure 2, along with the 
likelihood ratio for a single simulated data set. 

Assuming maximal mixing (sin 2^ = 1) ^^^d phase cos 77 = 1, we find a median 
sensitivity of Sc/c < 2.1 x 10~^^ at the 90% confidence level. Existing experimental 
hmits include the MACRO result of 5c/c < 2.5 x 10"^^ [11] and the hmit by 
Gonzalez-Garcfa and Maltoni using the Super-Kamiokande + K2K data, Sc/c < 
2.0 X 10-27 [8]. 

5. Conclusions and Outlook 

Using its large sample of atmospheric neutrinos, AMANDA-II is capable of de- 
tecting or constraining high-energy new physics effects in the neutrino sector. The 
Monte Carlo study presented here indicates a sensitivity to VLI effects competi- 
tive with existing limits, and a number of improvements (such as testing multiple 
observables) and optimizations (event selection criteria and the binning of the ob- 
servables) are forthcoming. We anticipate applying this analysis in the near future to 
the AMANDA-II data collected during 2000-2005. Furthermore, the same method- 
ology can also be applied to constrain other physics beyond the Standard Model, 
such as violations of the equivalence principle [13] or quantum decoherence resulting 
from interactions of neutrinos with the background space-time foam [14-16]. 

The next-generation IceCube detector, with an instrumented volume of 1 km^, 
will allow unprecedented sensitivity to these same effects. In 10 years of operation, 
IceCube will collect a sample of over 700 thousand atmospheric neutrinos and will 
be sensitive at the 90% confidence level to VLI effects at the level of Sc/c < 2.0 x 
10^2^ [12]. This high-statistics sample will also provide an opportunity to test other 
phenomenological models of physics beyond the Standard Model. 
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Abstract 

Neutrinos allow for a test of the hypothesis that the fermions of the Standard Model have 
Fermi-point splitting, analogous to the fermionic quasi-particles of certain condensed- 
matter systems. If present, the corresponding Lorentz-violating terms in the Hamiltonian 
may provide a new source of T and CP violation in the leptonic sector, which is not directly 
related to mass. 



1. Introduction 

The basic idea of this talk is to suggest neutrinos as a probe of radically new 
physics. Of course, this is a long-shot . . . but worth trying. 

One example of such new physics would be related to the concept of emergent 
symmetries [1~4]. Lorentz invariance, for example, would not be a fundamental 
symmetry but an emergent phenomenon at low energies. 

In order to be specific, we start from an analogy with quantum phase transitions 
in fermionic atomic gases or superconductors and consider the hypothesis [5-7] 
that the fermions of the Standard Model have tiny Lorentz-violating effects due to 
Fermi-point splitting (abbreviated FPS and explained below). 

If Fermi-point splitting would indeed occur for the quarks and leptons of the 
Standard Model, then neutrinos may provide a unique window to the underlying 
theory [8-11]. Specifically, there could be new effects in neutrino oscillations, pos- 
sibly showing significant T and CP violation (and perhaps even CPT violation). 
The aim of this talk is to sketch some of the potential FPS effects but we refer, in 
particular, to the contribution of M.C. Gonzalez-Garcia in these Proceedings for a 
more general discussion of non-standard neutrino oscillations. 

The outline of this write-up is as follows. In Sec. 2, some background on condensed 
matter physics is given and, in Sec. 3, a possible application to elementary particle 
physics is discussed. In Sec. 4, a simple but explicit neutrino model with both Fermi- 
point splittings and mass differences is introduced. In Sec. 5, some interesting results 
on neutrino oscillations from this model are reviewed. In Sec. 6, concluding remarks 
are presented. 

2. Fermi-point splitting in atomic and condensed-matter systems 

Ultracold quantum gases of fermionic atoms (e.g., ^Li at nano-Kelvin tempera- 
tures) are extremely interesting systems, especially as they can have tunable inter- 
actions by way of magnetic-field Feshbach resonances. In the so-called BEC-BCS 
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crossover region of these systems, a BCS-type condensate has recently been ob- 
served for s-wave pairmg [12]. As usual, BEC stands for Bose-Einstein condensate 
and BCS for the superconductivity triumvirate Bardcen, Cooper, and Schricffer. 

For the BEC-BCS crossover region in systems with p-wave pairing, there is the 
prediction [5,6] that a quantum phase transition between a vacuum state with fully- 
gapped fermionic spectrum and a vacuum state with topologically protected Fermi 
points (gap nodes) occurs. Here, we only give a simple illustration of this new type 
of quantum phase transition and refer the reader to Ref. [13] for an extensive review. 

The Bogoliubov-Nambu Hamiltonian for fermionic quasi-particles in the axial 
state of p-wave pairing is given by 

_/ |p|V(2m)-9 c_lP- (ei +ie2)\ , 
""^^ [c^p■{e,^^e,) -|p|V(2m) + g j ' 

with m the mass of the fermionic atom (considered is the direction of atomic spin, 
which experiences the Feshbach resonance), (ei, 62, 1) an orthonormal triad, 1 the 
direction of the orbital momentum of the pair, c± the maximum transverse speed, 
and q a parameter controlled by the magnetic field near the Feshbach resonance. 
The energy spectrum of this Hamiltonian is readily calculated: 

ElM = (|P|V(2™) - qf + cl |p xT|^ (2) 

Clearly, there are two regimes. For parameter g < 0, on the one hand, there is a 
BEC regime with mass gap, E ^ 0. For parameter g > 0, on the other hand, there 
is a BCS regime with two Fermi points in momentum space, 

bi = +PFl, b2 = -pFl, PF = y2mg, (3) 

at which the energy function vanishes, E{p) = for p = with a = 1,2. 

There is then a quantum phase transition at q = 0, with a mass gap for q < 
and a space-like splitting of Fermi points (Ab = bi — b2 7^ 0) for q > 0; see Fig. 1. 

This example also clarifies the concept of emergent relativity mentioned in the 
Introduction. Consider momenta close to one of the two Fermi points, for example, 
p = bi -|- k with |k|(pF- Then, the energy 2 becomes 

El^ - {PFlmf kl+clkl^l?{h;i-\-kl), (4) 

after the following rescalings: 

fcjj = k • 1 = {cui/pf) fc|| , fcj_ = I k X 1 1 = ( c/c_l) k± , (5) 

which would be appropriate for a local observer made of the same quasi-particles 
[4,13]. In terms of the rescaled momentum k, relation 4 corresponds precisely to 
the mass-shell condition of a massless relativistic particle. 

3. FPS hypothesis for elementary particle physics 

Based on the analogy with certain condensed-matter systems discussed in Sec. 2, 
the following hypothesis has been put forward [5,6]: perhaps the chiral fermions of 
the Standard Model also have Fermi-point splitting (FPS). Specializing to time- 
like splittings (A60 7^ 0) and vanishing Yukawa coupling constants (i.e., vanishing 
fermion masses), their dispersion relations would be given by: 

{E,,j{p)y=[c\p\+b[fjy, (6) 
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Fig. 1. Quantum phase transition at q = qc between a quantum vacuum with mass gap and one 
with topologically-protectcd Fermi points (gap nodes) . At q = qc, there appears a marginal Fermi 
point with topological charge N = (inset at the top). For q > qc, the marginal Fermi point has 
split into two Fermi points characterized by nonzero topological invariants N = ±1 (inset on the 
right). A system described by Hamiltonian 1, for 1 = (0, 0, 1), has critical parameter qc = 0. 

where a labels the 16 types of massless left-handed Weyl fermions (including a 
left-handed antineutrino) and / the iVfam fermion families (henceforth, we take 
-^fam = 3). The maximum velocity of the fermions is assumed to be universal and 
equal to the velocity of light in vacuo, c. Note that we still speak about Fevmi-point 
splitting even though the energy 6 for b^/J < gives rise to a Fermi surface. 
One possible FPS pattern is given by the following factorized Ansatz [6]: 

bii^=YaW\ (7) 

where Ya arc the known hypercharge values of the fermions and 6g three unknown 
energy scales. Independent of the particular FPS pattern, the dispersion relations 
of massless left-handed neutrinos and right-handed antineutrino would be 

(i?.„/(p))'= (c|p|+6(«)', {E,,jip))'^[c\p\+sbif^y, (8) 

where a value s = 1 respects CPT and s = — I violates it. 
More generally, one may consider for large momentum |p|: 

E{p) ^ c IpI ± 6o + m2cV(2c IpI) + 0( l/\p\^ ) . (9) 

The conclusion is then that the search for possible FPS effects prefers neutrinos 
with the highest possible momentum. 

At this point, two questions on energy scales arise. First, what is known experi- 
mentally? The answer is: not very much, apart from the following upper bounds: 

3 

< 1 keV, ^TO/<100eV, (10) 
/=i 

from low-energy neutrino physics [14] and cosmology, respectively. 
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Second; what can be said theoretically about the expected energy scale of FPS? 
The answer is: little to be honest, but perhaps the following speculation may be 
of some value. For definiteness, start from a particular emergent-physics scenario 
with two energy scales [7]: 

- -Blv of the fundamental Lorentz-violating fermionic theory; 

- B'comp as the compositcness scale of the Standard Model gauge bosons. 
Taking the LEP values of the gauge coupling constants, the renormalization-group 
equations for A'^tam = 3 give these numerical values: 

i^comp - 10^^ GeV , Eu, ^ 10'*2 GeV . (11) 

The speculation, now, is that perhaps ultrahigh-energy Lorentz violation re-enters 
at an ultralow energy scale: 

l&ol-^i^c'omp/^^LV-lO-^eV. (12) 
If correct, this motivates the search for FPS effects at the sub-eV level. 



4. Simple FPS neutrino model 

A general neutrino model with both Fermi-point splittings (FPS) and mass dif- 
ferences (MD) has many mixing angles and complex Dirac phases to consider (not 
to mention possible Majorana phases). In order to get an idea of potentially new 
effects, consider a relatively simple FPS-MD neutrino model [10,11] having 

- a standard neutrino mass sector with "optimistic" values for ^13 and S; 

— a FPS sector with large mixing angles, energy splittings, and Dirac phase uj. 
Specifically, the mass sector has the following mass-square-difference ratio, mixing 
angles, and Dirac phase: 

ATOoi 7Tlo — ?71? 1 IT n 1 TT 

Rrr.^-^^—^ ^ = ^, e21=^^32-J, siu^ 2^13 - , <^ - , (13a) 

Ato|2 mi, — 30 4 20 2 

and the FPS sector has energy-difference ratio, mixing angles, and Dirac phase: 

^=T17ii) -7(i)— (i) =1' X21=X32=Xl3=- = -. (13b) 

For later use, we also define two additional models. The first additional model 
is a pure FPS model [9] with trimaximal couplings (x2i = X32 = 7r/4 and xi3 = 
arctan •\/l/2 ), complex phase lo, and FPS ratio R. At sufficiently high energies, the 
model for i? = 1 and w = 7r/4 is close to the FPS-MD model mentioned above. 

The second additional model is a pure MD model with a mass-square-difference 
ratio Rm — 1 /30 and the following more or less realistic values for the mixing angles 
and Dirac phase: sin^ 26*23 — 1- sin^ 20i2 = 0.8, sin^ 26*13 = 0.2, and (5 = 0. 

In the rest of this contribution, these three models will be referred to as the 
FPS-MD model, the FPS model, and the MD model, respectively. 



5. FPS effects in neutrino oscillations 

Consider now a high-energy [E^, ^ c |p|) neutrino beam traveling over a distance 
L. Neutrino oscillations from the FPS-MD model of Sec. 4 are then determined by 
two dimensionless parameters: 

- ^^^""^ . (14.) 



'311^ \20GeVy\ L I \ 
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{ii^,sni-(2iJi.^), f>, sga Atti- ; H,siiL-(2xL^), .j,ssa At^} = {1/30, l/20,7r/2, + , 1, l,7r/4, +} 
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Fig. 2. Vacuum probabilities from the FPS-MD model (13ab) as a function of the dimensionless 
parameters p and t, defined by Eqs. (14ab). Top panels: P = P(t'^ — > fe)- Bottom panels: 
P" = P{ue — > i/fj,). If CPT invariance holds, also P = P{'Ue V^) and P" = — > Fe). Shown 
are constant-r slices, where the heavy-solid curves in the two left panels correspond to r = 
(pure mass-difference model) and the other thin-solid, long-dashed, and short-dashed curves for 
positive r correspond to r = 1, 2, (mod 3), respectively. 
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for numerical values of L and En appropriate to a neutrino factory [15]. Possible 
new effects in neutrino oscillations from FPS may occur as 

- energy dependence of the vacuum mixing angle 0i3 [10]; 

- novel source of T, CP, and perhaps CPT violation [11]; 

- modified flavor ratios for high-energy cosmic neutrinos [8,9]. 
In this contribution, we discuss only the last two effects. 

Figure 2 shows that, provided the FPS parameter A6q ' is large enough for given 
baseline L, the probabilities of time-reversed processes can be different by several 



tens of percents: Pivu 



20 % versus P{ve i^^) « 80 % at p ~ 1 and r ~ 3, 



for example. For the record, standard mass-diffcrcncc neutrino oscillations (r — 0) 
give more or less equal probabilities at p ~ 1: P{v^ Ve) ~ P{ve v^) w 0. 
In short, there could be strong T-violating (and CP-violating) effects at the high- 
energy end of the neutrino spectrum from FPS or other emergent-physics dynamics. 

Next, turn to the pure FPS model and also, for comparison, to the pure MD 
model, both defined in Sec. 4. Pion and neutron sources then give the averaged 
event ratios shown in Table 1, with the clearest difference between the two models 
for the case of a neutron source. In principle, these results may be relevant to 
high-energy cosmic neutrinos but it remains to be seen whether or not present 
experiments (e.g., AMANDA and IceCube) can access this type of information. 
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Table 1 

Averaged event ratios {Ne : Nfj, : Nr) from pion and neutron sources for pure Fermi-point- 
splitting (FPS) and mass-difference (MD) neutrino models as defined in Sec. 4. The MD event 
ratios are taken from Rcf. [16]. 





TT : initial ratios = (1:2:0) 


n : initial ratios =(1:0:0) 


FPS (cj) 
FPS (7r/4) 
MD 


(6:7-1- cos 2ui : 5 — cos 2ui) 
(0.33 : 0.39 : 0.28) 
(0.36 : 0.33 : 0.31) 


(1:1:1) 
(0.33 : 0.33 : 0.33) 
(0.56 : 0.26 : 0.18) 



6. Outlook 

From a phcnomcnological perspective, the Fcrmi-point-splitting (FPS) hypothesis 
suggests the following three research directions: 

- the possible energy dependence of the vacuum mixing angle from FPS, which 
can be tested by neutrino experiments at a superbeam or neutrino factory; 

- the possibility of a new source of leptonic CP violation, which impacts on the 
physics of the early universe (e.g., the creation of baryon and lepton number); 

- the possible modification of the propagation of high-energy cosmic neutrinos by 
FPS effects, which may be of relevance to present and future neutrino telescopes. 

From a more theoretical perspective, the outstanding issues are: 

- the precise nature of the conjectured re-entrance mechanism of Lorentz violation 
at ultralow energy from Lorentz violation at ultrahigh energy (condensed-matter 
physics can perhaps provide some guidance); 

- the explanation of the large hierarchies of basic scales (e.g., for mass or FPS). 
But apart from these theoretical ideas, experiment may, of course, suggest entirely 
different directions . . . 
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Abstract 

Models with extra dimensions and the fundamental scale at the TeV could imply signals 
in large neutrino telescopes due to gravitational scattering of cosmogenic neutrinos in 
the detection volume. Apart from the production of microscopic black holes, extensively 
studied in the literature, we present gravity-mediated interactions at larger distances, that 
can be calculated in the eikonal approximation. In these elastic processes the neutrino loses 
a small fraction of energy to a hadronic shower and keeps going. The event rate of these 
events is higher than that of black hole formation and the signal is distinct: no charged 
leptons and possibly multiple-bang events. 



1. Motivation: cosmogenic neutrinos and TeV gravity 

Cosmogenic neutrinos, produced in the scattering of protons off cosmic microwave 
background pliotons, liave access to TeV physics in interactions with terrestrial 
nucleons at center of mass energies ^/s = v^2mjv^E^ > 10 TeV. If the fundamental 
scale of gravity is ~ 1 TeV [1], which may happen in > 4 space-time 
dimensions, these vN interactions are transplanckian, ^/s > Md- 

The only consistent theory known so far in such a regime, string theory, tells us 
that the interactions arc soft in the ultraviolet. The scattering amplitudes vanish 
except in the forward region, an effect that can be understood as the destructive 
interference of string excitations [2] . The forward amplitudes are dominated by the 
zero mode of the string, corresponding to the exchange of a gauge particle of spin 

1, A ^ gs/t, for open strings, or a graviton of spin 2, A ^ (l/M|,)s^/t, for closed 
strings. Therefore, one expects that gravity dominates in transplanckian collisions. 

It must be noticed that present bounds on Mu from colliders (LEP, Tevatron) 
[3] or astrophysics and cosmology (supernovae cooling) [4] come from processes at 
energies below Md and are indirect, since they actually constrain the energy emit- 
ted to Kaluza-Klein gravitons of mass M oc in the n = D — A compact extra 
dimensions, which is a function of the compactification radius R. Those bounds rely 
on the assumption of all extra dimensions being large (the effective and fundamental 
Planck scales then relate through Mp oc R^M^"") that can be evaded in more so- 
phisticated compactification models [5]. In contrast, transplanckian collisions probe 
Md directly and independently of compactification details. 

2. Gravitational interactions 

We have shown that gravitational interactions arc the only relevant in the trans- 
planckian regime of energies. In impact parameter space, one must keep in mind 
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two critical values: the Planck length Xd ^ ^^d^ ^-^^d the Schwarzschild radius 
Rs{s) ~ i^/s/Moy^^"''^^'' . There are two types of interactions. 

Short-distance interactions, with impact parameter b < Rs, in which the colliding 
particles (a neutrino and a parton inside the nucleon) collapse into a black hole (BH) 
correspond to the exchange of strongly coupled gravitons of high momentum (non- 
linear gravity). The collapse involves strongly coupled gravity and is not calculable 
perturbatively. Most analyses are based on a geometric cross section [6] ctbh — 
7ri?|(s) for the partonic process, with s = xs. If a/I ^ Mjj, namely Rs ^ -^^_d^j 
one expects that this estimate will not be off by any large factors [7]. However, most 
of the BHs produced in the scattering of an ultrahigh energy neutrino off a parton 
are light, with masses just above Md, since the uN cross section is dominated 
by the low x region. In this regime the amount of gravitational radiation emitted 
during the collapse or the topology of the singularity are important effects that add 
uncertainty to the geometric estimate. 

Long-distance interactions, with 6 3> Rs, have to do with the exchange of weakly 
coupled gravitons of low momentum (linearized gravity) [8] . In transplanckian col- 
lisions quantum gravity acts inside the event horizon [Rs > Xd)- Therefore, these 
elastic interactions are due to classical gravity. They are characterized by a small 
deflection angle 9* in the center of mass (CM) frame, 

r ^ ( —] ^ a^v = cf/s = -n - cosr)(i . (i) 

The elastic collision of a neutrino and a parton that exchange D-dimensional gravi- 
tons is then described by the eikonal amplitude resumming an infinite set of ladder 
and cross-ladder diagrams in the limit in which the momentum q carried by each 
graviton is smaller than the CM energy or, in terms of the fraction of energy lost 
by the incoming neutrino, y = (E,^ — E'^)/ E^{1. In this limit the amplitude is inde- 
pendent of the spin of the colliding particles. Essentially, is the exponentiation 
of the Born amplitude in impact parameter space [9] : 

Aeik(s, = y y e^'i-b (e'^^'^"^ - l) = ATTsblF,,{hq) , (2) 
where x{s, b) is the eikonal phase, 

and a new scale be appears, 
■(47r)t-i /7i 



6c(s) 



-r 



2 / M 



D 



(4) 



The total partonic cross sections can be obtained from the amplitudes above using 
the optical theorem. They are i7eik oc &^ ~ s~, growing faster with energy than the 
BH cross sections ctbh i?| ~ S"+t. 

Therefore, in transplanckian collisions one may consider two types of processes 
[10]: elastic (long-distance) soft processes where the neutrino transfers to the par- 
tons a small fraction y < ymax of its energy and keeps going, and shorter distance 
(6 < Rs) hard processes where the neutrino loses in the collision most of its energy. 
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Fig. 1. Probabilities defined in the text for a 10 GeV neutrino reaching IccCubc from 9z = a,s 
a function of Mo for n = 2 and n = 6. 



possibly collapsing into a BH. Wc take j/max — 0.2, the typical inelasticity of a 
standard model (SM) interaction, but any value of the order of 0.1 yields similar 
results. On the other hand, there is a ?/min = E^hxes/ E^, determined by the threshold 
energy i?thrcs transfered to a parton that produces an observable hadronic cascade. 
The corresponding vN cross sections are obtained by convolution with the parton 
distribution functions f{q^q^g}{x, ^) with the appropriate energy scale ^ [9,10]. 

To estimate the relative frequency of both type of processes [10], consider a 
10^° GeV neutrino that scatters off a nucleon with E'tiucs = 100 TeV and M^j = 
1 TeV for n = 2 (6) extra dimensions. The number of eikonal interactions before the 
neutrino gets destroyed is the ratio of interaction lengths isnZ-^cik = cTcik/cBH = 
12.5 (1.64). For a SM interaction isM = 440 km while Lbh = 17 km (4 km) in ice. 
The total energy lost by the neutrino in these eikonal interactions and the energy 
lost to graviton radiation are relatively small: E^°^^ = 5.9 x 10'' GeV (1.2 x 10^ GeV) 
and El^g = 9.2 x 10^ GeV (1.2 x 10^ GeV) in 1 km of ice. 



3. Signals at neutrino telescopes 

The flux of cosmogenic neutrinos, yet unobserved, depends on the production rate 
of primary nucleons of energy around and above the GZK cutoff. It is correlated with 
proton and photon fluxes that must be consistent, respectively, with the number 
of ultrahigh energy events at AGASA and HiRes [11] and with the diffuse 7-ray 
background measured by EGRET [12]. We base our analysis on the two neutrino 
fluxes described in [13]. The first one saturates the observations by EGRET, whereas 
for the second one the correlated flux of 7-rays contribute only a 20% to the data, 
with the nucleon flux normalized in both cases to AGASA/HiRes. The higher flux 
predicts 820 downward neutrinos of each flavor with energy between 10^ GeV and 
10^^ GcV per year and km^, versus 370 for the lower one. The spectrum has a peak 
at neutrino energies between 10^ GeV and IQ^^ GeV. 

When a neutrino hits a nucleon it will start a hadronic shower. The total number 
of hadronic events in a neutrino telescope of cross sectional area A in a time T is 

iVevcnts = 27rAT j dE, ^ J d COS 6 , P^urv Pint , (5) 

where Pgurv is the probability that the neutrino survives to reach the detector and 
Pint the probability that it interacts inside the detector: 
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Fig. 2. Energy distribution (events per bin) of 
the eikonal (solid), BH (dashed) and SM (dot- 
ted) events in IceCube per year for the higher 
and the lower cosmogenic fluxes, Mj) = 2 TeV 
and n = 2, 6. 



Fig. 3. Contained events per year in IceCube 
and AMANDA for the higher (thick) and the 
lower (thin) cosmogenic fluxes and n = 2, 6. 
We show eikonal (solid), multi-bang (dashed— 
dotted) and BH (dashed) events. 



Pint(i?.)«l-e-^'''-^-<- , (6) 



with X the column density of material, 9^ the zenith angle and L the longitudi- 
nal detector size. When L is larger than the interaction length Lq, there may be 
multiple-bang events. Neglecting the energy lost by the neutrino in each interaction, 
the probability of N bangs and the average, and most probable, number of bangs 
are, respectively 

P^(L)=e-V^"i^^, (iV) = £ 7VP„ = L/Lo . (7) 

The different probabilities for a typical cosmogenic neutrino of 10^" GeV reaching 
vertically IceCube are shown in Fig. 1 for illustration. Double-bang events could also 
be produced by SM interactions (the decay of a tau created in a first interaction) 
or in the BH evaporation. For the double-bang tau event to be contained inside 
a detector like IceCube (1 km of length with 125 m between strings), the energy 
of the tau lepton must be between 2.5 x 10^ GeV and 10^ GeV. In this case, the 
probabihty is only 6.8 x 10^^. 

The energy distribution of the hadronic cascades and the total number of black 
hole and eikonal events at AMANDA (0.03 km^ and a length of 700 m) and Ice- 
Cube (1 km'^) for the neutrino fluxes introduced above are given in Figs. 2 and 3, 
respectively. In the SM we expect 1.32 (0.50) contained events per year in IceCube 
for the higher (lower) flux. Of those, 0.38 (0.14) would come from a neutral current 
and 0.94 (0.36) from a charged current. 
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4. Conclusions 



Cosmogenic neutrinos directly probe TeV gravity in transplanckian collisions 
with nucleons at Earth if there are D > A space-time dimensions. Two types of 
interactions take place. Hard processes, when the impact parameter is smaller than 
the Schwarzschild radius of the neutrino-parton system, produce mostly light black 
holes with theoretically uncertain cross sections. At larger distances, soft elastic 
processes occur in which the neutrinos lose a small fraction of energy to a hadronic 
shower, in a well known regime described by the eikonal approximation. 

The latter turn out to be dominant and produce a clear signal in large neutrino 
telescopes: contained hadronic showers without charged leptons. Furthermore, this 
signal cannot be confused with ordinary SM events due to an unexpectedly high 
neutrino flux because in the SM 24% of the events are accompanied by muons and 
multiple-bang events are very suppressed, in contrast to the elastic gravitational 
events. The values of the fundamental scale of gravity that IceCube could reach are 
comparable to those to be explored at the LHC. 
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Abstract 

Extremely high energy neutrinos propagating in the atmosphere or in the Earth can 
originate horizontal or up-going air-showers, respectively. We calculate the acceptances 
(event rate/flux) for detecting both types of events by fluorescence detectors, both space- 
based as with the EUSO and OWL proposals, and ground-based, as with Auger, HiRes 
and Telescope Array. We depict them as a function of the neutrino-nucleon cross section, 
crJJ^^, and show that from the ratio of these two classes of events, the inference of 
above 10^^ eV appears feasible, assuming that a neutrino flux exists at these energies. 
Our semi-analytic calculation includes realistic energy-losses for tau leptons and Earth- 
curvature effects. We also consider constraints on shower development and identification 
and the effects of a cloud layer. 



1. Introduction 

Above the Greisen-Zatscpin-Kuzmin (GZK) energy of Eqzk ~ 5 x 10^^ eV [1] 
ultra-high energy neutrinos are probably the only propagating primaries. Moreover, 
in contrast to cosmic-rays, they point back to their astrophysical sources carrying 
information not accessible with other primaries. The detection of ultra- high energy 
neutrinos also allows studies of the fundamental properties of neutrinos themselves, 
as for instance the neutrino-nucleon cross section at energies beyond the reach of 
our terrestrial accelerators. 

In this talk, based on Ref. [2], we study the potential for cosmic-ray experiments 
designed to track ultra-high energy air-showers by monitoring their fluorescence 
yield, to detect horizontal air-showers (HAS) and up-going air-showers (UAS) in- 
duced by a cosmic neutrino flux and show the ability of these experiments to infer 
the neutrino-nucleon cross section, <J^^, at energies above 10^^ eV, from the ratio 
of their UAS and HAS events. Such energies are orders of magnitude beyond the 
energies accessible to man-made terrestrial accelerators. From the point of view of 
QCD, such a cross section measurement would be an interesting microscope into 
the world of small-x parton evolution. Deviations from QCD-motivated extrapola- 
tions [3] could reduce the cross section due to saturation effects [4] or enhance it 
by the existence of new physics thresholds [5] . 

In Ref. [6] it was shown that by comparing the HAS and UAS event rates the 
neutrino-nucleon cross section may be inferred. The calculation of Ref. [6] gave an 
approximate result for the dependence of the UAS event rates on the neutrino- 
nucleon cross section. In this talk, following the results of Ref. [2], we improve upon 
Ref. [6] in several ways, as we show below. On the other hand, the prospects of 
inferring the neutrino-nucleon cross section at neutrino telescopes such as IceCube 
or at the Auger observatory have been studied in Ref. [7] . 
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2. Air-shower rates and constraints on shower-development 

Ultra-high energy neutrinos are expected to arise from the decay of pions and 
subsequently muons produced in astrophysical sources [8] (for the case of produc- 
tion from neutron decays see, eg, Ref. [9]). After propagating for many oscillation 
lengths and due to the maximal mixing between and Vr inferred from terrestrial 
oscillation experiments, all flavors are populated. Thus, a detector optimized for 
or or Vt can expect a measurable flux from cosmic neutrinos. 

The weak nature of the neutrino-nucleon cross section means that HAS begin 
low in the atmosphere, where the target is most dense, and thus that the event rate 
for neutrino- induced HAS is proportional to the cross section. Following Ref [2] , for 
the case of HAS event rates we will only consider Ve charged current interactions. 

For a neutrino- induced UAS, the dependence on the neutrino cross section is more 
complicated. The Earth itself is opaque for neutrinos with energies exceeding about 
a PeV of energy. However, "Earth-skimming" neutrinos, those with a short enough 
chord length through the Earth, will penetrate and exit, or penetrate and interact. 
In particular, there is much interest in the Earth-skimming process Vr t in the 
shallow Earth, followed by r decay in the atmosphere to produce an observable 
shower. In Ref. [6] it was shown that the rate for the Earth-skimming process 
z^T- — + T is inversely proportional to . The inverse dependence of UAS rate on 
(T^fj is broken by the r shower process in the atmosphere. As the cross section 
decreases, the allowed chord length in the Earth increases, and the tau emerges 
with a larger angle from the Earth's tangent plane. This in turn provides a smaller 
path-length in air in which the tau may decay and the resulting shower may evolve. 
This effect somewhat mitigates the inverse dependence of the UAS on . 

The main aim of the study in Ref. [2] was to provide a detailed and improved 
extension of the idea introduced in Ref. [6]. Hence, here we include the energy 
dependences of the tau energy-losses in the Earth, and of the tau lifetime in the 
atmosphere. For the energy-losses, we distinguish between tau propagation in rock 
and in water. In the case of the UAS, the path length of the pre-decayed tau may 
be so long that the Earth's curvature represents a non-negligible correction, that 
we include. We also consider the partial loss of visibility due to cloud layers. On the 
issue of shower development, we incorporate the dependence of atmospheric density 
on altitude and add some conditions for the showers to be observable. Shower 
detection will require that within the field of view, the length of the shower track 
projected on a plane tangent to the Earth's surface exceeds some minimum length, 
/,nin- In addition, a minimum column density, dmin, beyond the point of shower 
initiation is required for the shower to develop in brightness. On the other hand, 
after a maximum column density, dmax, the shower particles are below threshold 
for further excitation of the N2 molecules which provide the observable fluorescence 
signal. Therefore, visible showers end at dmax- Finally, the fluorescent emission 
per unit length of the shower will decline exponentially with the air density at 
altitude. We will take zthin = 24 km as the altitude beyond which the signal becomes 
imperceptible. Regarding the choice of dmin and dmax, they are inferred from the 
observed longitudinal development profiles of ultrahigh-energy cosmic ray showers. 
On the other hand, we assign a relatively small value to ^min to maximize the 
observable event rate. For a summary of the different values adopted to obtain the 
results, full details on the analytic description of the effects of these parameters on 
the event rates and comparison with prior work, we defer the reader to Ref. [2]. 
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Fig. 1. Acceptances with l^i^^ fixed at 5 km, rfmin at 400g/cm^, and dmax at 1200g/cm^. The 
curves correspond to HAS (dotted hne) and UAS over ocean with E„ = 10^^ eV (thick sohd hne), 
ocean with Ei, = 10^^ eV (thin solid line), land with Ei, = 10^^ eV (thick dashed line), and land 
with El/ = 10^" eV (thin dashed line). Panels are for (a) space- based (or ground-based) detectors 
in the absence of clouds with E^^ = 10^® eV; (b) ground-based detectors in the presence of a 
cloud layer at ^cloud = 2 km with E^^ = 10^® eV; (c) spaced-based detectors in the presence of 
a cloud layer at z^ioud = 2 km with E^^ = 10^® eV; (d) space-based (or ground-based) detectors 
in the absence of clouds with E^^ = 5 X 10^^ eV; (e) ground-based detectors in the presence of a 
cloud layer at ^cloud = 2 km with E^^ = 5 X 10^^ eV; (f) spaced-based detectors in the presence 
of a cloud layer at ^cloud = 2 km with E^^ = 5 X lO""^^ eV. 

3. Results 

In this section, we present the rcsuhs of our semi-analytical approach and take 
the product of area and solid angle ^ 10^ km^ sr, i. e. that of the EUSO design 
report [10]. 

In Fig. 1 are plotted UAS (solid and dashed) and HAS (dotted) acceptances in 
units of (km^-sr), versus fixed values of cr^^. Within the approximations followed in 
Ref. [2] , for the ideal case of a cloudless sky (panels a and d) there is no difference 
between the acceptances for ground-based and space-based detectors. However, 
there are significant up-down differences when the sky is covered by clouds (panels 
b, c, e and f). In this latter case, we model the cloud layer as infinitely thin with 
altitude Zdoud, but with an infinite optical depth so that showers are completely 
hidden on the far side of the cloud layer. 

The HAS acceptances depend on neutrino energy only via a^^{Ei/), and rise 
linearly with crj?^. Plotted against fixed cr^^, then, the straight-line HAS curves 
(dotted) are universal curves valid for any i?^ exceeding the trigger threshold E^^. 
The UAS acceptances have a complicated dependence on E^; it arises from the 
energy dependences of v propagation in the Earth, tau propagation in the Earth, 
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and path-length of the tau in the atmosphere before it decays, the latter also affect- 
ing the visible shower characteristics. We can clearly see that the UAS acceptance 
(and so also the rate) is typically an order of magnitude larger when neutrinos 
traverse a layer of ocean water, compared to a trajectory where they only cross 
rock. Thus, the UAS event rate is enhanced over the ocean relative to over land. 
The value of this enhancement depends on the shower threshold-energy E^^ of the 
detector (upper versus lower panels) and on the neutrino-nucleon cross section in a 
non-trivial way. On the other hand, quantitatively, the ground-based acceptances 
are quite reduced by the low-lying clouds, whereas the space-based acceptances are 
not, as one would expect. The suppression of the ground-based acceptance is most 
severe for small cross sections, for which the tau leptons emerge more vertically 
and disappear into the clouds before their eventual shower occurs and develops. 
Ground-based UAS acceptances are reduced by up to an order of magnitude over 
water, and even more over land. Ground-based HAS acceptances are reduced by 
an order of magnitude. For space-based detectors, the UAS acceptance is reduced 
little by clouds at 2 km. Larger neutrino cross sections lead to more tangential tau- 
showers which may hide below a low-lying cloud layer. We see that UAS reductions 
are a factor of 2 for the larger cross sections shown, and less for the smaller values 
of cross section. 

We obtain benchmark event rates by multiplying our calculated acceptances with 
a benchmark integrated flux of one neutrino per (km^ sryr). The result is a signal 
exceeding an event per year for an acceptance exceeding a (km^-sr). Thus we see 
that this benchmark flux gives a HAS rate exceeding 1/yr if cr^^ exceeds 10"'^^ cm^; 
and an UAS rate exceeding 1 /yr over water for the whole cross section range with 
E^ll ~ 10^^ eV, and over land if a^^ < 10~'^^cm^. When E'^^ is raised, however, 
the UAS signal over land is seriously compromised, while UAS rates over the ocean 
are little changed, HAS rates are unchanged, as long as E^y^ exceeds E^. 

We call attention to the fact that for UAS over both ocean and land, there 
is a maximum in the UAS acceptance at cross section values ^ (1 — 2) x 
10~^^ cm^ and ~ (0.3 — 0.5) x 10"'^^ cm^, respectively. For cross sections 
similar or smaller than those at the maximum, the acceptance for UAS is larger 
than that for HAS; conversely, for cross sections above those at the maximum, HAS 
events will dominate UAS events. The cross section value at the maximum lies just 
below the extrapolation of the Standard Model cross section, which for the two 
initial neutrino energies considered, 10^" eV and 10^^ eV, is 0.54 x 10~'^^cm^ and 
1.2 X 10~'^^cm^, respectively. If this extrapolation is valid, then one would expect 
comparable acceptances (and event rates) for UAS over water and for HAS. If the 
true cross section exceeds the extrapolation, then HAS events will dominate UAS 
events; if the true cross section is suppressed compared to the extrapolation, then 
UAS events will dominate HAS events. Importantly, the very different dependences 
on the cross section of the HAS and UAS acceptances offers a practical method 
to measure <t^^. One has simply to exploit the ratio of UAS-to-HAS event rates. 
Furthermore, the shape of the UAS acceptance with respect to establishes the 
"no- lose theorem" [2,6], which states that although a large cross section is desirable 
to enhance the HAS rate, a smaller cross section still provides a robust event sample 
due to the contribution of UAS. The latter is especially true over ocean. 
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4. Conclusions 

In this talk wc have presented a mostly analytic calculation of the acceptances 
of space-based and ground-based fluorescence detectors of air-showers at ultra- 
high energies. Included in the calculation arc the dependences of the acceptances 
on initial neutrino energy, trigger-threshold for the shower energy, composition of 
Earth (surface rock or ocean water), and several shower parameters (the minimum 
and maximum column densities for shower visibility, and the tangent length of the 
shower). Also included in the calculation are suppression of the acceptances by 
cloud layers and by the Earth's curvature. And most importantly, also included are 
the dependences on the unknown neutrino-nuclcon cross section. The dependence 
is trivial and linear for HAS, but nontrivial and nonlinear for UAS. 

The merits of the analytic construction are two-fold: it offers an intuitive un- 
derstanding of each ingredient entering the calculation; and it allows one to easily 
re-compute when different parameters are varied. While a Monte Carlo approach 
may be simpler to implement, it sacrifices some insight and efficiency. 

The differing dependences of HAS and UAS on enable two very positive 
conclusions: (1) the "no-lose theorem" is valid, i. e. that acceptances are robust for 
the combined HAS plus UAS signal regardless of the cross section value; (2) and 
an inference of a^^ above 10^^ eV is possible if HAS and UAS are both measured. 
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Abstract 

The spectrum of extra-galactic cosmic rays (CRs) is expected to follow the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff at about 5 x 10^*^ GeV which results from energy losses of 
charged nuclei in the cosmic microwave background. So far the confrontation of this feature 
with CR data is inconclusive. In the absence of close-by sources a power-law continuation 
of the spectrum might signal the contribution of new physics. We have investigated the 
statistical significance of a model where exotic interactions of cosmogenic neutrinos are 
the origin of super-GZK events. A strong neutrino-nucleon interaction is favored by CR 
data, even if we account for a systematic shift in energy calibration. 



1. Introduction 

The appearance of extremely high energetic (EHE) cosmic rays (CRs) is a mys- 
tery. Simple geometric arguments show that the energy of nuclei originating in 
cosmic accelerators should be limited corresponding to the size and magnetic field of 
the accelerating environment (Hillas criterion) [1]. Additional constraints arise from 
energy losses in the source that has to be balanced by the acceleration rate [2] . In 
particular, the energy loss length of synchrotron radiation decreases with the third 
power of the magnetic field (at constant gyro-radius) and limits the efficiency of 
small candidate acceleration sites. The very few sources that seem to be capable of 
accelerating protons up to 10^^ GeV include radio galaxy lobes and relativistically 
moving sources like jets of active galactic nuclei or gamma ray bursts. 

Not only the acceleration mechanism of these particles seems to be problem- 
atic, but also their propagation in the interstellar medium. Resonant pion photo- 
production in the cosmic microwave background (CMB) limits the range of EHE 
particles to a few 10 Mpc. It was first emphasized by Greisen [3], Zatsepin, and 
Kuzmin [4] (GZK) that this should be signaled by a cutoff in the CR spectrum at 
about 5 X 10^° GeV if extra-galactic sources dominate at these energies. So far, the 
confrontation of this effective cutoff with CR data is inconclusive. In particular, the 
observations of two large exposure experiments AGASA [5] and HiRes [6] seem to 
imply conflicting results. 

The "excess" of super-GZK events reported by the AGASA collaboration has 
led to speculations about a different origin of EHE CRs. Berezinsky and Zatsepin 
proposed that cosmogenic neutrinos [7] produced in the decay of the GZK pions 
could explain these events assuming a rapid rise of the neutrino-nucleon interaction. 
We have followed this idea in Ref. [8] and investigated the statistical significance of 
scenarios with strongly interacting neutrinos combining CR data from AGASA and 
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Fig. 1. Left panel: The cosmic ray spectra from AGASA, Fly's Eye Stereo, Hires-I/II (with energy 
shift given in Ref. [11]), and PAO. Right panel: The flux of extra-galactic protons and cosmogenic 
neutrinos from CMB interactions corresponding to the 99% CL of the fit. 

HiRes and neutrino limits from horizontal events at AGASA and contained events 
at RICE. 

The AGASA excess of super-GZK events could also be the result of a relative 
systematic error in energy calibration of ±30% [9]. This interpretation is fueled by 
a recent re-analysis of the AGASA data [10]. The preliminary results indicate that 
the energy of highly inclined showers has previously been over-estimated. If verified, 
this would result in a re-calibration of the spectrum according to a relative energy 
shift of -10% to -15%. 

Remarkably, this energy shift agrees with predictions from a matching procedure 
of CR spectra assuming an early onset in extra-galactic proton dominance at about 
10® GeV [11]. The calibration by the "dip" in these proton spectra resulting from 
electron-positron pair production on the CMB predicts a —10% and 20% energy 
shift of the AGASA and HiRes data, respectively. We have incorporated this re- 
calibration into our analysis and show also the effect of the preliminary data from 
the Pierre Auger Observatory (PAO) [12] on our results. 



2. Cosmic Sources of Protons and Neutrinos 

We assume that extra-galactic protons dominate the CR spectrum above 5 x 
10* GeV. A convenient parameterization of a spatially homogeneous and isotropic 
source luminosity is given by 

CcRiz,E) (X {1+ z)" E^'' e^l^, ^mln < ^ < ^max, (l) 

which accounts for an evolution with redshift z and an exclusion of nearby ( Zmin) 
and early (zmax) sources. In our analysis we have keep these parameters fixed at 
-Zmin = 0.012, corresponding to r,„in « 50Mpc, and Zmax = 2.0. The power-law 
injection is effectively limited to energies below -Emax, which we also fix at 10^^ GeV. 

The flux of protons originating at distant sources is subject to energy redshift and 
collisions with the interstellar photon background during propagation. The dom- 
inating interactions of UHE protons arc e^e^ pair production and meson photo- 
production in the CMB (see e.g. [13]). These effects are taken into account by means 
of propagation functions, which have been provided by the authors of Ref. [14]. The 
pions, which are resonantly produced at the GZK cut-off, decay into electron and 
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muon neutrinos via the reaction chain 7r+ fJ-^i^^i ^ i^ni^i^ii^e & ^ and the conjugate 
process. This constitutes a flux of UHE eosmogenic neutrinos which can be taken as 
a "guaranteed source" assuming extra-galactic protons as the highest energy CRs. 

Depending on the environment of the proton accelerators neutrinos might also 
emerge in the decay of pions produced by photo-production processes during accel- 
eration [15]. The flux of extra-galactic neutrinos have been estimated in Ref. [16] 
assuming a low cross-over to extra-galactic protons at about 5 x 10^ GeV. In con- 
trast to a transition at the "ankle" this flux of neutrinos from the source dominate 
the total neutrino flux below 10*^ GeV and exhausts current experimental neutrino 
bounds at energies of about 10^ GeV. The total neutrino flux beyond the GZK 
cutoff is however dominated by eosmogenic neutrinos and we will only use this 
contribution for our statistical evaluation. 



3. Exotic Neutrino Interactions 

The flux of eosmogenic neutrinos associated with a low cross-over to extra-galactic 
protons is comparable to the flux of protons at the GZK cut-off as can be seen from 
Fig. 1. A contribution of these neutrinos as super-GZK events requires a strong 
deviation from the feeble neutrino-nueleon interactions predicted by the Standard 
Model. In order to avoid large contributions in neutrino experiments the transition 
to a anomalously large cross section has to be very rapid. 

The realization of such a behavior has been proposed in scenarios beyond the 
(perturbative) Standard Model (SM), e.g. arising through compositeness, through 
eleetroweak sphalerons, through string excitations in theories with a low string and 
unification scale, through Kaluza-Klein modes from compactified extra dimensions, 
or through black hole and p-brane production, respectively. For details of these 
models we refer to some recent reviews, Refs. [17]. Skepticism about these scenarios 
has been raised in Ref. [18] based on limitations from s-wave unitarity or the natu- 
ralness of the couplings. However, these considerations do not necessarily apply for 
non-probative aspects like instant-ons. 

In the following we will use a flexible parameterization of a strong neutrino- 
nueleon inelastic cross section (cr"™) focusing on three characteristic parameters: 
(i) the energy scale E'th of the new underlying physics, (ii) the amplification A 
compared to the SM predictions, and (iii) the width B of the transition between 
weak and strong interaction. A mathematical convenient parameterization is given 

by 



logi 



•jN 



1 



1 + tanh ( logg ' 



(2) 



In general, experiments distinguish different CR primaries by their characteristic 
shower development in matter. For the sake of simplicity we will assume that the 
characteristics of the showers induced by strongly interacting neutrinos are indis- 
tinguishable from those induced by protons. In particular, we assume for both 
primaries (i) a complete conversion of the incident energy into the shower, and (ii) 
equal detection efheiencies at the highest energies. 



4. Quantitative Analysis 

Following the procedure of Ref. [8] we have evaluated the CR spectra of AGASA, 
Fly's Eye Stereo, HiRes-I/II, and PAO (preliminary data) assuming extra-galactic 
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Fig. 2. The range of the cross section within the 95%, 90% and 68% CL (outer to inner bands) ob- 
tained by the re-calibrated AGASA, Fly's Eye Stereo, and HiRes-I/II data, with (right panel) and 
without (left panel) the preliminary Auger data. Also shown are theoretical predictions of the neu- 
trino-nucleon cross section enhanced by electroweak sphalerons, p-branes, and string excitations 
(see Ref. [8] for details). 



proton dominance above and 5 x 10* GeV and a contribution of strongly interacting 
neutrinos as supcr-GZK events. In contrast to Ref. [8] we have re-calibrated the 
energy by a multiplicative factor 0.9 and 1.2 for the AGASA and HiRes data, 
respectively, following Ref. [11]. For consistency with neutrino experiments we have 
used the search results of horizontal events at AGASA and contained events at 
RICE [19]. Figure 2 shows the range of the cross section corresponding to the 68%, 
90%, and 95% confidence level (CL) from a goodness-of-fit test assuming cosmogenic 
neutrino fluxes from a low crossover scenario. The details of the statistical analysis 
and the approximations involved can be found in Ref. [8]. 

Compared to our previous results [8] the re-calibrated spectrum is in a much 
better statistical agreement with the combined flux of protons and neutrinos from 
exotic interactions. This is mainly due to the fact that the calibration method 
"optimizes" the CR data to extra-galactic proton spectra below the GZK cutoff. 
However, within this particular model the re-calibrated data still favors an addi- 
tional contribution from strongly interacting cosmogenic neutrinos. At the 95% CL 
the corresponding neutrino-nucleon cross section should exhibit a steep increase by 
an amplification factor of ^ > 10^. For A < 10^ the transition should be very rapid 
(B < 10) at about 5 x 10^^ GeV. The allowed parameter space further shrinks if 
we also include the preliminary PAO data. 

As an illustration of the result, we have considered three models of a rapidly 
increasing neutrino-nucleon cross section based on electroweak sphalerons [20], p- 
branes [21] and string excitations [22]. The details can be found in Ref. [8]. A 
separate fit of the source luminosity (Eq. (1)) with these cross section using the re- 
calibrated data gives a statistical acceptance at the 81%, 65%, and 60% CL (without 
PAO), respectively. If we also include the PAO data in the fit, these values change 
to 98%, 90%, and 88%, respectively, which is in reasonable good agreement with 
the results of Fig. 2. 
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5. Conclusions 

Wc have shown that in the absence of close-by sources exotic interactions of 
cosmogenic neutrinos might extend the cosmic ray spectrum beyond the GZK cut- 
off. Such a behavior is predicted in various extensions of the Standard Model, 
but might also be due to non-pcrturbative aspects like electro-weak sphalerons. In 
extension to our analysis of Ref. [8] we have presented here the results of a fit to 
the AG AS A and HiRes data, shifted according to the "dip" -calibration from extra- 
galactic proton spectra [11]. With these modifications a contribution of neutrinos 
is still favored by the data. The inclusion of preliminary data of the Pierre Auger 
Observatory weakens the results of the goodness-of-fit test. 
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